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Abstract. 

We use the largest complete sample of 64 galaxy clusters (Highest X-ray FLUx Galaxy Cluster Sample) with available high- 
quality X-ray data from Chandra, and apply 16 cool-core diagnostics to them, some of them new. In order to identify the best 
parameter for characterizing cool-core clusters and quantify its relation to other parameters, we mainly use very high spatial 
resolution profiles of central gas density and temperature, and quantities derived from them. We also correlate optical properties 
of brightest cluster galaxies (BCGs) with X-ray properties. 

To segregate cool core and non-cool-core clusters, we find that central cooling time, r coo i, is the best parameter for low 
redshift clusters with high quality data, and that cuspiness is the best parameter for high redshift clusters. 72% of clusters in our 
sample have a cool core (r coo i < 7.7 hZ\ 12 Gyr) and 44% have strong cool cores (f coo i < 1.0 hZ[ 12 Gyr). We find strong cool-core 
clusters are characterized as having low central entropy and a systematic central temperature drop. Weak cool-core clusters have 
enhanced central entropy and temperature profile that are flat or decrease slightly towards the center. Non-cool-core clusters 
have high central entropies. 

For the first time we show quantitatively that the discrepancy in classical and spectroscopic mass deposition rates can not 
be explained with a recent formation of the cool cores, demonstrating the need for a heating mechanism to explain the cooling 
flow problem. 

We find that strong cool-core clusters have a distribution of central temperature drops, centered on 0.47\,j r . However, the 
radius at which the temperature begins to drop varies. This lack of a universal inner temperature profile probably reflects the 
complex physics in cluster cores not directly related to the cluster as a whole. Our results suggest that the central temperature 
does not correlate with the mass of the BCGs and weakly correlates with the expected radiative cooling only for strong cool- 
core clusters. Since 88% of the clusters in our sample have a BCG within a projected distance of 50 kpc from the X-ray 
peak, we argue that it is easier to heat the gas (e.g. with mergers or non-gravitational processes) than to separate the dense core 
from the brightest cluster galaxy. 

Diffuse, Mpc-scale radio emission, believed to be associated with major mergers, has not been unambiguously detected 
in any of the strong cool-core clusters in our sample. Of the weak cool-core clusters and non-cool-core clusters most of the 
clusters (seven out of eight) that have diffuse, Mpc-scale radio emission have a large (>50 Kj\ kpc) projected separation between 
their BCG and X-ray peak. In contrast, only two of the 56 clusters with a small separation between the BCG and X-ray peak 
(<50 \C,\ kpc) show large-scale radio emission. Based on result, we argue that a large projected separation between the BCG 
and the X-ray peak is a good indicator of a major merger. The properties of weak cool-core clusters as an intermediate class of 
objects are discussed. Finally we describe individual properties of all 64 clusters in the sample. 



1. Introduction 

Early X-ray observations of galaxy clusters revealed that the 
intracluster medium (ICM) in the centers of many clusters was 
so dense that the cooli ng time of the g a s was much shorter than 
the Hubble time (e.g iLea et alJll973l: ICowie & Binnevl[l977 ; 
Fabian & NulseiJl977llMathews & Bregmanll978l) . These ob- 
servations led to the development of the cooling flow (CF) 
model. In this model the ICM at the centers of clusters with 



dense cores hydrostatically cools, so that the cool gas is com- 
pressed by the weight of the overlying gas. Hot gas from 
the outer regions of the ICM flows in to replace the com- 
pressed gas, generating a CF. Although early X-ray observa- 
tions seemed to corroborate this model and there was some 
evidence of expected Ha emission and UV emission, opti- 
cal observations failed to det ect the expected star formation 



rates, CO a nd molecular gas dMcNamara & 0'Connelllll989 



Edge! I2001I) . Later, discrepancies between spectrally deter- 
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mined mass deposition rates from ASCA and the classical deter- 
mination fro m density further brough t the classical CF model 



More recently, gratin 


g spectra from XMM (Peterson et al. 


2001 


Tamura et al. 


2001b 


a: KaastraetalJ 2001; Xuetal. 


2002 


Sakelliou et al.ll2002; Peterson et al.ll2003; Sanders et al. 


2008 



have revealed that the gas in central regions is not cooling 
at the rates predicted by the traditional cooling flow model. 
The classical CF model assumes that no mechanism exists 
to significantly heat the gas. Therefore, this discrepancy has 
inspired a search for heating models that can explain cur- 
rent observations . Among those mechanisms proposed are: 



conduction (e.g. IZakamska & Naravanl 120031) . central AGN 



heating via direc t cosmic ray-ICM interaction+conduction 



dGuo & Ohl 120081). AGN heating by bubble induced weak 



shoc ks (e.g. [Ma thews et al 



Dennis & Chandran 



cent studies (e.g. 



2006), soundwaves+conduction 



(e.g. Ruszkowski et al. 2004) and turbulence+conduction (e.g. 



20051). On the observational side 
20041 



Birzanetal 



Raff erty et al.l 120061 



Dunn & Fabianll2006[l2008[ lBirzan et alll2008l) have looked at 
the mechanical energy released by AGN driven bubbles and 
compared it to the energy needed to quench the CF. 

The failure of the classical CF model has changed 
the nomenclature of these centrally dense clusters from 
CF clusters to cool-core (CC) clusters as suggested by 



Mol endi & Pizzolatol (120011) . One problem with this nomencla- 
ture is that it is unclear what distinguishes a CC cluster from a 
non-CC (NCC) cluster. The name implies that gas in the center 
of the cluster is cool, but does that always imply a short cool- 
ing time? In fact authors define CC c lusters differently often 
based on central t emperature drop (e.g lSanderson et aDl2006al; 
Burns et al. 2008b. short central cooling ti me (e.g. lBauer et al.l 
20051: lO'Hara et al.ll2006l: lDonahuell2007h. or significant clas- 
sical mass deposition rate (IChen et al.l 120071) . There also is a 
question as to whether there is a distinct difference between 
NCC and CC clusters. That is, is there a parameter that un- 
ambiguously distinguishes NCC clusters from CC clusters? 
This is a nontrivial question since, when used as cosmologi- 
cal probes, clusters are often segregated into CC/NCC subsam- 
ples. Frequently CC clusters are chosen for mass determina- 
tion studies since they are believed to be dynamically relaxed. 
On the other hand, it requires a large amount of energy to 
quench a CF, which may stro ngly affect the entire ICM. For 
example, O'Hara et al. ( 20061) found that CC clusters (as de- 
fined as having a central cooling time (CCT) more than 3 cr 

— 1/2 

below 7.1 h 1Q ' Gyr) have more scatter about scaling relations 
than NCC clusters. Therefore it becomes important to unam- 
biguously differentiate between CC and NCC clusters before 
proceeding with determining their effects on scaling relations. 

It is worth noting that a discrepancy between mass deposi- 
tion rates determined spectroscopically and those from images 
in some clusters does not indicate in itself any breakdown of 
the simple cooling flow picture, since it may be that these clus- 
ters just did not have enough time to deposit the predicted level 
of mass. Here, we will show for the first time that this recent 
formation solution to the cooling flow discrepancy is very un- 
likely, based on statistical arguments. 



Recently, IChen et al.l (12007b investigated the cores of a 
sample of 106 galaxy clusters based on the extended Highest 
X-ray FLUx Galaxy Cluste r Sample (HIFLUGCS) sample 
dReiprich & Bohringerl Eo02h using ROSAT and ASCA, with 
their primary goal to study the effect of CC versus NCC clusters 
on scaling relations. They defined a CC cluster as one that has 
significant classically defined mass deposition rate (M c iassicai) 
and found 49% of their sample were CC clusters. Additionally 
they determined that the fraction of CC clusters in t heir sam ple 



d2006al) 



decreased with increasing mass. ISanderson et al. 
vestigated the cores of a statistically selected sample of the 20 
brightest HIFLUGCS clusters using Chandra. They find that 
nine (~41 %) of their clusters are CC clusters, where they define 
CC clusters as having a significant drop in temperature in the 
central region. Additionally they find that the slope of the inner 
temperature profiles has a bimodal distribution with NCC clus- 
ters having a flat temperature pro file and CC cluster s having 
a slope (dlog(7/)/dlog(r)) of 0.4. iBurns etall d2008l) investi- 
gated the properties of CC and NCC clusters from a cosmo- 
logic al simulation wh i ch pro duced both types of clusters. As 
with Sanderso n et al. I d2006al) . they define CC clusters based 
on the central temperature decrease (by 20% of T v ; r in their 
case). Their simulations show that an early merger determines 
whether a cluster is a CC cluster or not and suggest that CC 
clusters are not necessarily more relaxed than NCC clusters. 

Additionally, studies have begun on samples of distant 
(high-z) clusters. The goal of these studies is to determine the 
physical evolution of CC clusters and whether the CC fraction 
changes with redshift. The major obstacle for these high-z clus- 
ter studies is that the traditional indications of CCs are difficult 
to measure for distant clusters. Therefore authors offer prox- 
ies, based on studies of low-z clust ers, that can be measure d for 
distant clusters with limited signal. IVikhlinin et al.l d2007l) sug- 
gested cuspiness as a proxy for determining whether a distant 
cluster is a CC cluster or not. Cuspiness is defined as: 



dlog(n) 
dlog(r) 



0.04r 5l 



(HI- 



(1) 



where n is the electron density and r is the distance from the 
cluster center. For their sample of 20 clusters with z > 0.5, they 
find a lack of CC clusters. That is, they find no strong CC clus- 
ters (a > 0.75 by their definition) and few weak CC clusters 
(a = 0.5-0.75). Approximately 67% of their nearby sample (a 
subsample of 48 HIFLUGCS clusters) has a > 0.5 versus 15% 
for their z > 0.5 sa mple. On the other hand, using simulations, 
Burns et al.l (120081) find that th e CC fraction stays c onstant back 
to a redshift of z ~T. Recently lSantos et al.l (120081) compared a 
nearby sample of 1 1 clusters (z = 0.15 - 0.3) to a distant sam- 
ple of 15 clusters (z - 0.7 - 1.4). They suggest a concentration 
parameter as a proxy for identifying CC clusters: 



Csb = 



E(r < 40 kpc) 
Z(r < 400 kpc) ' 



(2) 



where E(r < 40 kpc) and Z(r < 400 kpc) are the integrated 
surface brightnesses within 40 kpc and 400 kpc respectively. 



1 The Coma, Fornax, and Centaurus clusters were excluded because 
of their large angular size. 
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As with IVikhlinin et al.l (120071) . they find no strong CC clus- 
ters in their high redshift sample. On the other hand they find 
roughly the same fraction of weak (or moderate) CC clusters 
in the distant as in th e local sample, in apparent conflict with 
Vikhlinin et alld2007l). However, s ince both the nearby and dis- 
tant samples o flSantos et al.l d2008l) are not complete, their frac- 
tions may be biased by selection factors, so a straightforward 
comparison is not possible. 

In this paper we investigate a statistically complete, flux- 
limited sample of 64 X-ray selected clusters and analyze their 
cores with the Chandra ACIS instrument. This is the first de- 
tailed core analysis of a nearby, large, complete sample with 
a high resolution instrument. With its ~075 point spread func- 
tion, Chandra is the ideal instrument for such a study. Likewise, 
the HIFLUGCS sample, which compris es the 64 X-ray bright- 
est cl usters outside the Galactic plane dReiprich & Bohringei 
2002b . is an ideal sample. Since the clusters are bright, they 
are nearby ((z) = 0.053) making it possible to probe the very 
central regions of the clusters. Moreover they have high sig- 
nal to noise allowing precise measurements to be made. It is 
worth noting that complete flux-limited samples are not nec- 
essar ily unbiased or representative with respect to morphology 
(e.g. Reiprich 120061) . The point is that, for flux-limited sam- 
ples, the bias can, in principle, be calculated (e.g., Ikebe et al. 
2002, Stanek et al. 2006). Also, such samples are representa- 
tive in the sense that their statistical properties, like cooling 
core frequency, are directly comparable to the same properties 
of simulated flux-limited samples. Samples constructed based 
on availability of data in public archives are, in general, not 
representative in this sense ("archive bias"). 

Our goal in this paper is to determine if there is a physi- 
cal property that can unambiguously segregate CC from NCC 
clusters. We use this property to examine other parameters to 
see how well they may be used as proxies when determining 
whether a cluster is a CC cluster or not. The article is organized 
as follows. We outline our methods of data reduction in Sect. [2] 
We present our results in Sect. [3] wherein we describe the var- 
ious parameters investigated for determining whether a cluster 
is a CC or an NCC cluster in Sect. 13.11 determine the best pa- 
rameter to segregate distant CC and NCC clusters in Sect. 13.21 
compare this parameter to others in Sect. l3.3l and determine the 
best diagnostic for cool cores in distant clusters in Sect. 13.41 
We discuss our results in Sect. [4] wherein we describe the basic 
cooling flow problem in Sect. 14. H and discuss the inner temper- 
ature profiles in Sect. 14.21 the central temperature drop seen in 
some clusters in Sect. 14.21 the relation between cluster mergers 
and the projected separation between the BCG and X-ray peak 
in Sect. I4.3l and the WCC clusters in Sect. 14.41 Finally we give 
our conclusions in Sect. [5] To aid clarity, we give in Table Q] 
some of the abbreviations used throughout the paper. 

In this work, we assume a flat ACDM Universe with Qm 
= 0.3, O vac = 0.7, and H Q = h lx 71 km s -1 Mpc" 1 . Unless 
otherwise noted, k is the Boltzmann constant and we use the 
following nomenclature for our coordinates: x is the projected 
distance from the cluster center, r is the physical 3-D distance 
from the cluster center and I is the distance along the line of 
sight. All errors are quoted at the 1 <x level unless otherwise 
noted. 



Table 1. Nomenclature 



Abbreviation 


ICM 


IntraCluster Medium 


CCT 


Central Cooling Time at 0.4%R 5 qq 


sec 


Strong Cool Core 


WCC 


Weak Cool Core 


NCC 


Non Cool Core 


CC 


Cool Core (SCCs+WCCs) 


CF 


Cooling Flow 


BCG 


Brightest Cluster Galaxy 


EP 


Emission Peak (X-ray) 


BCG-EP 


Projected separation between BCG and EP 



2. Observations and Methods 

All 64 HIFLUGCS clusters have been observed with Chandra, 
representing 4.552 Ms of cleaned data. For our analysis we 
used all unflared data taken after the CCD focal plane tem- 
perature was reduced to -120 °C (2001-Jan-29) and was pub- 
licly available as of 2007-May-07, with a few exceptions. (1) In 
cases where more than one data set was publicly available and 
one set was heavily flared, the entire flared data set was dis- 
carded, (2) for A7541 and A4O10, (3) in the cases of A2597 
and A2589 the original observations were heavily flared and 
the Pis (T. Clarke and D. Buote, respectively) of the newer pro- 
prietary data sets graciously provided them to us before they 
became publicly available (4) the Coma Cluster (A1656), has a 
recently released calibration observation (2008-Mar-23) which 
we used instead of the pre-2001-Jan-29 observations. Several 
other Pi's kindly provided data sets before they were publicly 
available, although they since have become publicly available 
(see acknowledgments). Details on the number and length of 
the observations can be found in Table Q] Note that we use 
redshifts compiled from NASA/IPAC Extragalactic Database 
(NEDf] and the values did not differ significantly (Az < 0.001) 
from the values used bv lReiprich & Bohringerld2002l) . 



2.1. Data Reduction 

The basic data reduction was done using CI AO 3.2.2 and 



CALD B 3.0 following the methods outlined in Hudson et al 
d2006h . Since we are interested in the cores of clusters, we 
took the cluster center to be the emission peak (EP) from 
the background sub tracted, exposure corrected image. See 
iHudson et alj d2006l) for details of our image creation tech- 
nique To avoid statistical fluctuations in determining the clus- 
ter peak, we smoothed all images with a Gaussian. In general 
we used a ~4" kernel, however in a few particularly noisy 



- There are later observations of A754, but they are all far offset and 
not useful for core studies and therefore not included. 

3 There is a post-200 l-Jan-29 observation of A401, however it is 
offset such that the cluster center is in the corner of a CCD. We there- 
fore included the centered pre-2001-Jan-29 observation as well. 

4 http://nedwww.ipac.caltech.edu 

5 For this analysis we did not apply adaptive smoothing. 
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used the radio measured hydrogen abs orption columns from 



the Leiden/ Argentine/Bonn Hi-Survey dKalberla et al.l I2005) 



13:50:00 



13:49:00 
a(J2000) 



13:48:00 



Fig. 1. An example of a mosaiced image. This is the back- 
ground subtracted, exposure corrected image created from the 
12 Chandra exposures of A 1795. 

cases (noted in Appendix we used a larger smoothing ker- 
nel. Fig.[T]shows an example of a mosaiced image. 

2.2. Temperature Profiles 

Using the EP as the cluster center we created a temperature 
profile for each cluster using the following method. In order to 
get good statistics for our temperature profiles, we estimated 
that 10000 source counts per annulus would give us a ~5% 
constraint on T. However, in several cases we had a small 
number of total source counts (~30000), so that using 10000 
source counts per annulus was not practical. We therefore ex- 
tracted annuli with at least 5000 source counts for any cluster 
with fewer than 100000 total source counts and 10000 source 
counts for any cluster with more than 100000 total source 
counts. We created pie shaped annuli, not extending beyond the 
boundary of the chips. Although this decreased our efficiency, it 
allows for unbiased density measurements when using the nor- 
malization of the thermal model (a full circle extending beyond 
the chip edge would generally cover more of the inner part of 
the annulus than the outer). In our subsequent studies it will 
also allow us to use the same spectra for deprojection mod- 
els as well as helping to constrain the residual cosmic X-ray 
background (CXB) in the outer cluster regions^. We then cre- 
ated spectra and fit the spectrum to an absorbed thermal model 
(WABS*APEC*EDGE3). For the majority of observations we 



6 The normalizations for the cluster emission for partial annuli of 
different observations can be scaled and tied when trying to determine 
the residual CXB. 

7 The edge component is a c orrection to underes timation 
Chandra's efficiency at ~2 keV (see lVikhlinin et alj2005l) . 



or, when available, pointed observations from E. Murphy (pri- 
vate communication). In a few cases, noted in Appendix|C] the 
X-ray measured value was significantly higher than the radio 
measured value and so we left the hydrogen absorption colu mn 
as a free parameter. See lBaumgartner & Mushotzkv d2006l) for 
an explanation of the discrepancy between radio and X-ray ab- 
sorption column density measurements. As noted in Table [2] 
and Appendix [C] there were 14 cases (~22%) in which we 
added a second thermal component to one or more of the cen- 
tral annuli in order to improve the fit. In these cases, the lower 
of the two temperatures was used for our annular analysis. 
There were also a few cases in which using non-solar metalicity 
ratios significantly improved the fit, however since the best-fit 
temperature was unaffected, solar ratios were used for simplic- 
ity. See Appendix ICl for notes on the clusters where this was 
the case. 



2.3. Cluster Virial Temperature, Radius and Mass 

The cluster virial temperature (r v i r ) is the temperature of the 
X-ray emitting gas that is in hydrostatic equilibrium with the 
cluster potential. Although T vir is not an imperative quantity 
for core studies, it is useful for scaling clusters for compari- 
son purposes (e.g. relative temperature drop) or estimating the 
cluster virial radius and mass. In this section we describe our 
method for determining T wn and how we use it to estimate the 
virial mass and radius. 

The temperature decline in the central regions of some clus- 
ters, if included, is the largest source of bias in the determi- 
nation of r v j r . The reason for this temperature drop is a cur- 
rent topic of debate and discussed in more detail in Sect. 14.21 
Here we simply discuss the removal of this central region, so 
it does not bias the fit to J v i r . In order to determine the size 
of the central region to be excluded, we fit our temperature 
profiles to a broken powerlaw. The core radius in the power- 
law was free and the index of the outer component was fixed 
to be zero. Additionally we removed the outer annuli where 
accurate subtraction of the blank sky backgrounds (BSB) be- 
comes critical a nd where clusters may have a decreasing tem- 
perature profile ( Ma rkevitch et al. 1998c De Grandi & Molendi 



of 



2002l;IVikhlinin et al.ll2005HBurns et aT]|2008l) . The core radius 
in the powerlaw was taken to be the radius of the excluded cen- 
tral region. Table Q] gives this core radius as well as kT wlI and 
overall cluster metalicity. Fig.|2]shows examples of the broken 
powerlaw fit to the temperature profiles of four representative 
clusters: (1) a long exposure of a cluster that has a central tem- 
perature drop, (2) a long exposure of a cluster without a cen- 
tral temperature drop, (3) a short exposure of a cluster with a 
central temperature drop and (4) a short exposure of a cluster 
without a central temperature drop. If the core radius in the 
powerlaw was consistent with zero and/or the temperature gra- 
dient (F in column 5 of Table[T} was positive, no central region 
was removed when determining T wu . This method allowed us 
to take full advantage of the largest possible signal, while being 
certain that T vu was not biased by cool central gas. Our results 
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Fig. 2. These figures show the broken powerlaw fit to four typ- 
ical temperature profiles representing long and short exposures 
of clusters with a central temperature drop and clusters with- 
out a central temperature drop. The dashed blue lines indicate 
the fit and the core radius. Top left: A2029 - a long exposure 
cluster with a central temperature drop. Top right: A3 158 - a 
long exposure of a cluster without a central temperature drop. 
Bottom left: A358 1 - a short exposure of a cluster with a central 
temperature drop. Bottom right: ZwC11215 - a short exposure 
of a cluster without a central temperature drop. 



suggest that using a fixed fraction of the virial radius is not ef- 
ficient and can be dangerous since the core radius, as defined 
above, does not scale with the cluster virial radius. 

Since we include all data outside of the core radius when 
determining 7\,j r , there are regions where residual background 
becomes important. This was especially true for some of our 
more distant clusters (e.g. RXJ1504, z = 0.2153). The resid- 
ual background results from three factors: (1) since the BSB is 
scaled to remove the particle background by matching the BSB 
9.5 - 12 keV rate to that of the observation, if the scaling fac- 
tor is much different from unity, the cosmic X-ray background 
(CXB) in the BSB will be significantly increased or decreased 
before it is subtracted, (2) likewise, the soft CXB component 
varies over the sky (e.g. lKuntjl200lh so the CXB in the obser- 
vation and BSB will be different and (3) in back-illuminated 
(BI) CCDs, there are sometimes residual soft flares. To com- 
pensate for these effects a second soft thermal component was 
included as well as, in the BI CCDs, an un-folded powerlaw 
(i.e. a powerlaw that is not folded through the instrument re- 
sponse). The temperature of the soft thermal component was 
not allowed to exceed 1 keV with a frozen solar abundance 
and zero redshift. The normalization of the soft component was 
also allowed to be negative in the case of oversubtraction of the 



CXB (e.g. for a large scaling factor). For all CCDs of the same 
type (i.e. front or back illuminated) in a given observation, we 
tied the residual background components with the normaliza- 
tion scaled by area. We emphasize that for most clusters this 
step was taken as a precaution since the cluster emission dom- 
inated over the residual CXB. Only in the most distant clusters 
of our sample, where the outer regions are beyond ~0.5 r va , is 
this correction essential. As an example the best-fit temperature 
for RXJ1504 without this correction is 7.1 - 7.9 keV versus 8.4 
- 10.9 keV with the residual BSB correction. On the other hand 
for A0085 (z = 0.055) no significant difference in temperature 
is found with or without the correction. 

Once we measured r v j r , we used it to determine a charac- 
teristic cluster radius and mass. I n the case of radius we used 
the scaling relation determined by Evrard et al. ( 19961) : 



h-\ MpcJ 



= Ni 



10keV 



(3) 



where j is the average overdensity, within rj, above the critical 
density (p c = 3//q/87tG) desired and Nj is the normalization. 
We use normalizations A^oo = 2.00 and A^so = 2.75, and de- 
fine r v ; r = r^Q. We note that these values may be overestimates 
of r=;()() and r^o since th ey are based on early simulations (e.g. 



Finoguenov et al 1l200ll) . For the purposes of core studies, how- 
ever, this approximation should provide an adequate approxi- 
mation of a characteristic cluster radius. 

To estimate a characteristic clust e r mas s (M500), we used 
the formulation of iFinoguenov et al] (1200 II) derived from the 
ROSAT and ASCA observations of the HIFLUGCS clusters, 



10 13 a- 1 m ; 



= (2.5 ± 0.2) 



kT vir 
1 keV 



1 .676±0.054 



(4) 



2.4. Central Temperature Drops 

To quantify the central temperature drop we calculated two 
quantities shown in columns 4 and 5 of Table [2] One is simply 
the temperature of the central annulus (To) divided by the virial 
temperature (7\,; r ). The second is the slope of the powerlaw fit 
to the central temperature profile. For clusters with a measured 
central temperature drop (Sect. I2.3l l. we fit the slope out to the 
break in the powerlaw. For those clusters with no central tem- 
perature drop, we fit the three innermost annulU The purposes 
of this measurement was (1) to check for a universal central 
temperature profile shape and (2) differentiate between a true 
systematic drop and random temperature fluctuations seen in 
merging clusters (e.g. A0754). 

2.5. Density Profiles 

As with the temperature profile, the surface brightness pro- 
file becomes uncertain at large radii due to uncertainties in the 
residual background. Since we are only interested in the cen- 
tral regions and the surface brightness falls off rapidly at large 



In a few cases there were clusters with central drops that did not 
have three bins out to the core radius. In these cases we also fit the 
powerlaw to the inner three annuli. 
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radii, we focused on fitting the central regions only. We argue 
that any uncertainty in the shape of the profile in the outer re- 
gions has a negligible effect when deprojecting the central re- 
gions. Since one of the parameters we are interested in, cus- 
piness (see Sect. 12.91 and Eqn. [TJ, is defined in terms of the 
derivative of the log of the density profile at 0.04 r5oo,we de- 
cided to use a slightly larger (20%) region for determining 
the density profiles. This way we would have a constraint on 
the derivative of the log of the density profile at 0.04 r$oo- 
Specifically, we extracted a spectrum from the projected cen- 
tral region (0-0.048 rsoo) and fit it with an absorbed thermal 
model (WABS*APEC*EDGE). We used the best-fit parame- 
ters of this model to create a spectrumftlM, which we used in 
turn to create a weighted exposure map 10 ! 

We then created a background subtracted, exposure cor- 
rected imag e, in the 0.5 - 7.0 keV range, similar to the method 
described in Hudson et al. I (|2006l). The difference in our method 
here and that described in lHudson et alj ( 120061) is that: (1) we 
only used a single weighted exposure map instead of many 
monochromatic exposure maps for different energy bands, 
(2) we created an error image (cn mg ) and (3) we created a back- 
ground subtracted image with no exposure correction. The er- 
ror image was created assuming the errors in the observation, 
background and readout artifact or out-of-time events (OOTs) 
were Poisson ( a/AO and could be added in quadrature. 

We used the background subtracted image with no exposure 
correction to determine annuli with at least 500 source countj^l 
per annulus and extracted a surface brightness profile and errors 
from the background subtracted exposure corrected image and 
the error image respectively. The surface brightness profile was 
then fit to a single f3 model: 

-3,8+1/2 

2 = S 



or a double (3 model: 



2 = So 



X C , 



-301 + 1/2 



+ En 



1 + 1 — 

Xc, 



-3/32 + 1/2 



(5) 



(6) 



( Cavaliere & Fusco-Femiano 19761) . where x c is the core ra- 
dius, So (So = Soi + S02 for a double f3 model) is the central 
surface brightness and f3 is constrained to be less than 2. 

From (3 and x c for a fit to a /? model (J3\, fa, x Cl and x Cl for 
a fit to a double f3 model), we extracted the shape of the density 
profile. 



1 + 



and 





i 








-3ft \ 


n = 


"01 


HI 


+ "02 



















(7) 



(8) 



for a single and double (3 model respectively, where «o is the 
central electron density, r c is the physical core radius associ- 
ated with x c for our cosmology and «o = ■Jr^ 1 + nf )2 - ^ n or( ^ sr 
to calculate the normalization («o) of the model, we used the 
best-fit spectral parameters from the 0-0.048 r^oo regions, as 
described below. 

For a single {3 model the central electron density, «o is: 



IO'Mtt D a D l £ N 
EI 



(9) 



where N is the normalization of the APEC model, D A is the 
angular diameter distance, Dl is the luminosity distance, f is 
the ratio of electrons to protons (~ 1.2) and EI is the emission 
integral divided by the central density: 



Inserting Eqn.|7]into Eqn.fTOl we get: 



EI = 2tt 



J-oo Jo 



1 + 



+ p 



2\-3s 



dxdl 



(10) 



(11) 



for the 0-0. 48 r 5 q o region fi tted with the APEC model 



dMewe et all 1 19851 1 19861: ISmith & Brickhousd l200d: 



Smith et alj|200lh and R = 0.048 



f"500- 



no = 



For a double [3 model the central electron density, «o is: 
" 10 14 4tt (S 12 LI 2 + LI0 D A DUN' 1 



S 12 LI 2 Eli + LIi EI 2 



(12) 



with the same definition of variables as Eqn. [9] Additionally, 
S12 is the ratio of the central surface brightness of model- 1 
to model-2, EI, is the emission integral for model-/ as defined 
by Eqn. QT| and LI,- is the line emission measure for model-;, 
defined as, 



v-2 



(13) 



LI, = f (1 + 3- 1 

J-tM \ X c 

See AppendixlAlfor details on this calculation. 
2.6. Cooling Times 

Once we created the density profiles, we used them, together 
with the cooling function of the best-fit temperature, to estimate 
the average cooling time of the gas using: 



3 (n e + rii)kT 

^cool — T " 



2 « e « H A(r,z)' 



(14) 



9 See http : / /cxc . harvard . edu/ciao/ahelp/mkinstmap . html 
for the definition of a spectrumfile. 

10 See http : // space . mit . edu/ CXC/ docs/expmap_intro . ps . gz 
for more details on this procedure. 

11 In a few cases, if the /? model fit was poor, this was reduced to 
>100 counts. 



where n Y and « e are the number density of ions and electrons 
respectively. The cooling function, A(T,Z), was calculated by 
spline interpolation on a table of values for the APEC model, 
for an optically thin thermal plasma, kindly provided by R. K. 
Smith. We define the central cooling time (CCT) as: 



fcool(O) - -L, , , 

2 nyMT 48 ) 



(15) 
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where « e o and n io are the central electron and ion densi- 
ties, respectively, and is the average temperature of the 
- 0.048 r 50 o region. 

We note here that in order to remove the bias introduced by 
different physical resolutions due to different cluster distances, 
we took any parameter calculated at r — 0, (e.g. no and CCT) 
to be the value at ro = 0.004 r^oa- The exception to this is To, 
which is defined as the average temperature of the central an- 
nulus, which in all cases had a radius < 0.004 rsoo- We define 
r ctr to be the central temperature of the cluster, either T4& or 
To depending on the parameter considered. In general us To 

«r () /r 48 > = 0.83). 
2.7. Entropy 

As is typically done in X-ray studies of galaxy clusters, we 
define entropy as K = kTn~ 2 ^ and central entropy as Kq = 
kTorC 2 ^, where no is the central electron density and To is 
the temperature of the innermost annulus from the tempera- 
ture profile. Additionally we define a purposely biased entropy 
(^bias) that takes advantage of the inherent difference in the 
density distribution of CC and NCC clusters to further sep- 
arate them. Specifically, taking into account the fact that CC 
clusters should have higher density cores, they will, in general, 
have a smaller central annulus (in the temperature profile, see 
Sect. 12.21 ) than NCC clusters (when fitting spectra). That is, the 
10 000 (or 5000) source count threshold for the central annulus 
will be reached at a smaller radius for CC clusters than NCC 
clusters, as the density is higher and typically the temperature is 
lower in this region for CC clusters than for NCC clusters. We 
then define the density as the average density for the innermost 
annulus, 



"annulus 



3 x 10 14 D A D L C,N 



(16) 



where the variables are defined as before and r annu i us is the 
physical radius of the central annulus. Since the central annu- 
lus for clusters with a central peak in density will be smaller, 
its value of ^annulus will be closer to «o than for clusters with- 
out a centrally peaked density. Unfortunately, since we define 
'"annulus by a constant number of source counts there are addi- 
tional biases due to length of the observation and redshift of 
the cluster. We argue, however that central density is a stronger 
effect. As an example we give the comparison of A2204 and 
A3 158. A2204 has a redshift of z = 0.1522 with a good obser- 
vation time of ~ 1 8 . 6 ks . A3 1 5 8 - a cluster without a bright core, 
has a redshift of z = 0.06 and good observing time of ~55.7 
ks. The innermost annulus[3 of A2204 has a radius of 7'.'0 = 
18.3 kpc, smaller than the innermost annulus of A3158 
which has a radius of 42'.'9 = 48.8 h~j kpc. Even in the ex- 
treme case of A3667 (z = 0.0566) that has a good observing 
time (~485.3 ks) more than 20 times longer than A2204, the 
innermost annulus has a radius (r annu i us = 15'.' 9 = 17.0 hz\ kpc) 
comparable to A2204. Based on this intentional bias, we define 



12 For comparison purposes this is the radius for 10000 counts. 
However when making the temperature profile for A2204 we used 
5000 count annuli. 



•^bias = kTon'., where as with Ko, To is the temperature of 
the central annulus from the temperature profile. 

2.8. Mass Deposition Rates 

We define three mass deposition rates: the classically deter- 
mined mass deposition rate (M c i ass i ca i), the spectrally deter- 
mined mass deposition rate (M spec ) and the modified spectrally 
determined mass deposition rate (M spec 2). M c i ass ; ca i is calculated 
from the gas density and tem perature assuming no en ergy in- 
put. Using the formulation of iFabian & Nulsenl I 1 19791) . within 
radius r, 



M c i assica i(< r) =* 



M gas (r) 



fcoolW - fcool(O) 



(17) 



Usually M c i ass j ca i is calculated at a certain radius r coo i, where we 

— 1/2 

have chosen r coo i as the radius at which t coo \ = 7.7 «_ ' Gyr. 
We chose 7.7 Gyr, since it corresponds to the look back time of 
z — 1, a representative period over which a cluster has time to 
relax and form a cooling flow ( IRaffertv et al. I l2006l) . However, 
for uniformity in comparison with our spectral mass deposi- 
tion rates, we took r to be the smaller of 0.048 rsoo and r coo i. 
Furthermore, to check the effect of this cut, we extrapolated 
our density profiles to calculate the classical mass deposition 
rates out to the cooling radius, assuming the average gas char- 
acteristics remain the same and the slope of the profile remains 
at our fit value. The mean cooling radius was ~0.07 rsoo, and 
the mass deposition rates at the cooling radius were indistin- 
guishable, within the errors in Mciassicai, from those calculated 
at 0.048 r 500 . 

M spec is a direct measurement of the amount of 
gas that i s cooling by fittin g the expected line emis- 



Peterson et al 



sion (e.g 

In order to determine M, 
region to an absorbed 
ing flow model 



2003J) from the multiphase gas. 
we fit the 0-0.048 r 500 
model with a cool- 



spec ' 

thermal 



1988) 



dMushotzkv & Szvmkowiak 
(WABS*[APEC+MKCFLOW]*EDGE). The higher tem- 
perature component of the cooling flow model was tied to 
the thermal model and the lower temperature component was 
frozen at 0.08 keV. Given the moderate spectral resolution of 
the Chandra ACIS, we checked the reliability of the measured 
values. Fig. [3] shows the comparison between the values 
obtained with the XMM RGS and the Chandra ACIS for the 
nine cluste rs in common with ou r sample and the 14 clusters 
reported bv lPeterson et alj ( 12003b . Fig.[3]shows that the values 
obtained with the Chandra ACIS are, in all b ut one casj^l 



consis tent with the upper limits obtained by iPeterson et al 



(2003). Based on this we argue that the values obtained with 
the Chandra ACIS are reliable to well within an order of 
magnitude for CC clusters. As we discuss later, this value may 
not be accurate for NCC clusters. 

M spec 2 is determined similarly to M spec , except that the 
lower temperature component is left as a free parameter. M spec 2 
represents the spectral fit to a cooling flow model that is 



13 In the case of A2052, the Chandra ACIS value is -30% higher 
(although only a difference of 3 M Q yr ) than the upper limit obtained 
with the XMM RGS. 
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A262 2A0335 A496 A780 A1795 A2052 MKW3S S1101 A4059 

Cluster Name 

Fig. 3. XMM RGS - Chandra ACIS spectroscopi c cooling rates 
Here we plot the upper limits measured by iPeterson et al 
d2003l) with the XMM RGS (black), compared to our Chandra 
ACIS measured values (thick red) for the nine clusters in com- 
mon. For eight of the nine clusters, the Chandra ACIS gives 
results consistent with the XMM RGS. For the one cluster 
(A2052) for which they are not consistent, the value mea- 
sured with the Chandra ACIS is ~3Q% higher (3 M Q yr _1 ) than 
the XMM RGS upper limit. This comparison corroborates the 
Chandra ACIS spectrally determined mass deposition rates. 



ity (and trimodality in some c ases) using the Kaye 's Mixture 
Model (KMM) algorithm (e.g. lAshman et al.lll994l) . For each 
parameter we used a single covariance for both (all three) sub- 
groups when determining the significance of the rejection of 
the single Gaussian hypothesis. We chose to do this because 
analytic errors are only statistically meaningful when the same 
covariance is used for each subgroup. However, since there is 
no reason to believe the bimodality in any parameter is sym- 
metric, we considered independent covariances for the sub- 
groups when determining the subgroup assignments. Moreover 
the analytic errors when using different covariances for each 
subg roup can still give a good guideline on fit improvement 
(e.g. lAshman et al.lll994l) . even if their true significance is un- 
known. 

Fig. |4] shows histograms of the 16 parameters with 
Gaussians (created from the KMM algorithm results) overplot- 
ted. For each parameter the CC and NCC subgroups were in- 
dependently determined. We constructed Gaussians using the 
means and covariances returned by the KMM algorithm and 
calculated the normalization of the Gaussians so that the in- 
tegral of the Gaussian was equal to the area of the bins (sum 
of the number of clusters in each bin times the width of the 
bin) within its relevant subgroup. Bins of clusters from the CC 
subgroup are colored blue and bins of clusters from the NCC 
subgroup are colored red. In the cases with a third subgroup 
between CC and NCC clusters, the bins are colored black. 



stopped at a lower-limit temperature 7i ow . The interpretation 
of this model is discussed in more detail in Sect. 14.11 



2.9. Cuspiness 

Vikhlinin et al. (2007) suggested cuspiness as a proxy for iden- 
tifying cool-core clusters at large redshift. Eqn.Q]gives the def- 
inition of cuspiness. For our low redshift sample, we can test 
the correlation of cuspiness with the parameters used to define 
a cool core. Since our density profile is based on a model, we 
define it in terms of the model parameters, with r = 0.04 r^QQ. 
For a single ft model 

3/3r 2 



(r 2 + r 2 ) ' 
For a double (3 model 



a = 3r z 



£ 12 LI 2 £i + Lhb 2 



(19) 



where S12, LI, are defined as before and b; = |l + ^ j 



and V. = 1 + 



3. Results 

3.1. Bimodality and Histograms 

In order to determine the best parameter to separate CC clus- 
ters from NCC clusters, we tested 16 parameters for bimodal- 



3.1 .1 . Central Surface Brightness (A): S 

The histogram for So appears to have two peaks on either side 



of ~10 photons cm s arcsec (0.5 - 7.0 keV). They are, 
however, not well separated and the KMM-algorithm does not 
reject a single Gaussian hypothesis for this data set (<85% con- 
fidence). If divided into two Gaussians, the KMM algorithm 
finds a partition at ~0.80 x 10~ 6 photons ctrT 2 s _1 arcsec -2 , 
which divides that sample into 37 CC clusters and 27 NCC 
clusters. The centers of the subgroups are 3.88 x 10~ 6 photons 
cm -2 s _1 arcsec -2 and 0.32 x 10~ 6 photons cirT 2 s _1 arcsec -2 
for the CC and NCC clusters respectively. 



(18) 3.1 .2. p Model Core Radius (B): r c (% r 50u ) 



In the case of a double ft model the inner core radius was used. 
The histogram of r c as a percentage of rsoo shows little evidence 
of bimodality, with four distinct peaks, however the KMM al- 
gorithm rejects the single Gaussian hypothesis at ~99% con- 
fidence. Using the KMM algorithm we find a split at ~1.3 x 
10~ 2 rsoo, which divides the sample into 30 CC clusters and 34 
NCC clusters. The centers of the subgroups are 0.35 x 10~ 2 r^oo 
and 4.31 x 10~ 2 rsoo for the CC and NCC clusters respectively. 



3.1.3. Central Electron Density (C): n 
The histogram for «o appears to have a single peak at ~ 2.5 



10" 
10" 



,1/- 



71 

l,U2 



cm 



-3 



with a lower, wider distribution below ~1 
ctrT 3 . The KMM algorithm rejects the single 
Gaussian hypothesis at >99% confidence. We separate the two 
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4.5 x 10 2 /j'( 2 cm 3 which divides the sam- 



distributions at 

pie into 37 CC clusters and 27 NCC clusters. The centers of 
the subgroups are 4.02 x 10~ 2 h\ 12 cm" 3 and 0.46 x 10~ 2 h\j 2 
cirT 3 for the CC and NCC clusters respectively. 



3.1 .4. Central Entropy (D): K 

The histogram for Kq looks like it may have a tri-modal distri- 
bution with peaks at ~ 10 h 7 \^ keV cm 2 , ~ 40 h~l n keV cm 

1/3 

and~ 250/i7," J keV cm-. The KMM algorithm rejects the sin- 
gle Gaussian hypothesis at more tha n 99% confidence an d the 
likelihood ratio test statistic (LRTS) (I Ashman et all 1994 sug- 
gests Kq is one of the most significant bimodal distributions. 
The KMM algorithm finds a division between the two sub- 

— 1/3 9 

groups at ~25 /z ?1 keV cm", which divides the sample into 27 
C C and 37 NCC clusters. The divi sion in Kq was first reported 
bv lHudson & Reiprichl f2007l) and iReiprich & Hudson! d2007h 
and fits well with the results o flVoit et alJd2008l) . Adding a third 
subgroup, improves the LRTS further and divides the distribu- 

— 1/3 

tion at ~22 and ~150 /? 71 keV cirr, which gives 24 strong 
CC (SCC) clusters, 22 weak CC or transition (WCC) clusters 

— 1/3 

and 18 NCC clusters. The centers of the subgroups are 1 1 h 7l 1 
keV cm 2 , 59 h~l' 3 keV cm 2 and 257 h~l' 3 keV cm 2 for the 
SCC, WCC and NCC clusters respectively. 



3.1.5. Biased Central Entropy (E): ^bias 

Naively one would expect n an nuius < "o (see Eqn.[T6]for the def- 
inition of n a nnuiusX however since n arm uius is not deprojected and 
«o is deprojected, it is possible that n annu i us will be larger than 
no, implying an appreciable amount of emission along the line 
of sight. This emission increases the apparent central density 
if it is not subtracted before determining the central density. In 
general, however, n an nuius will be smaller than no, especially for 
NCC clusters which require a large region to obtain our count 
criterion. The effect of this intentional bias is that NCC clusters 
end up with a large value of Kbias (as can be seen, the maxi- 
mum value of Kbias is larger than the maximum value of Ko). 
On the other hand CC clusters will have a value of K^ias sim- 
ilar to Kq, separating the distribution. The net effect of the bias 
seems to be to shift the 15 transition clusters in the Kq distribu- 
tion toward the NCC peak. The KMM algorithm rejects the sin- 
gle Gaussian hypothesis at >99% confidence. Surprisingly the 
LRTS suggests the bimodal distribution for Kbias is less sig- 
nificant than for Kq. The two CC/NCC subgroups are divided 
at 40 h~l /3 keV cm 2 , which gives 29 CC clusters and 35 NCC 
clusters. The centers of the subgroups are ll/z~ 1/3 keV cm and 
131 /i~ 1/3 keV cm for the CC and NCC clusters respectively. 



3.1.6. Cooling Radius (% r vir ) (F) 



We made two assumptions when calculating the cooling radius: 
(1) fcooi(rcooi) = 7.7 h~l' 2 Gyr (see Sect. 123 and (2) the gas 
at the cooling radius has the density extrapolated from the 0- 
0.048 fsoo density profile. Therefore, any cluster with a CCT 

—1/2 

longer than 7.7 h 7l Gyr, has a cooling radius of zero. Beyond 
that, there is a broad distribution above about 0.02 r^o. The 



KMM algorithm rejects the single Gaussian hypothesis at 98% 
confidence, however the algorithm does not converge if dif- 
ferent covariances are used for each subgroup (possibly be- 
cause of the large peak at zero). Using a common covariance 
for the two subgroups, the KMM algorithm partitions the two 
subgroups at 0.043 7-500 which divides the sample into 37 CC 
clusters and 27 NCC clusters. The centers of the subgroups are 
0.081 rsoo and 0.012rsoo for the CC and NCC clusters respec- 

1/2 

tively. We note that since a cluster with f coo i(0) > 7.7 » 71 will, 
by definition, have cooling radius of zero, it makes a poor pa- 
rameter for comparison studies. 



3.1.7. Scaled Spectral Mass Deposition Rate (G): 



M spec /M 500 



The histogram for M S pecM^500 shows strong evidence of a bi- 
modal distribution. There are 13 clusters with a best-fit value of 
M spec = and there is a second peak at ~10~ 14 hj\ yr _1 . Since 
the KMM algorithm was fitted to the log value of the data, 
the clusters with M spec = 0, were assigned to log[M spec /(M5oo 
10~ 14 h-ji yr~')] = -5. The KMM algorithm rejects the single 
Gaussian hypothesis at >99% confidence. The LRTS confirms 
that the bimodality is one of the most significant of the 16 
parameters. The KMM algorithm identifies 23 clusters with 
the first (NCC) subgroup. However, due to the large variance 
in the first subgroup, the two clusters with the largest value 
of M spec /M5oo are assigned to the first subgroup (see Fig. |4j- 
G). Since this is physically unreasonable, we take the partition 
value to be ~ 0.23 x 10~ 14 hj\ yr _1 , which divides the sample 
into 43 CC clusters and 21 NCC clusters. The centers of the 
subgroups are 1.3 x 10~ 14 h lx yr _1 and 0.045 x 10~ 14 h lx yr _1 
for the CC and NCC clusters respectively. As with cooling ra- 
dius, the existence of zero values makes M S pec/^500 problem- 
atic as a parameter for comparison studies. 



3.1.8. Scaled Classical Mass Deposition Rate (H): 

^classical/^500 

As with cooling radius, the 18 clusters with CCT > 

— 1/2 

7.7 h 7l ' Gyr, have M c i ass ; ca i = 0. Similar to M spec /M5oo, those 
clusters were assigned log[M c i ass j ca i/(M5oolO~ 14 /271 yr -1 )] = - 
5 when applying the KMM algorithm to the log of the data. 
Besides the peak at zero, there seems to be a second Gaussian 
distribution at -20 x 10~ 14 h lx yr -1 . The KMM algorithm re- 
jects the single Gaussian hypothesis at more than 99% confi- 
dence. The LRTS identifies M c i ass ic a i/A^500 as having the most 
significant bimodality. The partition between the two sub- 
groups is ~0.5 x 10~ 14 hj\ yr~', which divides the sample into 
43 CC clusters and 21 NCC clusters. The centers of the sub- 
groups are 20.5 x 10~ 14 h lx yr 1 and 0.0106 x 10~ 14 h lx yr 1 
for the CC and NCC clusters respectively. We note, as with 

— 1/2 

cooling radius, that since a cluster with f C ooi(0) > 7.7 ' will, 



by definition, have M c i ass ; ca i 
comparison studies. 



0, it makes a poor parameter for 
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3.1 .9. Cuspiness (I): a 

The histogram for cuspiness (defined in Eqn.[T} appears to have 
a Gaussian distribution around a mean at ~1.2. There seems to 
also be a broad distribution at ~QA and a sharp peak at ~0.9. 
The single Gaussian hypothesis is rejected by the KMM al- 
gorithm at ~99% confidence. The LRTS suggests that the bi- 
modal distribution is not as significant as for other parame- 
ters (e.g. M c i ass i ca i/M5oo). Using a trimodal distribution is a sig- 
nificant improvement, but produces an unphysical distribution. 
Basically four of the CC clusters centered around ~1.2 are sep- 
arated out of the CC subgroup (i.e. there are CC clusters on 
both sides of the four clusters, as opposed to them being be- 
tween SCC and NCC clusters). For the bimodal distribution, 
the partition of the subgroups is ~0.7 which divides the sam- 
ple into 35 CC clusters and 29 NCC clusters. The centers of 
the subgrou ps are 1.07 and 0-412 f or the CC and NCC clusters 
respectively. I Vikhlinin et al.l (120071) suggested a break of a=0.5 
for CC versus NCC clusters and a >0.7 for SCC clusters. Using 
these cuts there are 45 CC clusters (10 of which are weak CC 
clusters) and 19 NCC clusters. 

We note h ere that for many of th e clusters that we have in 
common with I Vikhlinin et alJ (120071) . we gene rally find larger 



are 0.72 x 10 30 h l/2 
ergs s 



value s of a than they do (based on Fig. 2 in IVikhlinin et al 
20071) . There are several possible explanations for this dis- 
crepancy. The most obvious is the different way in which we 
calculat ed ryio- We c alcula ted rsoo from 7\,j, using the for 



mula of Evrard et al . (1996). which, as we noted earli e r, may 



overestimate r^. On the other hand IVikhlinin et a l. (2007) 



estimate rsoo from their mas s model. If our value of r5oo is 
larger than the value used by IVikhlinin et al.l d2007l) . it makes 
sense that our values of a will be larger, si nce the profile 



should steepen around 0.04 r%t) (see Fig. 1 in I Vikhlinin et al 



(2007)). The second major difference is that our values of a 
are derived from surf ace brightness profile models, whereas 
Vikhlinin et al.l (12007b use direct deprojection and derive a 



from density profile models. Finally there are minor points that 
can contribute to the discrepancies, such as the center used to 
create the profile (we both use the X-ray peak, but for some 
NCC clusters this is not well-defined), the energy band for the 
surface brightness profile and the techniques used to create the 
profiles. In the end, we argue that these differences lead to an 
intrinsic scatter in the values of a obtained which are dependent 
on the method used to determine it. 

3.1.10. Scaled Core Luminosity (J): L x l\M m kT vil ] 

We define a scaled version of central luminosity L x - Scaled 
L x is taken from the projected spectral fit to the - 0.048r 5 oo 
region. This luminosity is then scaled by kT^ u and the gas 
mass of the region, calculated from the gas density profile 
(Lx/[A^gas kT v i T ]). The histogram of the scaled Lx, appears to 
be a single distribution with a peak at ~ 0.6 x 10 30 /z'( 2 ergs 
s _1 keV~' Mq 1 . In fact, the KMM algorithm does not reject 
the single Gaussian hypothesis (<40% confidence). If split into 
a bimodal distribution the subgroups are partitioned at 0.4 x 
10 30 h 1 ^ 2 ergs s _1 keV~' M^ 1 which divides the sample into 49 
CC clusters and 15 NCC clusters. The centers of the subgroups 



71 ergs s 





- 1 keV _1 M ' and 0.33 x 10 



30 jjl/2 



keV 1 M^ 1 for the CC and NCC clusters respectively. 



3.1 .1 1 . Central Temperature Drop (K): T /T vil 

The central temperature d rop {TqIT yu ) is often used as an iden- 
tifier for CC clusters (e.g. I Sanderson et aL 2006a; Burns et al 



I2008I) . In the histogram there seems to be a large peak centered 
at ~1 and a broad range of temperature drops that peaks around 
~0.4. The KMM algorithm rejects the single Gaussian hypoth- 
esis at >99% confidence and the LRTS suggests that T^/T^ has 
one of the most significant bimodal distributions. The partition 
between CC and NCC clusters is ~0.7, which divides the sam- 
ple into 24 CC clusters and 40 NCC clusters. The centers of the 
subgroups are 0.44 and 0.97 for the CC and NCC clusters re 



spectiv ely. For their sample of simulated clusters, iBurns et al 



d2008l) made a cut at To/T vl[ = 0.8, which also seems to be a 
plausible cut for our sample. 

3.1 .12. Slope of the Inner Temperature Profile (L) 

The histogram for the slopes of the inner temperature profiles 
does not seem to show any bimodality, suggesting that there is 
no universal central temperature profile. The fact that Tq/T^- 
shows some evidence of bimodality hints that the lowest tem- 
perature of the gas in the center of the clusters never cools much 
below ~0.4 r v ; r . However since the slopes show no evidence of 
bimodality, the gradient of this temperature drop from to 
0.4 r v i, is not universal. This would imply that the size of the 
region where the temperature drops below T vw is not universal 
(see Sect. I4.2I for a more detailed discussion). Usually, NCC 
clusters are thought to have no temperature drop, but the dis- 
tribution of log(ro/r v j r ) is continuous from to — 0.2, mak- 
ing it difficult to define a clean break. We also do not find any 
pe ak or distribution around the universal value of -0.4 found 
by ISanderson et al.l d2006al) . The KMM algorithm rejects the 
single Gaussian hypothesis at >99% confidence, but the LRTS 
suggests the bimodality is not as significant as in other param- 
eters (e.g. TqIT vu ). The CC subgroup has a very large variance, 
suggesting a broad distribution in the inner r -profile slopes for 
CC clusters. In fact because of the large variance in the CC sub- 
group, the KMM algorithm identifies the three clusters with the 
steepest centrally increasing temperature profiles as belonging 
to the CC subgroup (see Fig.|U-L). Since this is unphysical, we 
take the partition value to be -0.23, which divides the sample 
into 15 CC clusters and 49 NCC clusters. The centers of the 
subgroups are -0.24 and -0.08 for the CC and NCC clusters 
respectively. 

3.1 .13. Ratio of Central Temperatures in the Soft Band 
to the Hard Band (M): [T (0.5 - 2.0 keV)]/[r 
(2.0 - 7.0 keV)] 

We also considered fitting the same region to different energy 
bands in order to identify cool cores. The idea is that if there 
are many temperatures in the central region the soft band will 
be more sensitive to the cool gas, whereas the harder band will 
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be more sensitive to the hotter gas. We took our central annulus 
from the r-profile and fit it in the 0.5 - 2.0 keV band and then in 
the 2.0 - 7.0 keV band and found the ratio. The distribution of 
values does not look bimodal at all and appears to be a Gaussian 
centered on ~0.7. However, the KMM algorithm rejects the 
single Gaussian hypothesis at >99% confidence. On the other 
hand, the LRTS suggests the bimodality is not as significant 
as in other parameters (e.g. Tq/T^). As with Mspec/Msoo and 
the slopes of the central temperature profiles, the large vari- 
ance of one of the subgroups (the NCC subgroup in this case) 
causes some unphysical assignments (see Fig. @]-M). The four 
cases with the smallest ratio of [7b (0.5 - 2.0 keV)]/[r (2.0 - 
7.0 keV)] are assigned to the NCC subgroup. Using only the 
clusters with large values of [T (0.5 - 2.0 keV)]/[r (2.0 - 7.0 
keV)] that the KMM algorithm assigns to the NCC subgroup, 
the CC/NCC clusters are partitioned at ~0.86. This divides the 
sample into 50 CC clusters and 14 NCC clusters. The centers of 
the subgroups are 0.72 and 0.87 for the CC and NCC clusters 
respectively. 



3.1.14. Scaled Gas Mass within 0.048r 500 (N): 

M g JM 500 

The high density in the centers of CC clusters suggests that 
there should be relatively more gas in the center compared with 
NCC clusters. The distribution, however, does not look partic- 
ularly bimodal, suggesting that size of the core does not scale 
with mass. There seems to be a peak at ~0.4 x 10~ 3 , with sev- 
eral smaller peaks on either side. The KMM algorithm does 
not reject the single Gaussian hypothesis (<76% confidence). 
If split into two distributions they are partitioned at ~0.45 h' 3 ^ 2 
x 10~ 3 , which divides the sample into 38 CC clusters and 26 
NCC clusters. The centers of the subgroups are 0.87 /j~ 3 / 2 x 
10 3 and 0.28 h' 3 ' 2 x 10 3 for the CC and NCC clusters re- 
spectively. 



3.1.15. Scaled Modified Spectral Mass Deposition 
Rate (O): M svec2 /M 50 o 

The histogram for M spec 2M^500 is similar to the other M 
histograms although the upper peak is broader and not as 
pronounced. Interestingly, the scale of Mspec2/A^500 is sim- 
ilar to that of Mciassicai/A^soo about an order of magnitude 
above M spe cM^500 (see also Fig. [TUland Fig. [TBI ). As with the 
other M histograms, the clusters with M spec 2 = 0, were as- 
signed to log[M S p ec2 /(M 5 oo/10~ 14 hn\ yr 1 )] = -5 when apply- 
ing the KMM algorithm to the log of data. The KMM algo- 
rithm only marginally rejects the single Gaussian hypothesis 
(~90% confidence). The two subgroups are partitioned at 1.5 
x 10~ 14 h-j\ yr , which divides the sample into 44 CC clusters 
and 20 NCC clusters. The centers of the subgroups are 1 1.7 x 
10~ 14 h 71 vr" 1 and 0.39 x 10~ 14 h lx vr" 1 for the CC and NCC 
clusters respectively. 



3.1.16. Central Cooling Time (P) 

The histogram for CCT looks similar to the histogram of Kq, 
which is not surprising since they are calculated from similar 
quantities: T ct £j and rig. Like Kq, the histogram of CCT has 
two peaks with a smaller peak between the two. The KMM al- 
gorithm rejects the single Gaussian hypothesis at >99% confi- 
dence and the LRTS suggests the bimodality of CCT is the sec- 
ond most significant (after M c iassicai/A^50o) among the 16 param- 
eters. As with Kq, adding a third subgroup increases the LRTS. 
For a bimodal distribution the partition is ~5 Kj\ Gyr, which 
divides the sample into 42 CC clusters and 22 NCC clusters. 
For the tri-modal distribution the CC/NCC partition remains 
at ~5 h~\ 12 Gyr and the SCC/WCC partition is ~1 h~\ 12 Gyr. 

— 1/2 

There are four clusters with CCT between 5 « 71 Gyr and 

— 1/2 

7.7 « 71 Gyr, so that the difference between a partition at 
5 hZ\^ 2 Gyr and 7.7 h^ 2 Gyr is a matter of low number 

— 1/2 

statistics. Since 7.7 h ll Gyr corresponds to a look back time 
for z~l, about the time most clusters would have time to re- 
lax and form a cool core, we decided to take the partition 

— 1/2 — 1 11 

as 7.7 h 7l Gyr rather than 5 « 71 Gyr. Moreover we used 

— 1/2 

7.7 Gyr to determine M c i ass icai and cooling radius. Using 

1/2 

1 and 7.7 « 71 ' Gyr as cuts we divide the sample into 28 
SCC clusters, 18 WCC clusters and 18 NCC clusters .The cen- 

1/2 — 1/2 

ters of the subgroups are 0.45 « 71 Gyr, 1.91 « 71 Gyr and 
1 1 .2 /i 7l l/2 Gyr for the CC, WCC and NCC clusters respectively. 

3.2. The Defining Parameter 

The likelihood ratio test identifies M c iassicai/A^50u as having the 
most significant bimodality, however we did not choose it as 
the best method to separate NCC and CC clusters. There are 
two reasons for rejecting it as the best method. (1) Clusters 

— 1/2 

which have a CCT > 7.7 h 7l Gyr have M c i ass i ca i = 0, making it 
difficult to compare M c iassicai/^500 to other parameters. (2) The 
errors in M c iassicaiM^500 are quite large (see Fig.|6}G). In fact the 
average uncertainty in Mciassicai/Msoo is ~60% versus 15% for 
CCT, the parameter with the next most significant bimodality. 
This uncertainty is not accounted for in either the histogram 
or the KMM algorithm and therefore the significance of the 
bimodality and the cluster assignments are also quite uncertain. 

The parameters with the next highest significance of bi- 
modality are CCT and Kq. As noted earlier, Kq and CCT are 
both calculated from hq and T ctr (in general * Tq) and so 
it is not surprising that they have similar distributions. Fig. [5] 
shows the tight correlation between Kq and CCT. The dashed 
lines show the cut between CC and WCC clusters for each pa- 
rameter and the dot-dashed lines show the cut between WCC 
and NCC clusters. There are four clusters which are classi- 
fied as WCC clusters when using Kq and are classified as 
SCC clusters when using CCT. It is interesting to note that 
these four clusters have the lowest value of Kq for any of 
the WCC clusters (as classified by Kq) and there is also a 

— 1/3 2 

clear break (at ~30 h 7l keV cm ) between these clusters 
and the other WCC clusters (see Fig. [5}. As noted earlier this 



Tot = r 48 (see Sect. 12. 6t for CCT and T (x r 48 ) for K . 
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Fig. 4. Histograms of 16 parameters that may be used to distinguish between CC and NCC clusters. Blue bins represent CC 
clusters, red bins indicate NCC clusters and where they appear, black bins represent transitional or weak CC clusters (see text). 
Row-wise left to right and starting from the top row the histograms are: (A) central surface brightness (So), (B) (3 model core 
radius scaled by 7-500, (C) central electron density (no), (D) central entropy (Kq), (E) biased central entropy (A'bias), (F) cool- 
ing radius, (G) spectral mass deposition rate scaled by M500 C^spec/^50u)> (H) classical mass deposition rate scaled by M500 
C^ciassicai/A^50oX (I) cuspiness (or), (J) bolometric X-ray luminosity of the - 0.048 r^oo region scaled by the gas mass of the 
region and virial temperature (Lx/{M gas kT vlI )), (K) central temperature drop (To/TW), (L) slope of inner temperature profile, (M) 
central temperature measured in the soft band divided by central temperature measured in the hard band, (N) gas mass within 
0.048 rsoo scaled by Msoo,(0) modified spectral mass deposition rate scaled by M500 (M spe c2/^50o) an d (?) central cooling time 
(CCT). For the cuspiness histogra m, the dash-dot and d ashed vertical line indicates the CC/NCC cut and the weak CC, strong 
CC cut, respectively, suggested by IVikhlinin et al.l (120071) . 



break has also been repor ted by Hud son & R eiprichl (120071) . 
Reiprich & Hudson! d2007h and lVoit et al.l d2008l) . Additionally 



one borderline WCC cluster is classified as an NCC cluster 
when using Ko and one borderline NCC cluster is assigned to 
the WCC subpopulation when classified with Kq. Other than 



these six borderline cases, all the clusters are assigned to the 
same subpopluations whether CCT or ATo is used to classify 
them. 

Based on Fig. [5] there is not much difference between sort- 
ing the sample using Kq or CCT. Assuming bremsstrahlung 



Daniel Hudson et al.: What is a Cool-Core Cluster? 13 



Table 4. Summary of the KMM algorithm results for the 16 parameters. 



Parameter 


Break 




^2 


#of CC 


# of NCC 




IT 1 t 

Value 1 


CC 
Centroid" 1 " 


NCC 
Centroid 


clusters" 


clusters 










w) 


(v) 


z-o (iu pnotons cm s arcsec j 


A OA 


*2 QO 
J.OO 


A ^20 
KJ.jZ 


J / 


07 


■ . /« /I A-2\ 

rc/rsoo (1U J 


1 . j 


U. JJ 




1i\ 
JU 




/ 1 A— 2 7,1/2 3\ 

no (10 rt-j cm ) 


1.5 


A AO 

4.02 


A /I £ 

0.46 




O "7 


(«7i /3 keV cm 2 ) 


150 (22) 


11 (59) 


257 


46 (22) 


18 


^bias (h 7l ' keVcm 2 ) 


40 


11 


131 


29 


35 


Cooling Radius (% rsoo) 


4.3 


8.1 


1.2 


37 


27 


M spec /M 5 oo (A71 10" 14 yr" 1 ) 


0.23 


1.3 


0.045 


43 


21 


M c i assica i/M 50 o (« 7 i 10 14 yr ') 


0.5 


20.5 


0.0106 


43 


21 


a 


a -7/3 


1.U7 


A A 1 
U.41 


'j ^ 
ij 




ix/[M gas /tr v i r ] (10 JU n 7 'j ergs keV M ') 


0.4 


0.72 


0.33 


49 


15 


JO/J vir 


n 7 


V).y 1 


n ztzt 




/in 


Slope of Central Temperature Profile 


-0.23 


-0.24 


-0.08 


15 


49 


[T (0.5 -2.0 keV)]/[r (2.0 - 7.0 keV)] 


0.86 


0.72 


0.87 


50 


14 


M gas /M 50 o (n 7 f /2 10- 3 ) 


0.45 


0.87 


0.28 


38 


26 


M S pec2/M 50 o (h 7l I0- U yr 1 ) 


1.5 


11.7 


0.39 


44 


20 


1 /2 

Central Cooling Time (n ?1 ' Gyr) 


7.7(1) 


0.45(1.91) 


11.2 


46(18) 


18 



Columns: (1) the parameter, (2) the value of the parameter that separates CC and NCC clusters, (3) the KMM determined centroid value for the 
CC clusters, (4) the KMM determined centroid value for the NCC clusters, (5) the number of CC clusters determined by the KMM algorithm 
and (6) the number of NCC clusters determined by the KMM algorithm. 
1 The number in parenthesis, if any, is the cut between SCC and WCC (or transition) clusters. 

* The number in parenthesis, if any, is the centroid value for the WCC clusters. In these cases the other number is the centroid value of the 
SCC clusters. 

" This number includes SCC and WCC clusters where appropriate. The number in parenthesis, if any, is the number of WCC (or transition) 
clusters. 

" Using the cuts suggested bv lVikhlinin etai] J2007I) of a = 1 (WCC and SCC) and a = 0.5 (CC and NCC), gives 35 SCC clusters, 10 WCC 
clusters and 19 NCC clusters. 



emission (i.e A(r) oc r 1/2 ) and n, « n e = n, the dependency 
of the cooling time on the temperature and density is given by 
fcooi K T 1 ^ 2 n _1 . Therefore CCT is more dependent on n than 
on T, while Kq is more affected by T. Since the determination 
of T is more resolution dependent (i.e. it requires ~100x more 
photons to make a spectrum than one bin in a surface bright- 
ness profile), we argue that f coo i is a better parameter to use. 
Furthermore it is also a more traditional metric and short cool- 
ing times are the physical basis of the cooling flow problem. 

3.2.1. CC Fractions in Redshift and Temperature 

We further investigated the subpopulatio ns in the CCT 
by applying the Wilco xon rank-sum test ( Wilcoxonl 1945 



Mann & Whitnev 1947) to the redshift and distributions 



of the CC (as well as the SCC and WCC subsamples) and NCC 
clusters. The Wilcoxon rank-sum test is a non-parametric test 
that determines whether the values of a parameter of two sub- 
populations derive from the same parent distribution. For in- 
stance, it can be used to determine if the redshifts of the CC and 
NCC clusters come from the same parent distribution of red- 
shifts. The results are presented in Table [5] We find that within 
1 cr, the CC and NCC redshifts derive from the same red- 
shift population. That is, there is no evolution of the CC/NCC 



fraction within our redshift range. Likewise, when comparing 
the redshift distributions of WCC and NCC clusters and SCC 
and WCC clusters, the redshifts, within 1 cr, derive from the 
same population. However, for SCC and NCC clusters there 
is a marginally significant difference (~1.1 cr) in the redshift 
distributions, with NCC clusters having generally higher red- 
shifts. Naively, this seems to indicate that the fraction of SCC 
to NCC clusters is higher at lower redshift. However, when 
examining the T vu - populations, it is clear that CC (and SCC) 
clusters come from a lower temperature population than NCC 
clusters, with ~2 cr significance (~2.4 cr for SCC clusters). 
There is even a significant difference (~1.6 cr) between T vlI for 
SCC and WCC clusters. Only WCC and NCC clusters seem to 
derive from the same population of T vn -. Since our sample is 
flux-limited, low temperature clusters will only appear at low 
redshifts. Since SCC clusters come from a significantly cooler 
population than NCC clusters, it makes sense that they would 
also have a lower redshift distribution in a flux-limited sam- 
ple. Therefore we conclude that within our redshift range, there 
is no evolution in the CC fraction. We also find low tempera- 
ture clusters/groups are more likely to be CC clust ers than high 



tempe ratur e clusters, consiste nt with the results of iBurns et al 
d2008l) and lChen et al1(,2007l) . 
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Fig. 5. This plot shows central entropy (Kq) versus CCT. Since 
both quantities are derived from T ca (T^ as To) and no, they 
are, not surprisingly, tightly correlated. The dashed lines show 
the division between CC and WCC clusters. The dash-dot lines 
show the division between WCC and NCC clusters. The black 
line shows the best fit for all clusters, the blue line shows the 
best-fit line for the SCC clusters and the red line shows the best 
fit to the combined NCC and WCCs. The data points are color 
coded by virial temperature (see scale at right). 

Table 5. Results to the Wilcoxon rank-sum test applied to red- 
shift and T vir for the CC, SCC, WCC and NCC subpopulations. 



Sub samples 


Parameter 


Significance 


Which subsample 


being 


being 


(<r) 


has a larger avg. 


compared 


used 




param. value? 


(1) 


(2) 


(3) 


(4) 


CC-NCC 


redshift 


0.82 


NCC 


SCC-NCC 


redshift 


1.12 


NCC 


SCC-WCC 


redshift 


0.96 


WCC 


WCC-NCC 


redshift 


0.16 


NCC 


CC-NCC 


T 

1 vir 


1.98 


NCC 


SCC-NCC 




2.41 


NCC 


SCC-WCC 


T 

A vir 


1.62 


WCC 


WCC-NCC 


T 

A vir 


0.81 


NCC 



Strong cool-core clusters, owing to their high central densi- 
ties have enhanced central X-ray-emission. This may result in a 
higher chance of their detection and serve as an explanation for 
observing a higher fraction of SCC clusters in the HIFLUGCS 
and other flux- or luminosity-limited samples. Since we have 
a complete sample, this calculation can be done. We simulated 
samples of clusters which follo w the X-ray temp erature func- 
tion given by dN/dV ~ T 3 2 dlkebe et al.il2002l) , in the tem- 
perature range (0.001-15) keV and redshift range 0.00-0.25. 
From the above it is clear that SCC, WCC and NCC clus- 
ters come from the same parent redshift distribution within 
1 — cr standard deviation. Hence, we assigned to the clusters 
random redshifts conforming to the N cc D } L law, where Di is 
the luminosity distance. We calculated the luminosities using 
three different Lx~T vlI relations as determined for each of the 
three categories, the SCC, WCC, and NCC clusters, individu- 
ally (Mittal R. & Reiprich T., in preparation). 

In order to estimate the effect of imposing a flux-limit to 
a mixed sample of SCCs, WCCs and NCCs on their result- 
ing fractions, we applied the HIFLUGCS flux-limit, fx (0.1 - 
2.4) keV> 2 x 10~ n erg s" 1 cirT 2 , to the simulated sample. 
We tried two different input sets. In the first (simplest) case, 
we assumed the SCCs, WCCs and NCCs to have the same 
Lx-Tvjj slope (3.33). The normalizations were fixed to those 
found from the fits to the data. In the second case, we fixed 
the slopes for SCCs, WCCs and NCCs to the fitted values 
(Mittal R. & Reiprich T., in preparation). We find that in both 
the cases the output fractions are indeed biased. In particular, 
the output SCC fraction is higher than the input value. Thus we 
conclude that in reality the SCCs, WCCs and NCCs may oc- 
cur with similar fractions and due to the increased X-ray lumi- 
nosity in SCCs, their observed fraction is higher in the present 
sample. 

■3.3. Cooling Time Compared to Other Parameters 

In order to study underlying physical relations and correlate 
simple observables with those requiring high quality data, we 
compare 14 parameter^ from Fig. [4] to the CCT. The rela- 
tions are plotted in Fig. [6] The points are color coded by virial 
temperature with the color scale cropped at 10 keV. In or- 



Columns: (1) subsamples being compared, (2) the parameter (redshift 
or 7\,j,) being used for the comparison, (3) the significance that the two 
distributions are not consistent with the null hypothesis (they are from 
the same parent distribution), (4) the subpopulation that comes from 
the larger redshift or higher temperature distribution. 



3.2.2. Bias on CC Fractions due to Flux-Limited 
Nature of the Sample 

Flux-limited samples suffer from the well-known Malmquist 
bias, namely that brighter objects have a higher detection rate 
than fainter objects. In this section, we address how this bias 
may affect the observed fractions of SCC, WCC and NCC clus- 
ters. 



15 We already compared Kq to the CCT in Fig. 13.21 so it is not re- 
peated here. 



Table 6. Results of simulations done to investigate the impact 
of selection effects on the observed fractions of SCCs, WCCs, 
and NCCs. 



Case 


Category 


Input 


Output 






Fractions 


Fractions 


Same slope 


SCC 


0.34 


0.437+0.043 




WCC 


0.36 


0.281 ±0.041 




NCC 


0.30 


0.281 ±0.044 


Different Slopes 


SCC 


0.310 


0.445 ±0.041 




WCC 


0.335 


0.275±0.038 




NCC 


0.355 


0.280±0.041 
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der to quantitatively compare the parameters to CCT, we fit 
the relations to powerlawo Since there were errors on both 
the parameters and CCT, we used the bisecto r linear regres- 
sion routine, BCES dAkritas & Bershadvl[i996l) . to fit the data. 
Although some of the parameters (e.g. Kq) had correlated errors 
with CCT, for simplicity we assumed them to be independent. 
Additionally, since the relation between the CCT and the other 
parameters seemed to behave differently for the SCC clusters, 
we separately fit the SCC clusters and the non-SCC (WCC and 
NCC) clusters. The black line gives the best fit for all clusters, 
the blue line gives the best-fit line for the SCC clusters and the 
red gives best-fit line for the non-SCC clusters. Table [7] gives 
the fit values for all 14 parameters (along with Kq) for all three 
lines. 



3.3.1. So 

It has been suggested that the c entral surface bright ness, So, is 
an indicator of CC strength (e.g. lO'Hara et al hood) . We find a 
surprising large amount of scatter in the plot of So versus CCT, 
especially at short cooling times. Since no is derived from So 
and CCT is derived from no, one would naively expect a tighter 
correlation. Since S oc n 2 any scatter in the no versus CCT re- 
lation will be amplified in the So versus CCT relation. Not sur- 
prisingly, the outliers at the low end are the lower temperature 
objects. Since CCT depends both on no and To, cooler gas will 
have a shorter cooling time at the same density. Since cooler 
clusters are more likely to be CC clusters than hotter clusters 
(see Sect. I3.21 i, it is not surprising that there is more scatter for 
the SCC clusters than for WCC and NCC clusters. When using 
So to identify CC clusters, the scatter will lead to misidentifi- 
cation; however, adding information about T w i r may be used to 
help reduce the scatter. 

If So is used to separate the CC and NCC clusters, the WCC 
clusters are split into CC and NCC clusters, but not clearly on 
the basis of CCT. The WCC clusters with the longest CCTs are 
classified as NCC clusters, but the WCC clusters with interme- 



diate CCTs (~4 h~\ 12 Gyr) and short CCTs (-1-2 Gyr) are 
sometimes categorized as CC clusters and sometimes as NCC 
clusters. Fig.|7]shows how the correlation becomes tighter if the 
central surface brightness is scaled by kT w { r . Cutting T.o/kT Y1I at 
1.5 x 10~ 7 photons cirT 2 s keV~' successfully separates the 
SCC and NCC clusters, with five misclassified WCC clusters, 
all but one (A2634) with relatively long CCTs (>3 h^ /2 Gyr). 



-1/2, 



3.3.2. r c (%r 500 ) 

We checked to see if there was any correlation between CCT 
and r c as a percentage of rsoo- We expected a correlation since 
CC clusters usually require a double /3 model to fit their surface 
brightness profile, whereas NCC clusters do not. Therefore we 
expected rjrsoo to be much smaller for CC clusters than NCC 
clusters. Surprisingly the correlation appears to be tighter than 
for So, although there is still a lot of scatter. The KMM deter- 
mined cut would classify four SCC clusters as well as most of 
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16 Some parameters are more nearly linear in the log of the cooling 
time, in which case the parameter was exponentiated before fitting. 



Fig. 7. The scaled central surface brightness So/fcTVir versus 
CCT. The So versus CCT correlation tightens when the So is 
scaled by kT V n-. 



the WCC clusters as NCC clusters. If we put the cut at 0.03 rsoo 
(just after the third peak in the histogram counting from the 
left), however, we would more successfully separate the CC 
clusters and NCC clusters. On the other hand, this parameter 
does a poor job of separating SCC and WCC clusters. An inter- 
esting note is that the coolest clusters (groups) appear to have 
smaller values of r c jrsoo compared to hotter clusters. 



3.3.3. «o 

There is a tight correlation between «o and CCT. This is not 
surprising since, as noted earlier, no is the strongest factor in 
determining CCT. All the SCC clusters are CC clusters when 
determined with no and, similarly, both methods agree on the 
NCC clusters as well. The WCC clusters are split at CCT < 

— 1/2 

2 « 71 ' Gyr into nine CC and nine being NCC when determined 
by no. It is impossible to separate SCC and WCC clusters using 
«o- This suggests that the difference between SCC and WCC 
clusters with similar values of «o is that SCC clusters have a 
lower central temperature. This may be due to the central tem- 
perature drop seen in SCC clusters or simply that these SCC 
clusters are cooler than their WCC counterparts. Considering 
NCC and WCC clusters with similar values of n , the WCC 
clusters must be cooler than their NCC counterparts. 



3.3.4. i^BIAS 

As with central entropy, K mA s seems to form a tight correlation 
with CCT. The KMM break separates the SCC clusters from 
the other clusters (with one exception). Perhaps the most inter- 
esting aspect of Kbias is that the slope seems to break between 
the CC and NCC clusters. Unlike many of the other parame- 
ters this looks truly like a broken powerlaw, as fits to both the 
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Fig. 6. A comparison of CCT with fourteen of the parameters from Fig. [4] Row-wise left to right and starting from the top 
row the plots are: (A) central surface brightness So, (B) f3 model core radius r c as a percentage of r^oo, (C) central density 
«o, (D) biased central entropy /^bias, (E) cooling radius, (F) spectral mass deposition rate M spe c scaled by M500, (G) classical 
mass deposition rate M c i ass i ca i scaled by M500, (H) cuspiness a, (I) scaled central, bolometric, X-ray luminosity Lx/(M gas kT v { r ), 
(J) central temperature drop To/T vn - (K) slope of the central temperature profile, (L) central soft band determined temperature 
divided by central hard band determined temperature [To (0.5 - 2.0 keV)]/[7o (2.0 - 7.0 keV)], (M) central gas mass M gas scaled 
by M500 and (N) modified spectral mass deposition rate M S pec2 scaled by Msoo- The clusters are color coded by virial temperature 
with the colorscale shown in the bottom center. The dotted lines represent the division between CC clusters and NCC clusters as 
determined by the KMM algorithm for that particular parameter. The dashed lines represent the division between SCC clusters 
and WCC clusters. The solid black line is the best fit to all the data, the blue line is the fit only to the SCC clusters (as determined 
by CCT) and the red line is the fit to the WCC and NCC clusters (as determined by CCT). As noted in the text, the Fornax cluster 
is often a strong outlier. See notes on individual parameters and Sect. lC.lOl for specifics. 



SCC and WCC/NCC clusters have similar values at the tran- 
sition between them. This appears to be a good parameter for 
determining CC and NCC clusters, although ideally the redshift 
and observation length biases would need to be removed. The 
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merging cluster A3266 is the outlier with CCT~7.7 h 7 , and 
K mAS -900 ft" 1/3 keV cm 2 . 



3.3.5. Cooling Radius (% r vir ) 

Since the cooling radius is defined at ? cool (r cool ) = 7.7 h^\' 2 Gyr, 
the NCC clusters have no cooling radius by definition. These 
clusters were excluded when fitting the relation between cool- 
ing radius and CCT. There seems to be a trend of increas- 
ing cooling radius with decreasing cooling time for the WCC 
clusters. The SCC clusters, however, do not seem to have any 
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trend between cooling time and size of the cooling region (as 
a fraction of the virial radius). Therefore we can conclude that 
the density gradients for SCC clusters vary greatly. The large 
outliers are: SI 101 at the high end with a cooling radius of 
> 0.15 r 5 oo and the Fornax cluster at the low end. 



3.3.6. Scaled Spectral Mass Deposition Rate 

This plot explicitly demonstrates the danger of using M spec 
to divide the distribution. The error bars that are not visible 
in the histogram and not considered in the KMM-algorithm, 
show how much overlap there is within errors. There seems 
to be only a very weak trend with CCT, especially if the up- 
per limits are included. The clusters with the highest values of 
M spec /M 5 oo, however, are SCC clusters. As noted earlier and 
discussed in more detail later, several of the NCC clusters have 
non-negligible spectral mass deposition rates. This fact along 
with the spread in values for the SCC clusters makes it very 
difficult to use M spec to identify CC and NCC clusters, at least 
with the spectral resolution of the Chandra ACIS. 



3.3.7. Scaled Classical Mass Deposition Rate 

As with cooling radius, M c i ass icai is defined as zero for the NCC 
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clusters (clusters with CCT > 7.7 h 7l ' Gyr), so these clus- 
ters were omitted when fitting the comparison of M c iassicai/^500 
to CCT. There appears to be a weak trend of increasing 
Moiassicai/Msoo with decreasing CCT, albeit with a great deal 
of scatter. We attribute this scatter to the differences in the 
gradients of the density profiles. That is, clusters with a very 
steep density gradient will have a very short CCT but not a 
large value of ^classical* whereas clusters with relatively dense 
cores but flatter density gradients will have a longer CCT but 
a large M c i ass j ca i. Based on the KMM determined cut of 0.5 x 
10~ 14 h~\ yr~' for M c i assica i/M 5 oo, all the SCC clusters are clas- 
sified as CC clusters and all but three of the WCC clusters are 
classified as CC clusters. Two of these WCC clusters are bor- 
derline cases. 

There are six WCC clusters (IIIZw54, A 1650, A 1651, 
A2142, A2244 and A4038) that have M c i assica i/M 50 o > 15 x 
10~ 14 h^l yr , more consistent with SCC clusters, although 
toward the lower end of the SCC clusters. There are five SCC 
clusters (A0262, MKW4, NGC4636, A3526 and A 1644) that 
have Mclassicai/Msoo < 15 x 10~ 14 hz\ yr , more consistent with 
WCC clusters. There are four (one SCC and three WCC clus- 
ters) outliers that have very low values of M c iassicai/A^500 (i.e. 
< 4 x 10~ 14 h~{ yr 1 ). They are: NGC1399, A3266, A3667 
and A2634. It is interesting to note that all but two of the low 
A^dassicai/^500 SCC clusters are low temperature (IcT < 2.5 
keV) systems, in fact they are four of the nine coolest sys- 
tems. For the ten coolest systems, M c i ass j ca i/Mgoo is less than 
45 xlO 14 yr -1 , which is below the average for the SCC 
clusters. The WCC outliers are all merging systems suggesting 
current cool-core survival, but with some disruption of the cool 
core. 



3.3.8. Cuspiness 

Cuspiness shows a lot of scatter when compared to CCT. It 
seems to do a good job in identifying NCC clusters; only 
three NCC clusters have a > 0.5. All the SCC clusters have 
a > 0.5, but not all of them have a > 0.75 as suggested by 
Vikhlinin etail d2007l) . Using the KMM determined cut of a = 
0.7, several of the SCC and all of the WCC would be classified 



as NC C clusters. Likewise, using the cuts of I Vikhlinin et al 



(120071) . seven of the WCC clusters would be SCC clusters, 
three of the SCC clusters would be WCC clusters and four of 
the WCC clusters would be classified as NCC clusters. Looking 
at the fit to the data, there appears to be a shift in normalization 
and slope for the SCC clusters. That is, cuspiness increases 

— 1/2 

with shorter CCTs. However, below CCT~1 « 71 , the cuspi- 
ness drops suddenly (or at least there is a large range of val- 
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ues). Below a CCT of ~1 « 71 Gyr, the cuspiness increases 
more rapidly as the CCTs get shorter. 



3.3.9. Scaled Central Luminosity 

Scaled central luminosity (Lx/[M gas Wv&]) shows a rather flat 
but clear relation with CCT, with large scatter and larger un- 
certainties. Using the KMM determined cut, however, only one 
SCC cluster is misclassified as an NCC cluster. Some of the 
WCC clusters are classified as CC clusters and some as NCC 
clusters when using the KMM determined cut in scaled central 
luminosity. If the cut is r aised to 0.5 x 10 32 /j'f erg s" 1 keV" 1 
Mq , still only one SCC cluster is misclassified, and the num- 
ber of misclassified NCC clusters is reduced to two. Unlike, 
e.g. yiq, there seems to be no trend with the CCT of the WCC 
and whether it is classified as CC or NCC when using scaled 
central luminosity. As with many parameters the relation seems 
to get steeper for the SCC clusters. 



3.3.10. 7yr vir 

The most interesting feature in this plot is the break in the dis- 
tribution for the SCC clusters. Although there is a lot of scat- 
ter, there is a clear drop in central temperature for SCC clus- 
ters. However, the value of the drop seems to be independent 
of CCT. There are a few NCC clusters with central tempera- 
ture drops, most notably A2256. This may be due to a recent 
merger in which a core has not yet been completely destroyed 
(alt hough a possible dense c ore has been disrupted). 



Sanderson et al . ( 2006a) suggested a similar parameter for 



dividing CC and NCC clusters. Their parameter is roughly the 
inverse of Tq/T^ (7-ratio = 7\,i r /T core ). The major difference 
being that T COK is defined as the temperature within O.lrsoo, in 
general much larger than the region we used to measure Tq. 
They then define CC significance, which is the difference (in 
units of <x) between r -ratio and unity (i.e. r "^"°.~ - )■ They de- 
fine a CC cluster as any cluster with CC significance >3 (i.e. 
a significant drop in temperature). Since they use a bright sub- 
sample of HIFLUGCS, all 20 of their clusters are in our sam- 
ple. We find that all their CC clusters are SCC clusters in our 
sample. Their NCC clusters are either WCC or NCC clusters. 
Likewise, employing their method, using T vu -/To, we find good 
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agreement with their classification. On the other hand we argue 
the values of CC significance are very sensitive to the uncer- 
tainties in r () and r v ; r . That is if r -ratio < 1, then the longer 
the observing time the larger CC significance becomes and 
the more likely the cluster will be classified as a CC cluster. 
However, we concur with their result that in general SCC clus- 
ters have a central temperature drop and that WCC and NCC 
clusters have a shallow or no central temperature drop. 



3.3.1 1 . Slope of the Temperature Profile 

Like To/T v i r , this plot shows a break in the distribution at the 
SCC clusters. It is clearer in this plot than in the TolT vu - plot. 
This is most likely because the slope takes into account a sys- 
tematic decrease in temperature as opposed to random temper- 
ature fluctuations. That is, for an NCC cluster with some cool 
gas that has not been completely disrupted, there will be less 
of a gradient than for a cluster that has a systematic decrease in 
temperature in a dense, well-defined core. For the SCC clusters 
there seem to be many slopes independent of CCT, suggest- 
ing no universal central temperature profile. This along with 
To/T v i r demonstrates that studies that define CC clusters based 
on a temperature drop in the core are defining SCC clusters as 
CC clusters and WCC and NCC clusters as NCC clusters. The 
temperature drop and profiles are discussed in more detail in 
Sect. [£2] 



3.3.12. [T (0.5 - 2.0 keV)]/[r (2.0 - 7.0 keV)] 

[T (0.5 - 2.0 keV)]/[J (2.0 - 7.0 keV)] appears to be a very 
poor parameter for identifying CC clusters. As with M spe c, it is 
possible that multiple temperatures along the line of sight (e.g. 
from a merging system) can produce results that make an NCC 
cluster appear to be a CC cluster. In fact the best fit (albeit a 
very poor fit) shows a positive trend rather than the expected 
negative trend. Additionally, an NCC cluster (A1656) shows 
the smallest fractional value. 



3.3.13. Mgas/Msoo 

M gas /M 5 oo shows a very promising trend for the WCC and NCC 
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clusters; however, at CCT<1 h lx Gyr there appears to be no 
relation. In fact, the best fit to the SCC clusters predicts a de- 
crease in MgaJMsoo with decreasing CCT. The outliers appear 
to be cool clusters and groups, which are less dense in the core 
(short cooling time comes partially from low temperature) and 
possibly have steeper profiles (leading to less gas mass within 
0.048 r5oo than for rich clusters). Fig.[8]shows how dividing by 
kT v i r lessens the scatter for the SCC clusters, explicitly show- 
ing that the large scatter is caused by cool clusters. Of course, 
since the correlation is already tight for NCC and WCC clus- 
ters, dividing by kT vn results in more scatter for those clusters. 




Fig. 8. This figure shows that that M ga JM 5 oo versus CCT cor- 
relation tightens for SCC clusters if scaled by kT v \ r This sug- 
gests that cool clusters cause the scatter in this relationship. 
Unfortunately, but not surprisingly, the tight correlation be- 
tween Mgas/M 5 oo and CCT for WCC and NCC clusters be- 
comes worse. 

3.3.14. Scaled Modified Spectral Mass Deposition 
Rate 

As with MspeJMsm, the large uncertainty in M spec 2/^500 shows 
the limitations of the histogram. However, M spec 2/A^500 does 
show a trend with cooling time, albeit with a large scatter. The 
fact that it shows a trend is interesting since CCT and M spec 2 
are independently determined. See Sect. I4.ll for a further dis- 
cussion on M spec 2. 

3.4. What is the Best Diagnostic of CC Clusters? 

One large problem in studying the evolution of CC cluster frac- 
tion is that resolving the cores of clusters is difficult for distant 
clusters. In sect. I3.3l we compared various parameters to CCT. 
Not surprisingly «o and Ko provide the tightest correlation to 
CCT. In distant clusters no is slightly easier to determine than 
CCT and Ko since it requires no spectral information from the 
central region. On the other hand it is still dependent on one 
being able to resolve the very center of the cluster, so it is not 
the ideal parameter to determine whether a distant cluster is a 
CC or NCC cluster. 

Unfortunately, when compared to no and Ko, most of the 
other parameters show a lot of scatter in their relation with 
CCT. Of those other parameters /Tbias appears to show the 
tightest correlation. Unlike Kq, /^bias is n °t dependent on no, 
but on the average density in a somewhat larger region. It is 
therefore not limited by the need to resolve the central region. 
There are, however, two problems with using /^bias as a proxy: 
(1) it requires a large enough region to fit a spectrum (which 
can be physically quite large for distant clusters) and (2) the 
size of the region from which /^bias is determined differs be- 
tween observations (i.e. it is dependent mainly on observing 
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Table 7. The best-fit parameters for the 15 investigated CCT versus parameter. The parameters were fit to CCT = CCTo 
(parameter) r . 



Parameter 


r 


CCTo 


r 


CCTo 


r 


CCTo 




All 


All 

h-\ 12 Gyr 


sec 


sec 

h'l' 2 Gyr 


WCC+NCC 


WCC+NCC 
h^l 12 Gyr 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


K /(100 h~[ r> keV cm 2 ) 


1.085±0.036 


4.3±0.6 


0.843±0.109 


2.7±0.8 


1.278+0.073 


3.9+1.3 


Eo/(10~ 6 photons cirT 2 s - ' arcsec" 2 ) 


-0.894±0.061 


2.6±0.3 


0.991±0.737 


0.090+!!"? 

-0.063 


-0.830+0.073 


3.3+0.4 


r c /(r 5 oolO" 2 ) 


0.932±0.060 


1.6±0.2 


0.378±0.068 


0.63±0.06 


0.698+0.085 


2.8+0.5 


«o/(10~ 2 hi' 1 cm" 3 ) 


-1.194±0.051 


3.5±0.2 


-0.860±0.151 


1.8+0.4 


-1.016+0.046 


4.1+0.2 


A"ria<;/(100 /?;! /3 keV cm 2 ) 

ll DlAj / ^-i-V/V/ ""7 1 lVW T Will ^/ 


0.975±0.062 


4.7±1.2 


0.501±0.081 


1.5+0.3 


1.314+0.166 


4.4+3.3 


exp[Cooling Radius/(10~ 2 r v j r )]^ 


-0.350±0.025 


11.1+1.1 


-0.584± 1.055 


70.46+" 3 ;? 6 

-70.45 


-0.336+0.026 


12.2+1.0 


(Msnfr/M, m )/(10 -14 /ivi yr -1 ) 

\ L,± cjpec/ J " jUU // V A " ' 1 / 1 J 1 / 


-0.598±0.121 


1.5±0.4 


-0.996±1.313 


0.78+0.33 


-0.404+0.163 


4.0+0.9 


(Mcia S sicai/M 5 no)/(10" 14 fc, 1 yr _1 ) 


-0.816±0.742 


1.1 + " 

-1.0 


0.966±0.675 


0.018+° ° 4 ° 

-0.016 


-0.280+0.163 


4.1+1.3 




-3.726±0.198 


34±6 


-1.610±0.268 


2.8+0.9 


-3.589+0.345 


35+7 


L^/(10 32 ^( 2 ergs s" 1 keV" 1 Mz l ) 


-3.057±0.250 


0.41±0.07 


-1.114+0.167 


0.43+0.04 


-3.795+0.699 


0.24+0.14 


expfro/Tvi,] 1 ' 


4.931±0.365 


0.043±0.013 


-0.102±0.290 


0.76+0.10 


5.047+2.421 


04 i+°- 099 

"•"t-0.037 


exp[Slope of r -profile] f 


6.908± 1.209 


5.2±1.0 


-1.208±0.271 


0.32+0.04 


17.48+16.60 


8.4+4.2 


exp[T (0.5-2.0 keV)/T (2.0-7.0 keV)]+ 


-3.469±3.389 


28^ 


-1.027±0.586 


1.0+0.5 


-3.957+2.316 




M gas /(10- 3 h-] 11 M 500 ) 


-1.997±0.136 


0.58±0.10 


0.894±0.229 


0.49+0.06 


-1.832+0.171 


0.93+0.20 


(M Spec2 /M5 00 )/(10- 14 /t 71 yr 1 ) 


-0.375±0.028 


3.3±0.4 


-0.708±0.248 


4.3+3.5 


-0.347+0.058 


6.3+0.8 


E /(kT10~ 7 photons cm" 2 s" 1 keV" 1 ) 


-0.892±0.039 


5.5±0.4 


-0.630±0.221 


2.5+1.6 


-0.914+0.063 


5.7+0.4 


M gas /(10- 4 fcrM 500 keV-') 


-1.654±0.089 


2.9±0.2 


-0.805±0.138 


1.02+0.17 


-1.458+0.127 


3.6+0.4 



Columns: (1) the parameter name, (2) the fitted slope (T) including all clusters, (3) the fitted normalization (CCTo) including all clusters, 
(4) T for the SCC clusters, (5) CCT for the SCC clusters, (6) T for WCC and NCC clusters and (7) CCT for the WCC and NCC clusters. 
' For this parameter the fitting function was CCT = CCT (exp[parameter]) r 



time and redshift). As a solution to (1) we suggest that the 
normalization of the thermal model is not as difficult to con- 
strain as the temperature or metalicity (not as many counts are 
needed). In extreme cases the temperature and metalicity could 
be frozen at their overall cluster values so that only the cen- 
tral normalization would need to be constrained. As a solution 
to (2) we suggest that the region size be defined such that it 
is only dependent on the density profile and is independent of 
redshift, temperature and observing time. Specifically if two 
clusters have identical density profiles, then /^bias should be 
derived from physical regions of the same size. 

We suggest using a constant value of the projected emission 
measure in a cylinder, 



I(R) = In 



oo JO 



«h rdrdl, 



(20) 



to determine region sizes that depend only on the density pro- 
file. Here we define n e as the electron density and tin as the 
proton density. One can simply use the density profile or model 
for a cluster to directly determine R for a predetermined value 
of /. Alternately, for given spectral parameters and distance, the 
ratio of the count rate to normalization (as defined in XSPEC 
for example) is determined, so that the total observed count, 
C, can be determined and used to define R. That is, the count 
threshold for a region of size R that has a given value of / is: 



cnts UOOksA 



100 cnts cm 5 s 



( 2OO/17/ Mpc') 
D A (l +z) 



(21) 



where f Q b s is length of the observation, k is the ratio of the count 
rate to normalization and F is a constant that depends on the 
desired value of /, 

10- 11 ' 



F = 



16tt 



h n } Mpc 2 cm" 



(22) 



V "71 

A value of / = 8.6 x 10 65 cm -3 , corresponding to F — 
17910, produces a 10000 count region (in the 0.5 - 7.0 keV 
band) for a metalicity of Z = 0.25, a photoelectric absorption 
column density of Nh — 2 x 10 20 cuT 2 and our median virial 
temperature (4.3 keV), redshift (0.047) and observing time time 
(44 ks). See appendixlBlfor more details on these calculations. 

Although this method makes the Kbias region size consis- 
tent among clusters, there is an additional problem that a rea- 
sonable value of / at one redshift, will often be unreasonable 
for clusters at a significantly different redshift. For instance, for 
our proposed value of / — 8.6 x 10 65 cm 4 , which leads to 
a 10 000 count region for our median cluster parameters, would 
produce a 58 count region for a cluster with identical parame- 
ters but at a redshift of z = 0.5. However, if the region size is 
not constant for all studies (say the chosen value of / varies de- 
pending on the study), it makes comparisons between studies 
difficult. 

After Kbias, the next best parameter appears to be the 
scaled central X-ray luminosity (scaled Lx). The region from 
which we extracted scaled Lx is well-defined (0 - 0.048 fsoo), 
but suffers from two other problems. The scaling depends on 
M gas , which in turn depends on «q. In addition, as with ^bias, 
a spectrum is needed in order to determine the luminosity. Of 
course, as with /Tbias, in extreme cases the overall cluster tem- 
perature and metalicity can be used when getting an estimate 
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of the core luminosity. The major advantage in using scaled Lx 
over Krias is that it is better defined for this study. 

For distant clusters in which only a few thousand counts 
are available, it is necessary to use So or a. So has the ad- 
vantage that surface brightness profiles only need ~50 counts 
per bin, as opposed to spectra that require thousands of counts. 
Unfortunately there is the problem of resolution. For very dis- 
tant clusters, it may not be possible to resolve So in the very 
central regions. Also, generally if So is measurable then «o is 
also measurable and is a much better proxy. However, if So is 
the only parameter available, it should be divided by IcT vn be- 
fore being used as a proxy (Fig. |7). Cutting So/£7Vir at ~1.5 
x 10~ 7 counts cirT 2 s keV~' does a good job of separating 
NCC and CC clusters. 

This leaves cuspiness as the only viable candidate for dis- 
tant clusters with few counts. Unfortunately, there is a great 
deal of scatter in the a-CCT plot. Moreover, there seems to be 
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a drop in a at CCT ~ 1 Gyr, so that SCC clusters with 

1 II 

CCT~ 1 h~\"- Gyr have the same value of a as the WCC clus- 
ters. Unlike So there does not seem to be any trend with T wu 
and we have not found any parameter which can be used to 
tighten this correlation. Cutting a at 0.75 ensures for the most 
part that the sample above the cut will be CC cl usters, although 



not ne cessarily SCC clusters as suggested by IVikhlinin et al 



(2007). Cutting a at 0.5 seems to capture most of the CC clus- 
ters above the cut, with a slight contamination from the steepest 
NCC clusters. Finally cutting at a-l ensures that the clusters 
above the cut are SCC clusters, although there are still several 
SCC clusters below the cut. In short, it is impossible to separate 
out a sample, but it is safe to say a cluster with a > 1 is an SCC 
cluster and a cluster with a >0.5 is most likely a CC cluster. 



Recently, [SantoseLaLl 42008) suggested a concentration pa- 
rameter (Csb) as a proxy for determining whether a distant 
cluster is a CC or not (see Eqn.|2]i. Although this proxy is be- 
yond the scope of our analysis (the residual background must 
be carefully handled to find the surface brightness at 400 kpc), 
we used our (3 model fits to estimate Csb- Fig.|9]shows the re- 
sults of the fit, which seems very promising. The one strong 
outlier is NGC1399. As discussed in Appendix IC. 101 the sur- 
face brightness profile of the Fornax cluster shows a flattening 
at ~0.04-0.05 r 500 , probably due to NGC1404. This flattening 
will cause a severe overestimate of S(r < 400 kpc) when ex- 
trapolating the profile out to 400 kpc. To do a true comparison, 
we would need to measure the integrated surface brightness, 
rather than use models extrapolated from ~0.05 r^oo- 



4. Discussion 

4.1. The Cooling Flow Problem 

The cooling flow problem is the discrepancy between the spec- 
trally determined mass deposition rate and the classically de- 
termined mass deposition rate. Fig.[10]clearly demonstrates the 
cooling flow discrepancy. In all but three of the 46 CC clus- 
ters, the spectral mass deposition rate is less than the classical 
mass deposition rate, usually by at least one order of magnitude 
(dash-dot line). 
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Fig. 9. The conce ntration parameter Csb suggested by 
Santos et al.l (12008b versus CCT. These preliminary results 
seem very promising. As discussed in the text, the outlier 
(Fornax cluster) is probably due to a problem in extrapolating 
the surface brightness profile. The horizontal blue dash-dot line 



and bl ack dashed line show the suggested cuts of lSantos et al 



(12008b between SCC and WCC clusters (0.155) and WCC and 
NCC clusters (0.075), respectively. 



In two of these three cases, NGC1399 and A3266, the mass 
deposition rates are consistent within errors. Only in the case 
of A2634 is M spec > M classical . As noted in Sect. 13.3.71 all 
three of these clusters (A2634, NGC1399 and A3266) have 
anomalously low values of M c iassicai- The WCC clusters A2634 
and A3266 appear to be involved in mergers (although as dis- 
cussed later, probably not major mergers), the low values of 
^classical being then consistent with a disrupted core. Also as 
noted in Sect. 13.21 M spec can be over-estimated in merging clus- 
ters. Moreover in the case of A3266, it is on the line between 
NCC and WCC clusters with CCT -7.7 h^ /2 Gyr. As seen 
in Table [2] several of the NCC clusters have nonzero values of 



M spec , so it is not so surprising that A3266 would have M spec ~ 
Mclassical' The case of NGC1399 is not as clear. In many of the 
plots NGC1399 is an outlier. As discussed in Appendix |C. 101 
this may be due to a problem with extrapolating the surface 
brightness profile. Even if this is the case, M c iassicai is small 
(< 1 h~l M yr 1 ) and M c i assica i * M spec (see Appendix OO] 
for details). The anomalously low Mciassicai indicates that the 
density drops off quickly. In fact the cooling radius is only 
0.026 rsoo, which is among the four smallest. The other three 
(in order from smallest to largest) are A3266, A3667, A2634, 
all WCC merging clusters, two of which are also outliers in 
this plot. The SCC cluster with the next smallest cooling ra- 
dius is A0262 that has a cooling radius of almost twice that of 
NGC13 99. Although NGC 1399 shows some evidence of merg- 
ing (e.g. iDrinkwater et alJuOOll) . it appears to be just starting 
to merge so only the outer regions have been affected. It is also 
possible that the galaxy NGC 1404 has disturbed the core as it 
fell in. A more interesting scenario is that since NGC 1399 is 
dynamically young (Y. Schuberth, private communication), the 
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Fig. 10. Spectroscopic cooling rate M spec versus classical cool- 
ing rate M c i ass j ca i for the 46 CC clusters. The dashed-line shows 
parity between the two. The dash-dot line and the dotted line 
show one and two orders of magnitude difference respectively. 
In all but three cases, Mspec ^classical' 



cooling flow has just formed so no energy has been injected yet 
and M spec * M dassica i. 

With the exception of these three discrepancies, in the other 
43 CC clusters, M spec <K Mciassicai- Assuming no energy is input 
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to the CF and the CF formed ~1 .1 h 7l 1 Gyr ago, these two 
quantities should be equal. If the latter assumption is incor- 
rect, we can find the cooling radius such that M c i ass j ca i = M spec . 
We define the cooling time at this radius to be the CC forma- 
tion time (CCFT). We note that the CCFT is an upper limit 
to the actual time since formation because, firstly, as the gas 
flows to the center, its density increases causing it to cool more 
rapidly. Therefore the actual cooling time of the gas is shorter 
than f coo i, as defined in Eqn. [T4j Secondly, the spectroscopic 
mass deposition model, MKCFLOW, used to fit the line emis- 
sion from a multiphase gas often provides a better fit than a 
single-temperature thermal model. In other words, it may find 
a cool gas component even when there is none, yielding M spec 
to be greater than its actual value. Fig [TT] shows the histogram 
of CCFT for the 46 CC clusters. Ten of the clusters have only 
upper limits for M spec , for the others the distributi on peaks at 



around 0.5 Gyr, about five to ten times shorter than I Allen et al 



(2001a) found using ASCA and ROSAT data and assuming in- 
trinsic absorption. This means in order for the observations to 
match, in the absence of heating, the CFs would have had to 
form very recently. 

In order to check the plausibility that the cooling flow 
discrepancy is due to the assumption that cooling flows have 
formed very recently, we devised a rough test to estimate this 
likelihood. This test is based on the hypothesis that low- and 
high-z clusters are both drawn from the same underlying pop- 
ulation of clusters. Hence, the high-z clusters in our sample, 
when evolved uninterrupted to allow their cool cores to grow, 
should have mass deposition rates comparable to the low-z 
clusters in our sample. We show below quantitatively, for the 
first time, that the likelihood of this hypothesis is very small. 




Cooling Core Formation Time 



Fig. 11. This figure shows the histogram of cool core forma- 
tion time (CCFT) for the 46 CC clusters. Ten of the clusters 
only have upper- limits for M spec . For there to be no discrep- 
ancy between M c i ass j ca i and M spec , cooling flows would have to 
have formed very recently. 



There are two major limitations to this very simple test. 
Firstly, our sample is flux-limited and hence due to selection 
effects we are likely to pick clusters of increasing mass with 
increasing redshift. Secondly, even though our sample spans 
a limited redshift range, the cluster population changes with 
redshift. This is because, in a hierarchical Universe, the more 
distant clusters (high-mass) will actually never evolve into the 
mostly low-mass clusters we observe locally. These two rea- 
sons together make a comparison of properties between low- 
z and high-z clusters difficult. As an approximate solution to 
the first problem, we scaled the mass deposition rates by Msorj, 
even though this causes a bias against the high-z clusters. This 
is due to the fact that at a fixed gas density, the cooling time is 
an increasing function of the temperature, CCT oc T 1 ^ 2 oc M 1 ^. 
Thus, if clusters are self-similar, we should expect the cooling 
rate to scale with mass more slowly than oc M. Consequently, 
under the assumption that the cooling flows are recent phenom- 
ena, we should not expect the high mass clusters evolved to the 
lookback time of low-z clusters to show higher Mciassicai/Msoo 
than the lower mass ones. 

We divided our sample into four redshift bins, each with 
16 clusters. Although we have a local sample (median redshift 
for the 16 most distant clusters is z ~ 0.08), most of our clus- 
ters span a lookback time, fiookback ~1 GyiQ > (CCFT). 
We forward evolved each cluster in redshift bins 2, 3 and 4 by 
adding to its CC formation time, CCFT , the difference between 
its lookback time, LBT , and the mean lookback time of clusters 
in the first bin, (LBT) l . Hence, 



CCFT' = CCFT + LBT - {LBT) l 



17 This is the lookback time from z = to z ~ 0.08. 
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Fig. 12. This figure shows the forward-evolved M c i ass i ca i for 
high-z clusters in redshift bins 2, 3 and 4 (red open circles) 
and original M c i ass j ca i for low-z clusters in redshift bin 1 (blue 
open triangles). 

where CCFT' is the forward-evolved CCFT of that cluster. In 
order to calculate the cooling rate in these forward-evolved sys- 
tems, we determined the radius, R', at which the cooling time 
of the gas equals CCFT' and calculated M c i ass j ca i at this ra- 

diuS ' ^classical' Sll0Wn in Kg-Elis Classical SCaled M 500 Ver- 

sus redshift. Since the low-z clusters belonging to redshift bin 
1 (blue open triangles) are not evolved forward in time, these 



clusters, after having been evolved to match the physical state 
of low-z clusters, have systematically higher M c i ass j ca i than the 
low-z clusters. Thus, it is very unlikely that the discrepancy be- 
tween classical and spectral mass deposition rates in cool-core 
galaxy clusters can be explained away by invoking the recent 
cool-core formation hypothesis. 

We note there are 27 clusters with zero M', . ,. Of these, 

classical ' 

18 constitute the entire NCC cluster subsample, 6 are WCC 
clusters and 3 are SCC clusters. That the NCC clusters show 
no mass deposition rate even after evolving them forward is 
not surprising. Out of the 9 WCC+SCC clusters, the five high-z 
clusters are ones with extremely short or zero CCFT and evolv- 
ing them forward still results in CCFT' to be shorter than the 
cooling time at R = 0. 

Another interesting exercise is to check the emission mea- 
sure distribution of th e cool gas. In terms of the power radiated 
per unit temperature, IPeterson et alJ (120031) parametrized this 
as: 



#-!£» + i>(f) 

dT 2 pm p \T ) 



(24) 



For a steady cooling flow at constant pressure, we would have 
6-0, but if cooling is retarded then the amount of gas at lower 
temperatures is reduced (6 > 0). With our limited spectral reso- 
lution we approximate the fit in the following way. We assume 
that M spec 2 should be compared to the differential luminosity 
(^) between the two fitted temperatures kT^^ (To = Thigh) and 



have M\ , , = Mciassicai. Despite the fact, as noted above, kT lovl . We assume that M spec should be compared to § between 



that scaling M', { . by M500 very likely overcompensates the 
tailing-off of clusters due to selection effects, there seems to be 
a marked difference in the parent low-z and high-z (red open 
circles) samples. 

Of the 10 clusters in the highest redshift bin with A^ lassical 
greater than zero, 8 have ^Ciassicai/^ 500 > 10~ 13 yr , while of 
the 8 clusters in the lowest redshift bin with M'^ . ■ greater 
than zero, only 2 have *f lassical /M 5 oo > 10~ 13 yr . If the two 
subsamples are drawn from the same distribution, the probabil- 
ity, p, of any one cluster having ^ lassical /M5oo > 10~ 13 yr -1 
should be the same for both. The joint probability, P{p) - 
P(p;m\ < 2)P(p;ni2 > 8), of drawing one sample of 8 with 
no more than 2 clusters in this category and a second sam- 
ple of 10 with at least 8, is maximized for p = p max = 
0.5584, giving P(p max ) = P max < 0.00883. Repeating this 
calculation for all possible pairs of draws (always using the 
smaller value for the first of the pair), if the true probability of 
M'. ■ ^/Msoo > 10~ 13 yr~' is q, the probabilities of each pair 
for which P max < 0.00883, can be summed to compute the like- 
lihood of this outcome. This is maximized for q = 0.5, giving 
0.029 for the likelihood. Thus, the chance that both samples are 
drawn from the same population is no more than 2.9%. 

We also used the Kolmogorov-Smirnov (K-S) test to distin- 
guish between the two samples. Even though it is well-known 
that the K-S method exhibits poor sensitivity to the deviations 
that occur in the tails of any given two distributions, we find the 
probability of the null hypothesis, i. e. the probability that the 
low-z and the high-z clusters are drawn from the same distribu- 
tion, is only ~ 2.4%. From this we can conclude that the high-z 



kT\ 0VI and zero. Our argument for using kT\ ovl and not kT^h is 
that the small values of M spec (compared to M c i ass icai and espe- 
cially M spec 2) are due to the lack of lines for the coolest gas and 
therefore it gives a better estimation of the true mass deposition 
between kT\ ovl and zero than between £7hi g h and zero. The frac- 
tion of observed gas is M spe c2/M c i assica i and M SDec /M c i a S sic a i re - 
spectively. Fig.[13]is similar to Fig. 7 in Peterson et al.l ( 2003 ). 
For each cluster with measurements at both points (i.e. neither 
was an upper-limit) we did a fit to Eqn. [24] Our rough estima- 
tion of 5, based on the weighted average value of the fits, is ~ 
2.0 (with a very large v ariance). This value is on the steep end 
of the values found by IPeterson et al.l (120031) . One may think 
this slight difference is due to the fact that we are working with 
an unbiased sample and their sample was picked to be SCCs (in 
general), so that M spec is not far from M c i ass j ca i. However, there 
does not seem to be a separation of the SCC clusters (blue) and 
the WCC clusters (black). The three outliers in the upper left 
are: NGC1399, A3266 and A2634, which as mentioned earlier 
have M spec <; M c i ass i ca i. This model requires further study with 
higher spectral resolution instruments before any definite con- 
clusion can be drawn, but it is reassuring that we find consis- 
tent results using a complete, unbiased sample of clusters that 
is more than four times larger. 

4.2. Temperature Profiles 

Perhaps two of the most striking plots from Fig. [6] are subplots 
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J and K, which show a sudden break at ~1 « ?1 Gyr. In fact, 
we find that for the most part only SCC clusters have a central 
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Fig. 13. The fra ction of differential lum inosity vs Tq/T v 
ilar to Fig. 7 in 



sim- 



Peterson et al.1 (120031) . The weighted average 



value of 5 (see text) for the clusters with measurements at both 
data points is ~2.0 and is shown as the thick red-line. The blue 
crosses and lines are SCC clusters and the black are WCC clus- 
ters. 



temperature drop. All 28 SCC clusters have a central temper- 
ature drop and of the eight non-SCC clusters with a central 
temperature drop (see Table g), four (A1650, A2065, A2142 
and A2589) have CCT < 2 li^ 1 Gyr, putting them on the bor- 
der between SCC and WCC clusters. Additionally two (A2142 
and A3667) have well-known cold fronts. The other three ap- 
pear to be merging clusters in which a cool core has surv ived: 
A0400 ([Hudson et al.ll2006l) . A0576 dDupke et al.ll2007bl) and 
A2256 dSun et alj|2002b . It is also clear from Fig.|6}K that the 
inner temperature profiles of these eight non-SCC clusters are, 
for the most part, flatter than those of the SCC clusters (slope>- 
0.2). We argue tha t bas ed on this fact, studies su ch as those of 
Burns et al] (120081) and lSanderson et al] (I2006al) . which define 
CC clusters as clusters with a central temperature drop, are ba- 
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sically defining CC clusters as clusters with CCT < 1 h lx ' Gyr 



(SCC clusters in our sample). 

Some authors claim the existence of a universal inner tem- 
perature profile for rel axed clusters (e.g. Allen et al. 200 lbt 



Sanderson et al. 2006bl), whil e others find no such universal- 



ity (e.g. IVikhlinin et al.l 120051) . A universal inner temperature 
profile suggests that clusters either have a fiat central tem- 
perature profile or a drop that scales with r v ; r and r v ; r . It is 
clear from Fig. HJ-L that we do not see such a universal in- 
ner temperature profile. Fig. [14] shows all 64 of our tempera- 
ture profiles out to 0.1 r v i r , scaled by r v i r and T v j T . It is clear 
from this figure why we do not see a bimodal distribution in 
Fig. HJ-L. There is a continuous range of slopes from a 20% 
increase above T v [ r (MKW8) to a decrease down to 0.2 T wu 
(A3526 - Centaurus Cl uster). We note that some authors (e.g. 
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Fig. 14. The inner temperature profiles of all 64 HIFLUGCS 
clusters scaled by virial temperature and radius. The error bars 
on the measurements have been omitted for clarity. There ap- 
pears to be a range of slopes rather than a bimodal distribution. 



radii when constructing temperature profiles. We, however, 
used double thermal models for the inner annuli in 14 of our 
SCC clusters since the addition of the second thermal compo- 
nent improved the fit significantly (> 99% confidence accord- 
ing to the F test). In these cases the lower temperature was 
used, causing a steeper temperature profile than if a single ther- 
mal model had been used. 

It is interesting to note that in Fig. |4}K, there does seem 
to be a bimodal distribution in Tq/T^ (albeit with a large dis- 
persion). This suggests that while the radius of the gas with 
T < r v ; r does not scale with cluster size, the depth of the 
drop does. SCC clusters have a central temperature of ~0.4 
r v ; r , whereas NCC and WCC clusters have To ~ r v ;,. As we 
discuss in the next section, it is unclear whether this cool gas 
is associated with the CF or the central galaxy. In the case of 
the former this indicates the coolest major component of gas 
allowed by the feedback mechanism. That is, significant quan- 
tities of gas do not cool below ~0.4 T v i r . In the case of the 
latter, it is consistent with observing that the mass of the bright- 
est cluster galaxy (BCG) scales with the cluster mass, so that 
T o/T v \ r is constant. T his model would support the simulations 
of Burns et al.l (120081) who find that the CC grows with the clus- 
ter. 



Sanderson et al.ll2006al) only use single thermal models at all 



4.2.1 . What Causes the Central Temperature Drop? 

An old but still ongoing controversy is whether the temperature 
drop in CC clusters is related to the CF or is simply a reflection 
of the potential well around the BCG. In the case of the for- 
mer, the central temperature should be reflected by a modified 
cooling flow model and should simply depend on the energy 
factors (energy radiated away, along with input energy from 
conduction and any feedback mechanism). In the latter case 
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the temperature should reflect the mass of the central galaxy 
and should be independent of any cooling flow model. 

One check for possible multi temperature components is the 
need for a double thermal model in the central annuli. In the 
central annulus (annuli) for 14 of 64 HIFLUGCS clusters (all 
SCC clusters), a double thermal fit was a significantly better fit 
than a single thermal model. It is possible, since these are pro- 
jected temperature profiles, this is simply due to many temper- 
ature components along the line of sight. Even in this case, the 
projected temperature components that dominate are the ones 
near the cluster center where the emission is peaked. Whether 
it is the projected spectra of outer annuli or truly the need for 
a double temperature fit to the central spectra, we take this as 
evidence of many temperature components in the region of the 
emission peak for these clusters. Thus the change in tempera- 
ture is rapid in these regions. To explicitly demonstrate that a 
double thermal model is consistent with more than two tem- 
perature components we ran simulations of spectra with four 
and eight thermal models with temperatures equally spaced be- 
tween 1 and 2 keV. Fitting these spectra indicates that a dou- 
ble thermal model can, at Chandra's, energy resolution, fully 
describe the plasma emission of the underlying four- or eight- 
temperature components and provides the statistically best fit to 
the spectrum (j^ 2 /dof ~10 For the eight temperature plasma 
the returned temperatures were: 1.80±0.08 keV and 1.26+0.03 
keV. Therefore we interpret annuli that require two thermal 
models to contain, at minimum two thermal components, with 
a low temperature component % To. 

Fig. [T5l shows the modified spectral mass deposition rate, 
M spsc 2 versus M c iassicai- The striking thing about this plot is that 
the two values are very often consistent with each other, es- 
pecially for the SCC clusters. Based on this result, one may 
naively believe that the gas cools to T\ ovl (the lowest temper- 
ature for the modified cooling flow model - see Sect. 12. 8t at 
the rate predicted by the classical cooling flow model, but does 
not cool at a significant rate below this temperature. If this is 
the case there are two possibilities: (1) the gas is in thermal 
equilibrium throughout this range of temperatures or (2) the 
gas is rapidly heated from r[ ow to Thigh, when it reaches ri ow . 
An alternate explanation is that with the spectral resolution of 
the Chandra ACIS, we are unable to distinguish between gas 
cooling from Thigh to T\ ow and a rapid drop of temperature at 
the center. That is, the observation that M spec 2 ~ M c i ass j ca i could 
simply be a coincidence. 

Investigating further, we plotted kT\ ovl versus kTo (Fig. |T6T> 
for the CC clusters. The blue points are SCC clusters and the 
black points are WCC clusters. If the gas cools from Thigh to 
Ji ow at a rate consistent with a classical cooling flow (and is 
somehow stopped at Ti ow ), then we would expect ri ow < To. 
Fig.fTBIshows this to be the case for SCC clusters, however no 
correlation is found for WCC clusters. As noted earlier, WCC 



18 For example, for eight thermal components, the average reduced 
X 2 for a single thermal fit was 1.46, 0.976 for two thermal compo- 
nents, and 0.984 for three thermal components. The simulations were 
done for a lOOks observation with the background and thermal nor- 
malization taken from A1795's central annulus. All thermal compo- 
nents were given the same input normalization. 
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Fig. 15. This plot is similar to Fig. [10] M spec 2 is the spectrally 
determined mass deposition rate between Thigh and T\ ow , where 
the lower temperature is left free (unlike for M spec ). The lines 
are the same as for Fig.[l0] The mass deposition rates measured 
by both methods generally yield consistent results. 



clusters show little or no temperature drop at the cluster center. 
However, when fitted with a modified spectral mass deposition 
model, there appear to be significant quantities of gas down to 
~l/3 TV;, . The simplest explanation is that M S pec2 is not the cor- 
rect model and at the energy resolution of the Chandra ACIS, 
it is difficult to distinguish between multitemperature compo- 
nents along the line-of-sight and a cooling flow model down 
to ~l/3 TVi,.. With the spectral resolution of the Chandra ACIS 
we cannot draw any strong conclusions, however the results 
for the SCC clusters argue for further investigation with higher 
spectral resolution instruments to see if this trend of Ti ow < To 
persists. 

In order to check for the correlation of the central clus- 
ter gas with the central galaxy, we compared the position of 
the brightest cluster galaxy (BCG) to the position of the emis- 
sion peak. We identified the BCGs from visual inspection of 
the R-band images of the Second-Epoch Digitized Sky survey 
(DSS2, see e.g. http : //archive . stsci . edu/dss) and then 
extracted ma gnitud es and redshifts from NED, HvperLedJ^I 



( Paturel et al. 2003b or the c ompilation by Andernach & Tago 
(see lAndernach et al.l 120051 for a description). We measured 
the position of the BCG by fitting a bidimensional Gaussian 
on its image, using the NRAO program 

FITSviev@ 

The typi- 
cal uncertainty was ~075. Visual inspection of these images al- 
lowed the identification of usually one, sometimes two or three 
BCG candidates. In the cases in which there was more than 
one BCG candidate, we used the candidate closest to the X- 



19 http://leda.univ-lyonl.fr/ 

20 http://www.nrao.edu/software/fitsview; 
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Fig. 16. The lower temperature in the modified cooling flow 
model (kTi ov/ ) versus the cluster central temperature (kTo). Blue 
points are SCC clusters, while black points are WCC clusters. 
In general Ti ow < To, which is expected if there are multi- 
temperature components between Thigh and T\ ov/ . The fact that 
the WCC clusters show a greater discrepancy between ri ow 
and To suggests that at the energy resolution of Chandra's 
ACIS it is difficult to distinguish between the modified cooling 
flow model and multitemperature components along the line of 
sight. The three cases where ri ow is unconstrained, M spec 2 is 
consistent with zero. 



ray peak. Fig.[17]shows the histogram of distance between the 
X-ray peak and BCG. The histogram is color coded: blue for 
SCC clusters, black for WCC clusters and red for NCC clusters. 
The fact that all SCC clusters have a BCG at the center seems 
to support the idea that the central temperature drop is related 
to the presence of the galaxy. On the other hand, 78% of all 
HIFLUGCS clusters have a BCG within 12 kpc of their X- 
ray peak (including 61% and 50% of WCC and NCC clusters 
respectively) and ~88% have a BCG within 50 hz\ kpc. There 
is a clear discrepancy between the number of WCC clusters 
(92%) and NCC clusters (61%) that have a BCG at the X-ray 
peak (< 50 hz\ kpc) and those that have cool gas at their center: 
33% and 11% respectively. Also as noted earlier the tempera- 
ture drop in these clusters is generally smaller than it is for SCC 
clusters. 

Here, we do not present this as evidence that temperature 
drops are or are not related to the potential of the BCG. On the 
one hand, all but one cluster (A2256) with a central temperature 
drop have a BCG cospatial with the X-ray peak (< 12 hz\ kpc), 
suggesting that the temperature drop is related to the BCG. 
On the other hand, there are many clusters with a BCG at the 
center and no central temperature drop. It appears that regard- 
less of the mechanism that causes the temperature of the gas 
in the central region to be below r v ; v , it is easier to heat the 
gas (to ~ 7\,j r ) in the central region than to permanently sep- 
arate it from the dark matter potential well (i.e. the BCG). 
Simulations suggest that cluster mergers can destroy cooling 



100 , 
Dist between BCG and EP ( kpc) 

Fig. 17. Histogram of projected separation of the BCG and X- 
ray emission peak. The colors represent: blue - SCC clusters, 
black - WCC clusters and red - NCC clusters. Most of the clus- 
ters (~78%) have a BCG within 12 hz\ kpc of the X-ray peak. 
100% of the SCC clusters, 61% of the WCC clusters and 50% 
of the NCC clusters have a BCG within 12 tc\ kpc of the X-ray 
peak. The uncertainty in the X-ray peaks is generally 4", which 
corresponds to ~4 hz}. kpc for our median redshift. 



flows (e.g. iBurns et al.ll2008. and references therein). It is not 
clear though that such heating need be done by something as 
energetic as a major merger. For instance, A1650 has an almost 
flat tem perature profile and y et no evidence for a recent major 



merger (IDonahue et al ] |2005|K Other examples of such clusters 
are WCC clusters A2244 (IDonahue et all Eooll) . A1651 and 



A 1060. 

To explore this further we compared the central stellar ve- 
locity dispersion of the BCG to the central temperature. The 
BCG central stellar velocity dispersions were collected from 
HyperLeda, using mean reported values and standard deviation 
of the values for the uncertainty. In a few cases, if only one 
measurement was reported, we used the reported error on the 
measurement. 16 of our 64 HIFLUGCS clusters (including 8 
SCC clusters) did not have any data available. Fig. 1X81 shows 
BCG central stellar velocity dispersion versus kTo- The points 
are color coded: SCC-blue, WCC-black and NCC-red. Clusters 
with a central temperature decline are marked with a circle and 
clusters without a central temperature drop are marked with a 
square. The 16 clusters with no data available are omitted, but 
have the same general distribution in kTo as the 48 plotted clus- 
ters. Even when looking at just the SCC clusters, there seems to 
be only a weak correlation between the BCG's central velocity 
dispersion and kTo- 

Assuming that our above conjecture is correct, that the gas 
is relatively easy to heat without removing it from the cen- 
ter of the potential, then it appears that the temperature of the 
central gas is not generally determined by purely gravitational 
processes in the local potential. That is, some heating mech- 
anism heats the gas either slightly above the local virial tem- 
perature in SCC clusters or up to ~ T vu in WCC and NCC 
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Fig. 18. BCG's central velocity dispersion versus kTo. Blue 
points are SCC clusters, black points are WCC clusters and red 
points are NCC clusters. Circles signify clusters with a central 
temperature drop and squares signify clusters without a central 
temperature drop. The 16 clusters with no BCG velocity disper- 
sion data available span the same range. There appears to be no 
correlation between the BCG's central velocity dispersion and 
kT . 



clusters. We emphasize that the mechanisms do not have to be 
the same in all cases (e.g. possible AGN outbursts or sloshing 
for SCC and WCC borderline cases and mergers for NCC clus- 
ters). This would explain the large scatter in Fig. Q~8] It would 
also predict that the central gas in clusters such as A 1650 will 
eventually return to the local virial temperature in the center. 
In such a case, in order to maintain pressure equilibrium with 
the surrounding gas, the density would increase dramatically, 
shortening the CCT and returning the cluster to the SCC sub- 
sample. It is not clear in this model which clusters would be 
capable of returning to their SCC status and which ones not. 
Pe rhaps the early maj or mergers of NCC clusters, as proposed 
bv lBurns et alj (120081) . prevent the formation of the large cen- 
tral galaxy needed to reform the CC. Perhaps once the entropy 
is increased enough, the gas cannot cool to the local virial tem- 
perature in the time available. Even when suppressed by mag- 
netic fields, thermal conduction is still likely to increase rapidly 
with temperatu re. Given that its effect is alr eady argued to be 
significant ("e.g. lZakamska & Narayanll2003l) . it may well play 
a role in preventing cool cores from being re-established after 
disruption by mergers in hot clusters. We suggest that simu- 
lations of the formation of CCs considering different central 
potentials could help understand which ones can re-form and 
which cannot and what physical effects are important in deter- 
mining the fate of the central gas. 



4.3. BCG Separation and Mergers 

Several works show a special relationship between the cooling 
activity in cluster cores and the brightest cluster galaxies lo- 
cated within a certain projected distance from the X-ray peak, 



typically 50 /j 71 ' kpc, to the X-ray peak (see lMittal et al. I I2009L 
for more details). Given that only eight of our 64 clusters have a 
significant (> 50 hz} kpc) projected separation between the X- 
ray peak and BCG suggests that separating the BCG from the 
gas is more difficult than simply heating the gas and/or disrupt- 
ing the cooling flow. Looking more deeply, the eight clusters 
with BCG-X-ray peak separations >50 tCj\ kpc ordered from 
largest to smallest are: (1) A3376 - -939 kpc, (2) A0754 



-714 h-\ kpc, (3) A1367 
-642 h-\ kpc, (5) A3667 - 
428 h-{ kpc, (7) A2256 



-666 h-\ kpc, (4) A1736 



-155 h-{ 
110 h~l 



kpc, (6) A2163 - 
, kpc and (8) A2255 - 
72 hjl kpco All eight of these clusters have been identified 
as merging clusters, A3667 is the only CC cluster (and it is a 
borderline CC/NCC cluster) and all except A 1736 have been 
identified as having a radio halo and/or relic(s) (see Appendix- 
[C]for individual references). The appearance of such diffuse, 
Mpc scale non-thermal emission is thought to be powered by 
major mergers. Of the other 56 clusters, only eight (15%) 
have been identified as possibly containing diffuse non-thermal 
emission on large scales: A1656, A3562, A85, A133, A401, 
A2152, A4038 and MKW8. Of these eight, for two (A401 
and A 2142) the detections seem unlikely dGiovannini & Feretti 
20001) . three (A85, A133 and A4038) are small scale (a few 



tens of kpc) relic s associated with nearby radio galaxies (e.g 



Slee et al.ll2001al) and one (MKW8) as identified as a possible 



relic as se en in the VLA Low - Frequency Sky Survey (VLSS) at 
74 MHz dCohen et al.ll2007l) . iKempner et all d2004l) identified 
three classes of radio relics that were fundamentally different. 
We argue that in these four cases, the relics are not the large 
scale Gischt associated with mergers, but AGN-relics. This 
leaves only two (4%) unambiguous detections of large scale 
diffuse radio emission: the Coma cluster and A3562. Based on 
the large discrepancy between the number of clusters with large 
scale (~Mpc) radio structure that have a large BCG-EP separa- 
tion and those that do not, we argue that a large separation of 
the EP and BCG is a very good indication a major merger and 
therefore could be applied as a useful method for discovering 
radio halos and radio relics-G/sc/if. 

4.4. WCC Clusters 

The WCC clusters are an interesting set because they seem 
to occupy a transition between NCC and SCC clusters. 
They are defined as having short to moderate CCTs (1.0- 

— 1/2 

7.7 « 71 ' Gyr), and generally have flat or shallow central tem- 
perature drops and a central entropy that is enhanced com- 
pared to their SCC cousins. Due to their flat or shallow tem- 
perature profile, they are mostly classified as NCC clusters in 
studies that de termine CC cluste rs by a central temperature 
drop. Recently iBurns et al. I ri2008l) suggested that CC clusters 
were not necessar ily more relaxed than NCC clusters. Since 
Burns et al.ld2008l) define CC/NCC classification based on tem- 
perature drop, they would classify the WCC clusters as NCC 
clusters. Could these WCC clusters be the relaxed NCCs? 



21 Additionally the peculiar velocity of A2255's BCG is much larger 
than that of any other cluster (greater than A2255's velocity disper- 
sion). 
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Donahue et al.1 (120051) studied two radio quiet, relaxed WCC 
clusters - A 1650 and A2244. They conclude that these objects 
are CC clusters that have had a major AGN outburst, raising 
the central entropy and temporarily disrupting the CF. 

In an effort to distinguish between these two possibilities 
we examined the X-ray morphology of the WCC and NCC 
clusters. We divided these 36 clusters into four types: (I) re- 

— 1/2 

laxed with short CCT (< 3 h 7l Gyr), (II) relaxed with mod- 
erate/long CCT (> 3 h~l /2 Gyr), (III) disturbed with short CCT 
(< 3 h'l 12 Gyr), (IV) disturbed with moderate/long CCT (> 
3 h^ 2 Gyr). We classify relaxed versus disturbed by visual in- 
spection. We define a relaxed cluster as having: (1) round or 
elliptical isophotes, (2) the EP at the center of the isophotes 
and (3) little or no offset between the centers of the different 
levels of isophotes (little or no sloshing). We note, based on 
our classification scheme, that all our NCC clusters will be of 
type II or IV. Our physical interpretation of the four types are: 
(I) CCs that have been disrupted, but will re-form - their CCTs 
are too short to be relaxed NCCs, (II) clusters that have not had 
a recent merger but have had their CC severely disrupted (or 
destroyed) by a previous major merger (III) clusters that are 
merging with the core surviving (although it may be destroyed 
in the future) and (IV) merging clusters in which there was no 
cool core or the core has been destroyed. 

In order to check our categorization based on visual inspec- 
tion, we plotted the distance between the EP and the X- ray 
emission weighted center (EWC) (see Hudson et al.ll2006l for 
details on determining the EWC). This method gives a simple 
way to classify clusters, since disturbed clusters will generally 
have a larger distance between their EP and EWC than a re- 
laxed clusters. Fig. [19] shows the results for the 36 WCC and 
NCC clusters. Black points indicate the clusters identified as 
being relaxed and the red points indicate the disturbed clusters. 
Fig. [19] confirms that, in general, the disturbed clusters have a 
larger separation between the EP and EWC than the relaxed 
clusters. 

We find seven clusters of type I, six clusters of type II, five 
cl usters of type III an d 1 8 clusters of type IV. By the definition 
of iBurns et al.l ( 120081) . these 36 clusters are clusters that will 
never form a CC. We argue, however, given the short cooling 
time of the clusters of type I, it is unlikely that their CCT will 

—1/2 

never drop below 1 /j 71 Gyr, which in turn implies they will 
also have a central temperature drop. 39% of the WCC clusters 
are of this type and are perhaps a special type of SCC cluster 
that experienced an anomalous event that has temporarily dis- 
rupted the strong cool core (e.g. raised the central temperature, 
entropy and cooing time). 

We acknowledge that a more sophisticated method should 
be employed when determining whether a cluster is relaxed or 
not. However based on preliminary results, unless WCC clus- 
ters (especially those with short CCT) are fundamentally differ- 
ent from SCC clusters, some of them will become SCC clus- 
ters. That is, unless there is a process that keeps their CCT > 

— 1/2 

1 h lx Gyr and their temperature-profile flat, the gas in the core 
will cool below T vn - and the CCT will drop below unity. 
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Fig. 19. Distance between the X-ray emission peak and emis- 
sion weighted center for the 36 WCC and NCC clusters. The 
black points are clusters that appear relaxed and the red points 
are clusters that appear disturbed. The vertical dashed line sep- 
arates clusters with short CCT and long CCT. The dash-dot line 
divides the clusters between WCC and NCC. It is clear that in 
general the separation of the EP and EWC is larger for the ap- 
parently disturbed clusters. The error on the separation is 4", 
which is the typical size of the smoothing kernel used when 
determining the EP. 



5. Conclusions 

We provide the most detailed systematic view into X-ray cores 
of galaxy clusters to date. We find that the best method to de- 
termine whether a cluster is a cool-core (CC) cluster is with the 
central cooling time (CCT). We divide clusters into three types: 
strong cool-core (SCC), weak cool-core (WCC or transition) 
and non-cool-core clusters. SCC are defined as having very 

—1/2 

short CCT (< 1 h 7l 1 Gyr) and are characterized by low central 
entropy (<30 luV^ keV cm 2 ), systematic central temperature 



drops (with 7o/r v ; r ~0.4) and a brightest cluster galaxy (BCG) 
at the X-ray peak. WCC or transition clusters are defined as 
having moderate CCT ( CCT between 1 - 7.7 h^ /2 Gyr) and 

— 1/3 

are characterized as having an elevated entropy (>30 /j 71 keV 
cm 2 ), flat or slightly decreasing central temperature profiles 
and having a BCG at or near (<50 hz\ kpc) the X-ray peak. 

— 1/2 

NCC clusters are defined as having long CCT (>7.7 h ll Gyr), 

— 1/3 

and are characterized by large central entropies (Kq > 110 « 71 
ke V cm 2 ) and temperature profiles that generally are flat or rise 
towards the center. 

Based on the above classification, our main conclusions 

are: 

1. In our flux-limited statistically complete sample, we find 
72% of the clusters are CC clusters with 44% of the clusters 
being SCC clusters and 28% being WCC clusters. 

2. For intermediate redshift clusters (where radii as small as 
0.4% rsoo cannot be resolved), we find that Kbias and 
scaled Lx are the best proxies for CCT. In general A^ias 
is better, but needs to be standardized before it is used. We 
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TC J «e«H rami with a possible 
value of / = 8.6 x 10 65 h~\ cm" 3 . 

3. For high redshift clusters (clusters with very few counts), 
we find that cuspiness is the most useful proxy. However, it 
suffers from a large scatter in the relation to CCT, so precise 
categorizations can be difficult. Preliminary results suggest 
that a concentra tion parameter (£(< 40 kpc)/Z(< 400kpc)), 
as suggested by lSantos et al.l (120081) . may serve as the best 
proxy for distant clusters. 

4. Dividing our representative sample into four redshift bins, 
we find it unlikely that the discrepancy between spectral 
mass deposition rate (M spec ) and classical mass deposition 
rate (M c iassicai) is due to the fact that CCs formed very re- 
cently. The probability that the lowest- and the highest- 
redshift clusters with non-zero forward-evolved mass de- 
position rates come from the same population is extremely 
low (< 3%). This requires the discrepancy to be explained 
by some heating method. 

5. There is no evidence found in our work for a universal cen- 
tral temperature profile as claimed previously for smaller 
samples. This suggests that the radius of the cool gas is not 
universal. As seen in Table Q] the radius within which the 
temperature falls below T vlI differs from cluster to cluster. 

6. We find the majority (-78%) of HIFLUGCS clusters, in- 
cluding 100% of the SCC clusters, 61% of the WCC clus- 
ters and 50% of the NCC clusters, have a BCG within 
12 \c\ kpc of the X-ray emission peak. This number in- 
creases to 88% for a distance of 50 hz\ kpc. We find that 
seven out of eight clusters with the BCG-EP separation > 
50 tCj\ kpc have Mpc scale radio emission (halo or relic - 
Gischt) versus only two of the remaining 56 with a BCG 
very close to the EP. 

7. There is a weak correlation between the SCC central tem- 
perature and the cooling of the gas as predicted by the cool- 
ing flow model. We also find no correlation between the 
SCC central temperature and the central velocity disper- 
sion of the BCG. We interpret this, along with the result 
that so many clusters have a BCG at their peak, as indi- 
cating that the central temperature of the gas is influenced 
by heating, which occurs differently in different clusters. 
Therefore, the central temperature is not simply reflective 
of the central potential well or the expected cooling rate 
from Thigh to ri ow , but requires more complex physics. 

8. We find that <39% of the WCC clusters which are relaxed 

— 1/2 

also have relatively short cooling times (f coo i < 3 /? 71 ). 
We argue that these clusters are similar to SCC clusters but 
have had an event that temporarily disrupted the cool core, 
raising the core temperature, entropy and cooling time. 

Results obtained here on cluster cool-core properties and 
their correlations, using a complete well-controlled sample 
with simple selection criteria, can be taken as benchmark 
for next generation cosmological hydrodynamical simulations 
aiming at describing detailed properties of cluster cores. 
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Appendix A: Calculating Central Density for a Double fi model 

Starting from the definition of the normalization of the APEC model (IMewe et al.l I1985L Il986t ISmith & Brickhouse 1200(1 
Smit h et al . 2001) and taking n{r) = n e with ( - —, (calculated individually, but generally ~1.2), 



N = 



i(r) 2 dV, 



io- 14 r 

I n( 

AnD A D L U 

where the terms are defined as in Eqn.|9] For a double /3 model the expression for n{r) is: 
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The unabsorbecQ surface brightness at a projected distance, x from the center over an energy range between E\ and E2 is 

A x (T,Z,E)dE r°° 
E(jc) = nirfdl, 



4tt(1 



(A.3) 



where r 2 = x 2 + I 2 and Ax{T,Z,E) is the emissivity function for a plasma of temperature T and metalicity Z at energy E. This 
can be rewritten in terms of «oi and «02 as: 



C 2 A x (T,Z,E)dE r °° 



2 ^ ;2\\-3ft 



A x (T,Z,E)dE 



00 



l + 



JC 2 + / 2SS-3^ 



4^(1 

Solving the integral gives the standard expression for the double /? model in terms of surface brightness: 
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where LI, and are as defined in Eqn.Q~2] Using this relation along with the fact that no = n(0) = ^Tn^, we find: 



"01 = 



and 



E12 LI2 2 
-n , 



S12 LI 2 + LI, 



LI, 



nl. 



02 S12 LI 2 + LI, "° 

Inserting these values into Eqn. IA.2l to find an expression for n(r) in terms of no, we get 
V2i2 LI 2 + LI, 

Inserting this expression of n{r) into Eqn. lA.ll and solving for no, we recover Eqn. 



(A.8) 



(A.9) 





< 1 \2> 


-3/3, 






£12 LI 2 




+ LI, 

















(A. 10) 



22 In our case Eqn. lA.ll alreadv takes into account the absorption and any absorption in the subsequent calculations cancels. 

23 Note, we have explicitly assumed Ax(T,Z,E) is the same for both components. In the case it is not, £12 can be redefined as 
A X (T 2 , Z 2 , £*)S i /Ax(Ti , Z, , £)S 2 and the calculations follow identically. 
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Appendix B: ^bias Calculations 

From the definition of surface brightness (Eqn. IA.3l i. a cluster at redshift z, of a region with an angular radius x, has an integrated 
surface brightness (or Flux T) between energies Ey and Ey. 

I A x (T,Z,E)dE r°° rx 

r=J ^F, 3 n e n H xdxdl, (B.l) 

2(i + zr J-ooJo 

where n e is the electron density, nn is the proton density, Ax(T, Z, E) is the emissivity function as defined in Eqn. |A.3l To remove 
the redshift dependence of the projected region size, we convert the projected region of angular radius x to a cylindrical region 
of physical radius R, such that R as xD A {z). Eqn. lB.ll becomes: 

C A x (T,Z,E)dE 

f ° M , M u,)- W (B2) 

where D A is angular diameter distance, D L is the luminosity distance and I(R) is defined as in Eqn.[20] Therefore the total counts 
C collected by a telescope for an observation of length f t, s , in an energy band from E\ to Ei, of a cylindrical region of physical 
radius R is: 

ff 2 a(E)A x (T, Z, E)A eS (E)dE 

C = t obs — = I(R), (B.3) 

4nD A D L (1+z) 2 

where a(E) and A e g(E) are the absorption from Galactic hydrogen and the effective area of the telescope at energy E, respectively. 
We can calculate J E ' a(E)Ax(T, Z, E)A e ff(E)dE for an absorbed thermal model using XSPEC with an appropriate ARF and 
RMF Specifically, since normalization N oc CR = C/f Q b s , XSPEC can be used to find the constant of proportionality k. From the 
definition of N (see Eqn. IA.ll i: 



JE 



2 a{E)A x (T, Z, E)A eS (E)dE = (B.4) 
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Using an on-axis Chandra ARF and RMF, we determined k (122.3 photons cm 5 s~') for an energy band from 0.5 - 7.0 keV, with 
Z = 0.25 solar, Nh — 2 x 10 20 cirT 2 and our median observation time (44 ks), redshift (0.047) and virial temperature (4.3 keV). 
Inserting our determined value of k into Eqn. lB.5l and solving for /, such that C = 10000 counts, yields / = 8.6 x 10 65 cirT 3 . 
Therefore using the criterion that our median observation would have /Tbias determined by circle with 10 000 counts, equivalent 
regions from other observations would have: 



JL = 17 910(-^)( 
cnts \ lOOks A 



(200 At. 1 Mpc) 



100 cnts cm 5 s 
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Appendix C: Notes On Individual Clusters 

C.1. A0085 



This cluster appears to have two subclumps, one near the center and one further to the south (IKempner et al ,l2002l) . In determining 



the temperature profile and global cluster temperature the latter w as excluded. This SCC cluster hosts a well-studied radio relic, 
which is close to but not connected to the central radio galaxy (e.g Slee et al1l2001bh . The central region of this cluster requires a 



double thermal model out to ~11" (~12 h-,\ kpc). 



C.2. A0119 



This is possibly a merging cluster, which shows elongation towards the northeast. The X-ray peak of this NCC cluster, which 
does not dominate the surface brightness, has a cD galaxy cospatial with it. T he cluster contains three wide-angle-tailed (WAT) 
radio galaxies which may be interacting with the ICM (e.g. Feretti et all 19991) . 
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C.3. A0133 



Centr al regions of this cluster show an east- west elongation. An in-depth study with XMM and Chandra bv lFuiita et al.l ([2004, 



20021) revealed an X-ray ton gue extending nor thwest. This SCC cluster hosts a radio relic, that is close to but not connected to 



the central radio galaxy (e.g Slee et al. 2001bl) . The central region of this cluster requires a double thermal model out to ~15' 



46 /i 7 j kpc). 



C.4. NGC0507 Group 

The overall X-ray spectrum of this group shows a suspicious hard tail. An additional powerlaw component was included in the 
overall temperature fit. It is possible that the hard tail is due to unresolved low mass X-ray binaries (LMXBs), although the 
central region (~L8 = 35.7 h^ 1 kpc) was removed and no evidence of a hard excess is seen in the spectra of the central annuli. 
It is possible that (given the redshift of NGC0507 z = 0.0165) the LMXBs are only strong enough to be measured in a large 
region and are insignificant compared to the group emission in the central region. On the other hand the powerlaw has a steep 
photon index (Fx = 2.3 - 3.0) that is not consistent with LMXBs, which usually have a photon index of Tx=1.6. The component 
has a total flux of ~4 x 10~ 12 ergs cm~ 2 s _I corresponding to a luminosity of 2 x 10 42 ergs s _1 (over 0.4 - 10.0 keV). It 
is also possible that the hard excess is related to an insufficiently subtracted particle background, which is visible in this cool 
cluster. Both models (additional particle background or powerlaw) give an identical overall temperature. We also note that a more 
detailed analysis of the residual background in outer cluster regions shows no residual particle background. The central region of 
this cluster requires a double thermal model out to ~62" (~20.6 kpc). 



C.5. A0262 

The spectral fits to the inner regions are poor ( ^ 2 /dof ~ 1 .4) even with a double thermal model. Using non-solar abundance ratios 
significantly improves the fit, but does not change the best-fit overall temperature. We therefore used solar ratios for simplicity. 
The central region of this cluster requires a double thermal model out to ~43" (~14 hz\ kpc). 



C.6. A0400 



This cluster hosts the double radio source 3C75 within its center and shows evidence of merging dHudson et al .12006]). As noted in 



Kalberla et al 



Hudson et al. I d2006l) . the hydrogen column density is higher than measured with the radio (Nh = 0.85 xlO cm 
20051) and therefore we left it free for all spectral fits. We find a hydrogen column density of Nh = 0.98-1.22 xlO 21 crrr 2 for our 



fit to the overall cluster. 



C.7. A0399 

This cluster is near to A401 and shows evidence of interaction with A401 (e.g. Sakelliou & Ponmanl 2004 ). The temperature 
profile of this cluster peaks at the X-ray center. 



C.8. A0401 

See also A0399. This cluster may host a radio halo dGiovannini et al. I ll999l) . We included an early observation (before 2001), 
since the later observation was offset, with the cluster center in the corner of a CCD. The BCG closest to the X-ray peak is 
~ 34 /i^j 1 kpc away, making it one of fourteen clusters with the BCG >12 hz\ kpc from the X-ray peak. 



C.9. A3112 



Although the background flaring seen in some observations was removed, the effect seems to have broadened a fluorescence 
line. This can be seen in the fit to the overall cluster spectrum (at ~7.5 keV). This effe ct seems simply to ma ke the fit poor 
(X 2 /dof ~1.6), but it does not affect the best-fit values whether the line is removed or not. iTakizawa et al. I d2003h first presented 
the Chandra data of A3 1 12, interpreting the radio active central cD galaxy as interacting with the ICM. iBonamente et al.l (120071) 
claim a soft excess and hard excess in this cluster that may be related with the central radio active BCG. We do not see a similar 
effect, however we do not separately fit the l'-2'5 annulus that lBonamente et al.l d2007l) fit. We do confirm that the l'-2'5 annulus 
is isothermal in our AT-profile so that their result is not due to a temperature fluctuations in the cluster. This SCC cluster is one 
of sixteen clusters for which no data exist for the BCG central velocity dispersion. 
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C.10. NGC1399 Group/Fornax Cluster 

This nearby SCC cluster has two X-ray peaks cospatial with NGC1399 and NGC1404. The X-ray peak is taken to be cospatial 
with the BCG NGC1399. The peak on NGC1404 was removed for spatial and spectral analysis. Fornax appears to be an outlier 
in six of the plots of parameters versus CCT in which it has an anomalously low value for its CCT These parameters are: (1) So, 
(2) rjrsoo, (3) /Tbias, W co °li n g radius, (5) M c i ass i ca i/M5oo and (6) M gas /M5oo- Additionally it is the only SCC cluster in which 
^classical ^ M spsx . One possible explanation is that NGC1404 is about 0.04 rsoo from the X-ray peak and due to the extended 
emission around it, the surface brightness profile severely flattens. The extrapolated outer profile therefore overestimates the 
projected gas lowering the central density (and altering associated parameters). To check how much this influenced the Fornax 
cluster as an outlier, we fit only the central part of the surface brightness that could be fit well to a single /?-model. This model 
most likely underestimates the projected gas, thereby providing the largest possible values for central density. In the case of 
Fig. [10] this model raises M c i ass i ca i to 0.75+0.04 M Q yr' 1 making it larger than M spectra i = 0.52+0.02 M Q yr' 1 . We emphasize 
that this result overestimates M c i ass icai and in any case M spec » M c iassicai ! making it an odd SCC cluster. In all other cases the 

— 1/2 

Fornax cluster remains an outlier. In the case of the CCT the value falls to -0.4 Gyr moving the Fornax cluster to the left 
in Fig[6] (1) So, not surprisingly, remains almost the same, ~0.21 x 10' 6 photons cm' 2 s' 1 arcsec' 2 . The shift to the left in this 
case makes the Fornax cluster even more of an outlier. (2) r c /r5oo, likewise, remains almost identical, ~4 x 10' 4 . Unlike, So, the 
shift to the left in this case makes the Fornax cluster more similar to other SCC clusters. (3) Kbias is, of course, unaffected by 
the density model, however moving the Fornax cluster to the left makes it more consistent with the other SCC clusters. (4) The 
cooling radius remains almost identical, increasing to 0.03 rsoo. Additionally moving the Fornax cluster to the left makes it even 
more of an outlier. (5) M c i ass i ca i/M5oo increases slightly to ~ 1.8 x 10' 14 yr' 1 . Moving the Fornax cluster to the left makes the 
trend of decreasing Mctessicai/Msoo w i tn CCT worse, but it is more consistent with other groups with short CCT. (6) Mgas/Msoo, 
is raised slightly, ~0.2 x 10' 3 . Similar to M c i ass i ca i/M5oo, moving the Fornax cluster to the left makes the trend with CCT worse, 
but makes it more consistent with the other groups. We emphasize that the Fornax cluster is still an outlier in all six cases, and 
that these values are the other extrema, with the true value somewhere between these and where they are plotted in Fig [6] The 
physical interpretation is th at it is possible that Fornax is in the process of forming a cool core (i.e it has a nascent core in the 



terms of lBurns et al.l (12008b ) and therefore is dynamically different from the other SCC clusters that have well-established cores. 
The central region of this cluster requires a double thermal model out to ~170" (~ 15.9 kzl kpc). 

C.11. 2A0335+096 

This cluster, along with A0478 and NGC1550, has a sig nificantly higher hydrogen column density than measured at radio 



wavelengths (Nh = 1 .68 xlO 21 cm' 2 iKalberla et al.l 120051) . We fit all spectra with the column density free. For the fit to the 
overall spectrum we find Nh = 235-2 Al xlO 21 cm' 2 . This cluster has two major galaxies near the X-ray peak, which resides 
between the two of them (~ 10 h~\ kpc from the closest). Of the 16 clusters in which no information about the BCG central 



velocity dispersion is available, this cluster has the shortest CCT. The central region of this cluster requires a double thermal 

-1 
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model out to ~38" (~26 h 7 } kpc) 



C.12. IIIZw54 

IIIZw54 is a pair of galaxies near the center of a poor galaxy group. We used a 6" kernel when smoothing the image in order 
to determine the emission peak. This cluster does not have a bright core, although it appears to be quite round and relaxed. The 
brighter of two galaxies in the galaxy pair IIIZw54 (a cD galaxy) is cospatial with the X-ray peak. 

C.13. A3158 



Lok as et al. I d2006h report A3 158 as a relaxed cluster based on the velocity dispersion of the galaxies. The X-ray emission appears 



to be elliptical and there are two cDs near the cluster center, one of which lies at the X-ray peak. This cluster definitely does not 
have a bright core, with a central density of only ~5 x 10' 3 cm' 3 . The temperature profile peaks in the center at ~5.7 keV. 

C.14. A0478 

This cluster, along with 2A03 35+096 and NGC15 50, has a significantly higher hydrogen column density than measured in the 
radio (Nh — 1-64 xlO cm'" Kalberla et al ] |2005l) . We fit all spectra with the column dens ity free. Our fit to the o verall cluster 



yields a column density of (Nh = 2.89 - 2.96 xlO 21 cm' 2 ), consistent with the value found bv lSanderson et al.l (12005b . This cluster 
has the highest spectral mass deposition rate of any HIFLUGCS cluster, making it an i deal candidate for a grating observation. 
Unfortunately the RGS data from a long XMM-Newton exposure was virtually unusable ( de Plaa et alj|2004 ). This SCC cluster is 



also one of sixteen clusters in which no data for the BCG central velocity dispersion are available. 
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C.15. NGC1550 Group 

This cluster, along with 2A0335+096 and A04 78, has a significantly higher hydrogen column density than measured in the radio 
(Nh = 0.981 xlO cm 2 lKalberlaetaDl2005h . We fit all spectra with the column density free. Our fit to the overall cluster yields 
a column density of Nh = 1-34 - 1.41 xlO 21 cirT 2 . The column density appears to peak towards the center of the cluster, having 
a value of Nh = 1.9 - 2.2 xlO 21 cirT 2 in the innermost annulus. 



C.16. EXO0422-086/RBS 0540 

The short observation of this SCC cluster indicates a round, centrally peaked clust er with a mod e rate c entral temperature drop. 
The BL Lac object EXO 0423.4-0840 at the center of this cluster was studied by iBelsole et alj d2005l) . This is one of sixteen 
clusters for which no data about the BCG central velocity dispersion are available. 



C.17. A3266 

This cluster has a very low background scaling factor; therefore an additional unfolded powerlaw component was included in the 
spectral fits to account for any residual particle background. Reading in the background as a corfile (i.e. a second background 
component with an adj ustable scaling factor), the o verall best-fit temperatur e is found to be consistent with our result including 
an unfolded powerlaw. iHenriksen & Tittleyl(l2002l) : lFinoguenov et al.l(l2006h presented detailed analyses of this merging system. 



C.18. A0496 

The high abundances in the central region of this cluster are better fit with a VAPEC model, however since this did not change the 
best-fit values of temperature, solar ratios were used for simplicity. A double thermal model greatly improved the fits to spectra in 
annuli out to Q'.3. However, the high temperature component for annuli between Of 18 and 0'3 is » kT wu - ( kT ~ 8 keV). Although 
this may be evidence of very hot gas near the cluste r core, the investigati on is beyond the scope of this paper. Therefore for annuli 
between Of 1 8 and 03, we used a single thermal fit. iDupke et alj d2007al) studied the longest of the three Chandra observations of 
this cluster in depth. They argue that there is a cold front at the center of this cluster, which is caused by an off-center passage of 
a smaller dark matter halo. 



C.19. A3376 

This cluster was fit with an unfolded powerlaw component to account for possible low-level flares in both obse r vation s. This 
cluster appears highly disturbed in the X-ray with a strong elongation along the east- west directio n. Bagc hi et al. I d2006l) report 
the existence of double relics, one to the east of the cluster center and one to the west. iNevalainen et al.l d2004l) found a diffuse, 
hard excess with the BeppoSAX PDS at 2.7 '<x significance. The BCG of this cluster is ~ 1 Mpc from the X-ray peak, the most 
distant of any cluster in the sam ple and one of eigh t clusters with a separation of >50 kpc. There is a radio galaxy with bent 



jets very close to the X-ray peak dMittal et al ]l2009h . Optically it is clearly much fainter than the BCG and is most likely an AGN 
that may have been activated by the merger. The jets are bent in the opposite direction to the elongation of the cluster, possibly 
bent from ICM ram pressure. 



C.20. A3391 

The short observation of this NCC cluster shows an elliptical shaped ICM with a BCG cospatial with the emission peak. 
Tittley & Henriksenl ( 2001 ) discovered a filament between A3391 and the nearby cluster A3395. 



C.21. A3395S 



This c luster is very close to and may be interacting with A3395e. A3395e was excluded from all extended analysis. Donnelly et al 
(2001) claim A3395s and A3395e are near first core passage. 



C.22. A0576 

Kempner & Davidl (120041) originally presented an analysis of the Chandra data. lDupke et al. (2007b) presented a detailed analysis 



of the XMM-Newton and Chandra data suggesting that it is a line-of-sight merger. The X-ray image seems s omewhat perturbed 
with e lliptical isophotes with alternating NW-SE shifted centers, reminiscent of sloshing, already noted by iKempner & David 



(2004). The BCG is ~24 h 1{ kpc from the X-ray peak, making it one of fourteen clusters with the separation >12 h lx kpc. 



There is, however, a slightly fainter galaxy closer (<12 /i 7l ' kpc) to the X-ray peak that is radio active, whereas the BCG is not. 
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The peculiar velocity of the BCG is one of five clusters that is more than 50% of the velocity dispersion. This WCC cluster is 
one of the three WCC/NCC clusters with CCT»1 h^ /2 Gyr (i.e. not on the border between SCC and WCC) and a systematic 
temperature decrease at the cluster center. 



C.23. A0754 



This irregularly shaped cluster hosts a halo and relic (K assim et al.ll2.Q01l) . iHenrv et al.l ( 120041) presented a detailed analysis of 
this merging system using the XMM-Newton observation. Only the pre-2001 Chandra observation is used, since it was the only 
one that contained the cluster core. More recent observation have been made but do not cover the cluster center and therefore are 
not useful for core studies. The BCG for this cluster is ~714 hz\ kpc away from the X-ray peak, making it one of eight clusters 
where this separation is >50 tCj\ kpc. 



C.24. A0780/Hydra-A Cluster 

This cluster is known to have a massive central AGN outburst (Nulsenetal ]|2005l). 



C.25. A1060 



This WCC cluster is also known as the Hydra cluster. ISato et al.l (120071) recently presented a Suzaku observation of this cluster. 
This cluster has two brigh t galaxies near the cor e, one of which is cospatial with X-ray peak. Both galaxies have a clearly visible 
diffuse X-ray component ( lYamasaki et al 1120021) . 



C.26. A1367 



Due to the short exposure time and lack of a bright core, we used a 12" kernel when sm oothing to determine the X-ray peak . 
This is a very well-studied merging cluster. This cluster h as an infalling starb urst group ( Sun et al. 2005 ; Cortese et al. 20061) 
and optical evidence suggests that t his is a merging system (ICortese et al.ll2004h . The X-ray image appears rather disturbed with 
several off-centere d bright sources. Sun & Vikhlininl (12005b studied the survival of galaxy coronae in this system. This cluster 
hosts a radio relic dGavazzi & Trinchierilll983l) . The BCG of this cluster is ~666 hz) kpc from the X-ray peak making it one of 



eight clusters where this separation is >50 h lx 
the cluster velocity dispersion. 



kpc. It is also one of five clusters where the BCG peculiar velocity is > 50% of 



C.27. MKW4 

A single thermal model is a poor fit to this high metalicity center. Although a second thermal model does provide an improvement, 
freeing the ratio of elements for a single thermal model provides the best-fit. Since freeing the abundance ratios does not change 
the overall best-fit temperatures of the annuli, solar ratios with a single thermal model were used for simplicity. 



C.28. ZwCI1215 

The short observation of this NCC cluster, shows a round cluster with no bright central peak and an elliptical BCG located at the 
X-ray emission peak. This is one of the sixteen clusters for which no data about the BCG central velocity dispersion are available. 
The BCG of this clusters is also ~18 hz\ kpc from the X-ray peak, making it one of fourteen clusters where this separation is 
>12 fCj\ kpc. 

C.29. NGC4636 Group 

This nearby group contains extended nonthermal emission in the central region extending out ~122" (~9. 19 hz\ kpc). The lumi- 
nosity of this emission (Lx ~10 40 ergs s~') is consistent with the expected unresolved LMXB population for NGC4636. In 
addition to a powerlaw component, the central region of this cluster requires a double thermal model out to ~35" (~2.6 tCj\ kpc). 



C.30. A3526/Centaurus Cluster 



This is a well-studied, prototypical CC cluster, with a central temperature drop (having the largest fractional drop, IcTq ~ 0.2 
&7W) and enhanced central metalicity. An arc-like X-ray feature n ear the center has been identified as most likely being a cold 
front associated with sloshing in the core ( I Sanders & Fabianll2002l) . The central region of this cluster requires a double thermal 
model out to -72" (-16.5 hr\ kpc). 
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C.31. A1644 



As with A0085, this SCC cluster shows evidence of merging, with the existence of a double X-ray peak. lReiprich et al.l d2004l) an- 
alyzed the XMM EPIC observation of this cluster. They found the flux of the northern (smaller) subclump is below the HIFLUGCS 
flux limit whereas the flux of the southern (larger) subclump is above the flux li mit. Therefore for pur poses of this analysis the 
smaller subclump was excluded from spatial and spectral analysis. Additionally iReiprich et al.l (I2004I) found evidence that the 
smaller sub-clump was being stripped as it passes through the ICM. This is one of sixteen clusters in which the central velocity 
dispersion of the BCG is unavailable. The central region of this cluster requires a double thermal model out to ~32" ( 
kpc). 



30 h-\ 



C. 32. A1650 

This CC cluster hosts a radio quiet cD galaxy dDonahue et al .l2005l) . lMittal et al.l d2009b find an upper-limit to the b olometric radio 



lumin osity of ~9 x 10 38 /i 71 2 ergs s l . The original short Chandra observation showed a flat temperature profile dDonahue et al 



2005). However, th e longer, mosai c ed ob servations show a slight temperature decrease in the central region. Due to the elevated 



entropy in the core. lDonahue et al.l (12005b concluded a major AGN outburst had disrupted the cooling flow. This cluster is one of 
four clusters on the border between SCC and WCC. Its CCT (~1.2 h~^ 2 Gyr) is slightly longer than the 1 Gyr cutoff. Moreover 
this cluster shows a central temperature decrease typical of SCC clusters. This is one of sixteen clusters in which the central 
velocity dispersion of the BCG is unavailable. 



C.33. A1651 

As with A 1650, Donahue et al. (2005) claim this is a radio q uiet CC cluster, howev er Mittal et al. ( 2009t) detect central radio 



TO 40 h'l ergs s _1 . 



Gon zalez et al. d2000l) fit the optical light out to ~670 h,L kpc 



emission with a bolometric luminosity of 
over one quarter of r v i r . The X-ray structure looks quite round and shows no evidence of external interaction. However, the X-ray 
peak does not dominate as much as in SCC clusters and there is no evidence of a central temperature drop. This WCC is one of 
sixteen clusters in which the central velocity dispersion of the BCG is unavailable. 



C.34. A1656/Coma Cluster 

This well-stud ied NCC cluster appears to be involved in a merger with a group. This cluster hosts the first detected radio halo 
dWillsonlll970l) . 



C.35. NGC5044 Group 

The spectra for the inner regions of this group are not well fit by a single thermal model (^ 2 /dof > 2). After trying several 
different models to fit the residuals, we found that the statically best model which is also physically motivated is a thermal model 
that allows oxygen, silicon, sulfur, and iron to vary from solar ratios and a powerlaw to account for a the clear hard tail (most 
likely due to LMXBs). We note that a double thermal model with the above elements not constrained at solar ratios provides 
the statistically best-fit (in the innermost annulus ^ 2 /dof = 1.04 vs.^ 2 /dof =1.13 for the thermal plus powerlaw model). In this 
model, however, the hotter thermal model has a temperature of kT = 1.4 - 3.0 which is hotter than any gas found in the outer 
annuli and the measured virial temperature (kT w \ r = 1.22*q|S). Unless there is a hot halo of gas extending from the center of the 
group out to ~16 h^l kpc, this model is unphysical. Finally, we note that adding a powerlaw to the second thermal model does not 
improve the fit and similar high temperatures are found for the hotter thermal component as for the simple two thermal model. 



C.36. A 1736 

This NCC cluster is a member of the Shapley Supercluster. Due to the short exposure time and lack of a bright core, the X-ray 
peak was found by smoothing the image with a ~10" kernel. The X-ray morphology shows an irregular shape with no well- 
defined core. A preliminary temperature map shows heating to the east and west of the emission peak with cool gas extending to 
the south. The BCG of this cluster is ~642 hz^ kpc from the peak, making it one of eight clusters where this separation is >50 h~\ 
kpc. This is the only cluster with a separation of >50 \c\ kpc that does not have any known associated radio halo or relic. 



C.37. A3558 



This WCC cluster is located in the core of the Shapley Supercluster. The observation was heavily flared, an d even after a 
conservative cleaning of the light curve there was evidence of some low-level flaring in the back-illuminated chips. iRossetti et al. 
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( 120071) presented the XMM and Chandra analysis of this cluster, concluding that it had a cool core that had survived a merger. 
We find evidence of a slight temperature drop in the core of this WCC cluster. 



C.38. A3562 

This WCC cluster is located in the core of the Shapley Superclus ter. The X-ray emission from this cluster app ears to be elongated 
along the northeast-southwest direction. iGiacintucci et al.l (120051) re port the detection of a rad io halo (also see lVenturi et al . 2000) 
and argue for a merger scenario between A3562 and SC 1329-313. iFinoguenov et al.l (120041) presented a detailed analysis of the 
XMM observation of this cluster. The BCG of this cluster is ~31 hz\ kpc from the X-ray peak, making it one of fourteen clusters 
where this separation is >12 }c\ kpc. However, the BCG is located in a chip gap, so the separation may simply be an instrumental 
effect (i.e. the peak on the BCG may not be detected due to the chip gap). The XMM observation also shows an offset between 
the X-ray peak and BCG but on a scale of only ~23 \C,\ kpc (Y. Y. Zhang, private communication). 



C.39. A3571 

This WCC cluster is a member of the Shapley Supercluster. 



C.40. A1795 



The core of this well-studie d SCC cluster has a lar ge filament seen in X-rays and Ha (jCrawford et al. 2005 ). Early core studies 



with Chandra were done bv lFabian etai]d200ll) and lMarkevitch et alJj200ll) . lFabian et alJd200ll) found a CCT of -0.4 h^ L Gyr, 



approximately the same age as they estimate for the filament. The difference between their measurement for CCT and our 
measurement is probably due to the different values of r coo i used to determine CCT. In order to keep consistency between clusters 
we determined the CCT at 0.004 rsoo, however at the red shift of A 1795 we a re able to determine the CCT at an even smaller 

— 1/2 I I j I 

radius which gives a CCT ~0. 5 h nx Gyr, cons i stent w ith lFabian et alJ (|2001). Moreover, their technique for determining CCT 
is sl ightly differen t from ours. iMarkevitch et alJ (1200 II) found a cold front in the core of A1795, which they attribute to sloshing 



gas. 

with our measurement of ~15 M yr~ 



— 1/2 

Oeg erle et al.l (120011) studied FUSE observations and found an upper limit for M spec (20 kpc) < 28 h 7l M q yr" 



consistent 



C.41. A3581 

This SCC cluster is a member of the Shapley Supercluster. It is the only Shapley Supercluster member from the HIFLUGCS 
sample that is an SCC cluster. It is also the most distant of these clusters from the center of the Shapely Su percluster. The central 
region of this cluster requires a double thermal model out to ~40" (~18 }c\ kpc). I Johnstone et al.l d2005l) analyzed the Chandra 
data from A3581. They find a point source coincident with the powerful radio source PKS 1404-267 at the cluster center. They 
find a central temperature drop to ~0.4 kT v i r at the cluster center, similar to our measurement of ~0.5 kT v i r . 



C.42. MKW8 

This NCC cluster shows little substructure in the X-ray image. The X-ray isophotes are elliptical with the major axis along the 
east-west direction. The isophotes seem to have a common center (i.e. no evidence of sloshing), however the X-ray peak appears 
to lie to the east of the center of the isophotes. There are two bright galaxies at the center of the cluster. The brighter of the 
two corresponds to the X-ray peak (which unfortunately falls in a chip gap). The second galaxy is to the east, corresponding to 
the direction of the elongation of the surface brightness. This cluster shows a possible radio relic at 74 MHz in the VLA Low- 
Frequency Sky Survey (VLSS) data. The extended r adio emission is n orthwest of the X-ray peak and extends southwest to 
northeast ~ 165 hz\ kpc at the resolution of the VLSS dCohen et al 1 l2007l) . 



C.43. RXJ1 504.1 -0248/RBS 1460 

RXJ1504 is the cluster with the highe st redshift and X-ray luminosity in HIFLUGCS, and shows the largest classical mass 
deposition rate. iBohringer et al.l (120051) reported the results to the Chandra observation of this cluster. This cluster was origi- 
nally not included in HIFLUGCS because its X-ray flux is only slightly (<20%) above the flux limit. RXJ1504 appears only 
margi nally extended in the ROSAT All-Sky Survey. Additionally the galaxy at the center of the X-ray emission is classified as 
AGN dMachalski & Condonl 1999b and its optical spectrum shows emission lines. It was assumed that even if there is only a small 
AGN contribution from the central AGN to the total X-ray flux (~20%), the cluster would fall below the flux limit. However, 
the C handra image reveals that there is actually no significant point source emission at the center of this cluster ( Bohringer et al.l 
2005), which argues against any significant contamination by AGN emission. The BCG features a compact and flat-spectrum 
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radio source, ( Mittal et al.ll2009 ). Therefore, this cluster is included into HIFLUGCS. This SCC cluster is one of sixteen clusters 
in which the BCG's central velocity dispersion is not available. 



C.44. A2029 

The spectra of the inner annuli fit best t o non-solar metalic ity ratios, but freeing ratios does not change the best-fit temperatures, 
so solar ratios were used for simplicity. IClarke et alj (120041) studied the core of this cluster in detail with Chandra. 



C.45. A2052 



32 h~\ kpc). iBlanton et al.l (120011) found 



The central region of this cluster requires a double thermal model out to ~45" ( 
prominent X-ray cavities in the original Chandra observation. They determined these cavities to be cospatial with radio lobes 
from the central radio source. 



C.46. MKW3S/WBL 564 



This SCC cluster shows some disruption in the core and bubbles to the south dMazzotta et al.ll2004l) . MKW3S is one of sixteen 
clusters in which data about the BCG's central velocity dispersion are not available. This cluster is a member of the Hercules 
Supercluster. 



C.47. A2065 

A2065 is a member of the Corona Borealis Supercluster, in projection close to the Hercules Supercluster but twice as distant. 



— 1/2 

This cluster is one of four clusters on the border between SCC and WCC clusters. Its CCT is (~ 1 .3 « 71 ' Gyr) is slightly longer 
than the 1 Gyr cutoff. This cluster shows an inwardly decreasing temperature profile as seen in the SCC clusters. This cluster 
is one of sixteen clusters in which the BCG's central velocity dispersion is not available. It is one of five clusters where the 
BCG peculiar velocity is more than 50% of clu ster velocity dispersion . This offset suggests possible sloshing which may have 
disrupted the CC. Based on the Chandra data, Chatzikos et al. ( 20061) suggest tha t the cluster is involved in an unequal mass 
merger and that one cool core has survived the merger. iFeretti & Giovanninil(ll994l) identified a WAT ~19' (1.6 hz^ Mpc) south 
south-west of the c luster center. The jet s of the WAT are bent away from the center of the cluster. In the NRAO VLA Sky Survey 
at 1.4 GHz (NVSS lCondon et al.|[l998l) . there appears to be a diffuse radio source -91" (-124 h~\ kpc) to the southwest of the 
cluster center. It is unclear whether this source is associated with the central radio source. 



C.48. A2063 

This WCC cluster appears to have a very regular morphology in X-rays, with some hint of an elongation to the northeast. The 
BCG resides at the X-ray peak. The NVSS shows three bright radio sources in a line along an axis from southwest to northeast 
but only the center source is associated with the BCG, while the other two are cospatial with two neighbouring galaxies. As with 

— 1/3 9 

many WCC clusters this cluster shows a flat central temperature profile and a raised central entropy Kq > 50 keV cm . This 
cluster is close to MKW3S. 



C.49. A2142 



This cluster has a double cold front ( Ma rkevitch et al .l2000h . The separation between the BCG and the X-ray peak is ~23 h n \ kpc 
for this cluster, making it one of fourteen with this value >12 h^ kpc. This is one of sixteen clusters in which no data is 
available for the BCG's cen tral velocity dispersion. It is possible this cluster hosts a radio halo, but the evidence remains dubious 
dGiovannini & FerettfeOOOt) . 



C.50. A2147 

A2147 is a member of the Hercules Supercluster. Due to the short observing time combined with the lack of a bright core, 
we used a 10" kernel when smoothing to determine the X-ray peak of this NCC cluster. There are three bright galaxies in a 
line near the core, of which the northernmost (the BCG) is located at the X-ray peak. This is one of the six NCC clusters in 
which M sp ec > 0. The X-ray morphology indicates that it is a merging cluster. The X-ray emission extends toward the south 
from the peak following the line of the three bright galaxies as well as extending to the southeast. There is, additionally, a sharp 
drop in the X-ray emission to the northwest. We argue that the observed M spec is not due to cooling, but results from multiple 
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temperatures along the line of sight caused by the m erger. Although the cluster has been labeled as a CC cluster in the past (e.g. 
Henriksen & Whitdl 19961) . [Sanderson et al.l (l2006al) found it to be an NCC cluster and likely merger system. 



C.51. A2163 



This well-known merging cluster contains the largest known radio halo dFeretti et alj |2001). The separation between the BCG 
and X-ray peak is ~ 158 \c\ kpc for this cluster, making it on e of eight clusters where t his va lue is >50 hz\. Our measurement of 
(~16keV) is higher than the value of ~12 keV found by Markevitch & Vikhlinin (2001) with data from the original, shorter 
Chandra observation. However, a recent measurement by Vikhlinin et al. ( 20081). using the same C handra as we, finds kJ vlr ~ 15 
keV, more consistent with our result. The difference between our result and I Vikhlinin et al" ( 2008 ) is barely inconsistent within 
1 a and is probably due to differences in the techniques used to determine kT v [ r in this extremely hot cluster. This is the second 
most distant and hottest cluster in the HIFLUGCS sample. This is one of sixteen clusters for which data on the BCG's central 
velocity dispersion are not available, however, since the BCG is not cospatial with the X-ray peak so this information is not 
important for our analysis. 



C.52. A2199 

The central region of this cluster requires a double thermal model out to ~29" (~17 h^\ kpc). 



C.53. A2204 



Recently iReiprich et al.l (|2008) determined the temperature of this cluster out to ~ r2oo using Suzaku. They find that the temper- 
ature declines all the way from 0.3 r2oo to r2oa, consistent with predictions of simulations. This is one of sixteen clusters where 
data on the BCG's central velocity dispersion is not available. 



C.54. A2244 



As with A1651. lDonahue et al.l(l2005l) claim it to be a radio quiet CC cluster, but lMittal et al. (2009) detect central radio emission 



with a bolometric luminosity of ~7 x 10 39 h n ^ ergs s 1 . Although this is not particu larly luminous, it is c onsistent with radio 
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activity in other CC clusters (Mit tal et al.H2009l) . Due to elevated entropy in the core. IDonahue et alJ (120051) concluded a major 
AGN outburst had disrupted the cooling flow. Like many WCC clusters, this cluster shows a flat temperature profile. However, 
we point out that the same was true of A1650 until a deeper observation revealed a slight temperature drop in the core. This is 
one of sixteen clusters in which the central velocity dispersion of the BCG is unavailable. 



C.55. A2256 



This well-known merging cluster is the only one of two NCC clusters that shows a systematic temperature decrease in the center. 
The temperature decrease is the largest of any NCC or WCC cluster. Surprisingly, the separation between the BCG and X-ray 
peak is ~ 132 hz\ kpc for this cluster, making it one of eight clusters where this value is > 50 hz\ kpc. Since this separation 
is quite large, the cool gas is not associated with the BCG. It is most likely this gas is the remnant of a C C (perhaps from a 
merging group) that has been s tripped its central galaxy. This cluster hosts both a radio halo and relic (e.g. IClarke & Ensslin 
2006tlBridle & Fomalontll 19761) 



C. 56. A2255 

Due to a short exposure time and lack of a bright core, we used a 10" kernel when smoothing the image to determine the X- 
ray peak. The separation between the BCG and X-ray peak is ~72 h^ kpc, making it one of eight clusters where this value is 
>50 /i^j 1 kpc. This cluster BCG also has by far the largest pecul iar velocity of any cluster; almost twice the velocity dispersion of 
the cluster. This cluster hosts both a radio halo and a relic (e.g. lFeretti et al 1ll997h . 



C.57. A3667 



This w ell-known merging cluster shows a very sharp cold front ( Vi khlinin et al kOQlblah and two radio relics (e.g Roettiger et al 



19991) . The separation of the BCG and X-ray peak is 

- WC 
Gyr. 



455 h-\ 



kpc for this clusters making it one of eight clusters where this 
value is >50 h~) kpc. This is the only WCC cluster with such a larger separation; however, A3667 is on the border between NCC 



and WCC clusters, with CCT 
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C.58. S1101/Sersic 159-03 



Kaastra et al.l Q2001) provided a detailed analysis of the XMM-Newton RGS and EPIC data. Recently I Werner et alj d2007l) have 



claimed discovery of a diffuse soft excess seen by XMM-Newton and Suzaku and suggest it is of non-thermal origin. This is one 
of sixteen clusters in which the central velocity dispersion of the BCG is unavailable. 



C.59. A2589 



This WCC shows a systematic temperature drop to wards the center, a lbeit r ather fiat (IcTo/kT = 0.93 and F=-0.079). Like A1650, 
it is on the cusp between SCC and WCC clusters. IZappacosta et al.l (120061) studied this cluster with a radio-quiet BCG in detail 
with XMM-Newton. They find the cluster to be exceptionally relaxed with a gravitating matter profile that fits a NFW profile with 
c v i r = 6.1 ± 0.3 and M v j r = 3.3 +0.3 x 10 14 M Q (r v j r - 1.74 ± 0.05 Mpc). They conclude that processes during halo formation act 



against adiabatic contraction. Additionally iBuote & Lewisl (120041) studied the original short Chandra observation that suffered 
from flaring. 

Following a method to determine the residual CXB (similar to what is described in Sect 12.3b . we measured the surface 
brightness profile out to 750 h^ kpc (~0.5 r v j r ). We fit this surface brightness profile to a double-/? model and the temperature 
profile to a broken powerlaw. The slope of the inner kT profile was fixed at zero and the outer kT profile fit well to a powerlaw 
of slope -0.36 with a break radius of 4'.2 (204 h^ 1 kpc). Using the fit to the temperature profile and double y8-model, we find a 
viriaj 24 ! mass and radius of M v j, = 2.7+0.8 x 10 14 hz) M Q and r v i r = 1.6 +0.2 fu. 1 Mpc respectively, consistent with the results of 



Zappacosta et al.l Q2006) 



C.60. A2597 



McN amara et al.l d2001l) analyzed the original, short, flared observation of A2597, noting the ghost bubbles. iMorris & Fabian 
d2005l) found high spectral mass deposition rates from the XMM-Newton EPIC and RGS consistent with ~100 M yr 1 down 
to almost keV, although the improvement to the spectral fits of the RGS data from the addition of a cooling flow model is 
marginal. The long Chandra ACIS observation shows a mass deposition rate of ~150 M yr~' down to ~1.3 keV and dropping 
to ~10 M Q yr 1 down to ~0 keV. 



C.61. A2634 

This WCC cluster contains the WAT source 3C465. There are a pair of galaxies (NGC7720) located at the X-ray peak. An 
extended bright X-ray halo (radius = ~ 12"), much brighter than the ICM emission, is cospatial with the galaxy pair. The halo 
seems to be associated with the larger southern galaxy. A2634 is the only CC cluster in the sample with M spe c > Mdassicai- The 
temperature profile shows a sudden drop at ~2.'7 (~ 100 \c\ kpc). Other than NGC7720, there is no obvious core in A2634 and 
the elongation of the ICM to the southwest is consistent with a merging cluster. Moreover the inverted temperature profile is 
more common in NCC clusters than in WCC clusters. We interpret the short cooling time and low M c iassicai as a cool core that 
has either been disrupted or is in the process of being destroyed by a merger. M spec may reflect the original mass deposition rate, 
but probably is strongly affected by multitemperature components along the line of sight in a merging system. 



C.62. A2657 

This WCC cluster has a slight increase in temperature in the central region. The Chandra image shows a cluster similar to e.g. 
A1650 and A2244. The central emission peak is clearly visibly but is not as sharply peaked as in SCC clusters. The overall 
ICM appears to be quite round, with some sloshing features (differently centered X-ray isophotes at different radii) in the central 
region. 



C.63. A4038 

The distance between the BCG and X-ray peak is ~12.4 tCj\ kpc for this cluster making it one of fourteen cluster s where this 
separ ation is >12 \c\ kpc. This cluster hosts a radio relic, close to but not connected to the central radio galaxy dSlee & Roy 



1998). 



C.64. A4059 

The central region of this cluster requires a double thermal model out to ~22" (~20 hz\ kpc). 



In this case r vil is defined for an overdensity of 104.7, as used bv lZappacosta et all J2006h . 
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Table 1. Observational Parameters (1) Cluster name, (2) Right Ascension and (3) Declination of the X-ray emission peak, (4) Number of 
observations used, (5) The total good time, (6) The fitted cluster virial temperature, (7) The fitted clustre metalicity and (8) The size of the 
region with T < T V1I (see section- l2~3l l. 
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14 h 40 m 43 s .08 
15 h 04 m 07;52 
15 h 10 m 56^06 
15 h 16 m 43 s .51 
15 h 21 m 51 s .75 
15 h 22 m 29532 
15 h 23 m 05*ll 
15 h 58 m 20:65 
16 h 02 m 16 s :78 



-09°18'11'.'0 
-01°15'2076 
-21 o 52'56'/0 
+33°15'21'/5 
+36°09'1479 
+06°01'37'.'4 
+ 13°01'52'/8 
+ 13°34'3978 
-44° 14' 18'.' 3 
-35°27'0079 
+09°58'01'/9 
+ 15°23'37'.'l 
-53°37'55'.'5 
+10°27'53'/8 
+02°24'36'.'2 
-08°33'37'.'9 
-61°27'11'.'0 
-13°15'3979 
-39°56'48'.'5 
-53°41'37'.'0 
-54°32'33'.'6 
+55°45'50'.'6 
-09°41'20'.'8 
-12°05'45'.'0 
-27°31'42'.'0 
+19°40'3072 
+01°53'46'.'l 
+03°39'18'/4 
+02°41'12'.'6 
-41°18'43'.'8 
-17 24'32'/3 
-01°45'42'.'7 
-04°11'49'.'2 
+27°57'3479 
-16 o 23'0678 
-27 o 10'26'.'8 
-31°29'43'.'2 
-31°40'17'.'0 
-32°51'57'.'5 
+26°35'33'.'l 

-27°oi'io75 

+03°27'57'.'7 
-02°48'16'.'8 
+05°44'41'.'4 
+07°01'1978 
+07°42'28'.'7 
+27°42'22'/2 
+08°36'26'/9 
+27°13'4972 
+ 15°58'25'.'6 



9 
1 
3 
1 
1 
1 
1 
2 
5 

12 
1 
1 
2 
8 
4 
1 
1 

2 
2 
1 
1 
1 
1 

2 
1 
2 
1 
1 

2 
7 
2 
7 
1 
1 
1 
1 
1 
1 
1 

12 
1 
1 
2 
3 
2 
1 
1 
2 
1 
1 



114.7 
11.4 
120.5 
43.5 
28.0 
21.5 
49.0 
29.6 
90.5 
429.7 
19.7 
23.3 
55.7 
101.4 
107.6 
9.8 
29.5 
66.4 
64.2 
17.1 
21.0 
27.0 
44.1 
186.7 
31.4 
68.7 
30.1 
11.9 
150.3 
288.6 
70.3 
225.3 
9.3 
29.8 
15.8 
14.8 
9.9 
19.4 
25.9 
173.5 
6.2 
23.1 
52.5 
107.0 
163.9 
56.7 
49.5 
28.5 
44.8 
17.0 



m 
-0.10 

^ 77+0.31 

4 82+^ 

4 09+ 08 
^■ u ^-0.l0 

11 13 +0 ' 39 

11.1J_0.43 
3 4 c+0.08 

3 16 +004 

t co +0.06 

2.01™* 

6.27+8:s 

9Q+^ 
-QQ2 
3.92 + °!»2 

5.09 

5.81 

6.34 

9.15 

1.22 

3.12 



+0.0 



M 

-0.07 
1 +0.27 
-0.27 
;+0.17 
-0.16 
>+0.03 
--0.04 
j +0.11 
--Q-12 
4QC+0.13 
H -^ J -0.15 
443+0.21 
H - HJ -0.16 

7 00 +013 

-Q 12 

6.08+™, 



w 
m 



1.97 
3.00 

9.53+ 1 
8.26 +a 



3 3 c +0.02 
-0.02 

QA+0.09 

J - 7 -0.09 

5 40 +0 - 20 

7 77+O.O6 
-'■ / -0.06 

8-40!] 

4____iS 



U - J5 -Q.Q5 

0-48+cuj 

o.4i!:? 

43+o.n 

"■^ -0.08 

60 +004 

u - uu -0.03 

0.50+ 007 
3Q + ^ 

a O4+0.06 

-QQ6 

a OO+0.05 
-R-R5 

27+™ 1 
"■^ ' -0.00 

0.54!; 
0.37^ 

041 +0.04 

-R-R3 
31+" 03 

-Q-Q3 
25+™ 1 

-RR1 

-Q.Q7 

26 + ™ 6 
-RR6 

u - uu -0.03 

31+0.04 

+0.13 
-0.09 

09 +01 ° 
59+ 007 

30 +0 05 
-RR6 

27 + <> 03 
-Q-Q3 

47+ 003 
' -0.03 

35 +0 04 

U-J -Q.Q2 

65 +m 
25 + ™ 9 

2O +0 02 

-RR2 

60 +m 
u.ou _ QQ1 

0.37 4 



0.33 



-0.04 

37+ 03 

u - J/ -0.02 

47+o.k) 

' -0.09 

33 +003 
-Q-Q3 

25 +0 04 

U-Z,J -0.02 

34+ 07 

-RR6 

a 97+O.O6 
' -Q.Q6 

44+0.07 

u - -Q.Q7 

0.45+!!", 



o.28!"-; 

53 + " : ' 
0.45 +B 
0.38_ 
0.384: 
0.66 
0.41; 
0.29; 
0.60; 
0.32; 
0.33; 



-0.02 

+0.05 
-0.04 
1+0.05 
-0.04 
1+0.04 
-0.04 
+0.24 
-0.23 
+0.06 
-006 



°^-0.8 


4 s+0.6 

-Q-5 

3 4 +0 - 3 
J - -0.2 

2 +01 




5 2 +08 
2 O^' 
5 5 +0 - 4 




5.2+ 07 

2.0!^ 

t 9+5:6 



4 8 +03 

^•°-0.3 







13.0!^ 



1 Q+0.4 
L8 -0.2 



2.0! 01 
3.16% 



5 6 +0 - 6 

2 2 +L2 





5 4 +L3 
J -^-o.i 



c 0+6.6 
3.8_ 23 





6 .9+0.6 

3 + ^ 
J - -0.3 



12.1+J 7 . 3 

3.7 
2.55 

6.2 



0\05 

+2.4 
-1.8 

5 1 +40 


>24.0 
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Table 1. continued. 



Cluster 


RA 


DEC 


#of 


GT 




(J2000) 


(J2000) 


Ob 


ks 


(1) 


(2) 


(3) 


(4) 


(5) 


A2163 


16 h 15 m 46 s :69 


-06 o 09'0073 


2 


80.3 


A2199 


16 h 28 m 38 s .32 


+39°33'01'/2 


1 


16.6 


A2204 


16 h 32 m 46:94 


+05°34'31'.'3 


2 


18.6 


A2244 


17 h 02 m 42*68 


+34°03'3973 


1 


56.0 


A2256 


17 h 03 m 14 s .26 


+78°38'59'.'9 


1 


8.0 


A2255 


17 h 12 m 34 s .15 


+64°04'11'.'5 


1 


39.6 


A3667 


20 h 12 m 42 s .66 


-56°50'48'/6 


7 


485.3 


S1101 


23 h 13 m 58*40 


-42°43'31'.'0 


1 


7.9 


A2589 


23 h 23 m 57 s :40 


+ 16°46'37'.'9 


3 


76.7 


A2597 


23 h 25 m 19 s .93 


-12°07'27'/5 


2 


112.0 


A2634 


23 h 38 m 29 s .25 


+27°01'54'/2 


1 


48.7 


A2657 


23 h 44 m 57 s :48 


+09°11'31'.'0 


1 


16.1 


A4038 


23 h 47 m 43 s .18 


-28°08'31'.'2 


2 


40.0 


A4059 


23 h 57 IT W93 


-34°45'33'/3 


2 


109.9 



kT vir 
keV 
(6) 



15 91 +tm 
1JJ -Oil 

4.37 + »S 



8.92 4 
5.78 4 



Z 
Solar 

(7) 
32 +U Uf> 
52+^ 
Ixf 
43+ 004 

-0.03 

26 +018 

n 97+0.06 
' -0.07 

35 +l) l)1 
-Q Qi 
n o "i +0.06 

-Q-Q5 

0.80!g 

0.38!g 
0.34+005 
n qs+&8o 

50 +m 

-003 



Core 

%r vir 
(8) 



i&fj 

-*— Q-ll 

7 61 +0 ' 65 

<r R1+0.19 
J -° -0.20 

6 39+ 04 

9 C7+O.I2 

-Q-13 

3 89+0-05 

4 05 +007 

o 1Q+0.11 

-Q-U 
3 52+0-12 

3 14+8:o3 

-003 






14.4^ 


L2 +0.7 



6.4 

1.' 
5. 



#5 

-3.1 
-1.2 
-0.6 
■1.0 
-0.6 








9 4+0.1 



' Abbreviated Cluster Names: 
2A0335: 2A 0335+096 
EXO0422: EXO 0422-086 
ZwC11215:ZwCl 1215.1+0400 
RXJ1504: RX J 1504. 1-0248 



' Alternative Cluster Names: 
Coma: A 1656 
Fornax: NGC 1399 
Centaurus: A3526 



Table 2. Derived Parameters: (1) cluster name, (2) central electron density, (3) central cooling time, (4) fractional central temperature drop, 
(5) slope of the central temperature profile (positive indicates a declining temperature with radius, see text), (6) cuspiness defined as jff^ at 
r = 0.04r5oo, (7) mass deposition rate within 0.048?5oo, estimated for a classical cooling flow (see text), (8) mass deposition rate estimated from 
the spectrum within 0.048f"5oo, (9) mass deposition rate down to a free lower temperature (7i w)> estimated from the spectrum within 0.048t"5oo, 
(10) the lower temperature found for M spcc2 , (11) central entropy and (12) cool-core classification - SCC = strong cool-core, WCC = weak 
cool-core and NCC = non-cool-core (based on the central cooling time, see sect, |3.2| l. 



Cluster 



(1) 



no 

10- 2 
(2) 



CCT 

h-\ 12 Gyr 
(3) 



6 92 +0 - 3i 

"' -0.52 

30 +0 08 
4 77+0.25 
^- ' -0.25 
9 43+0.56 



51+ 04 

14 03 +s 95 

0.48^ 
43 +0 (M 

"■^ -0.03 

8.04+;^ 
12.13+;^ 

8 - 81 -L08 
"""-0.05 

0.693" 

"'■"-001 

5-48!»;S 

8-22^ 
43 +010 

"•^ -0.07 
n 93+O.04 
u ""-0.03 

47+ 07 

"•^ ' -0.05 

7 62 +2M 

'•"^-1.60 

47+ 02 

"•^'-0.02 

16 47+ 3 10 
12 46 +249 

3 62 +a84 

j.uz, _ 059 

9 53 +237 

'"'-'-1.64 

0-41!^ 
2.871^ 

97 Q7+1.90 
-1.77 

0.281^ 
10.99 +2 £ 9 

o 2r° 6 ' 

"• xl -0.01 

0.42^J 
0-84^ 2 

| 95+O.29 
^ -0.20 

3 63 +043 

15 97+ 4 52 

1J " ' -2.95 

o 2r 004 

"• ZI -0.03 

16.59+i 3 ,, 7 



(4) 
34+0 0^ 

"■•'^-0.07 

90 +0 06 

"■'"-0.06 

41+ 004 

"•^ x -0.04 

65 +0 04 

o.uj _ 004 

35 +a01 
"•-'-'-o.oi 
n q-j+o.06 

"■'-'-0.06 

1 09 +006 
1 05 +a07 

56 +0 01 

u.ju _ 001 

0.52tS 
0.28!-! 

1 i i +0.04 
1- 11 -0.04 

1 15 +a05 

1,1 -0.05 

"■^ x -o.oi 
75 +a01 

"• /J -0.01 

89+ 03 

"•°'-0.03 

1 07 +a08 

I -"'-0.08 

21+0.01 
"•^-o.oi 

1 04 +0 - 04 

*■" M1.04 
i i r +0.08 
i - il -0.08 

1 oi +007 

0.91!-: 

0.86!^ 
76 +a02 

"■ ' -0.02 

99+0.03 

"■"-0.03 

90+ 003 

"■'"-0.03 

74+ 01 

"•'^-0.01 

1 06+ 008 

47+ - 03 

"•^'-0.03 
1QQ+0.002 
u - lyy -0M2 

0.39!-| 
0.84!- 3 

1 oi +007 

99 +0 07 
"•"-0.07 

60 +0 01 
"•""-o.oi 

1.06 + °°! 



(5) 



n 9QQ+0.022 

-\j.^jy_ 0Q22 

-0 103 +0 105 
-0 370 +0 039 

"••""-0.038 

-0 240 +0 039 
-0 540 +0 036 

"■ J -0037 
f) 1 9Q+0.088 
-0.092 

009 +0 091 
0.028 + °| 2 ? 
-0 226 +<l012 

""""-0.012 
f) 993+0.008 
-w.z,z.j„ 008 

-0 521 +H{)24 

" ,J -0.025 

-0 015+ 057 

"■" 1J -0.057 

0.044!- 77 
-0.2813;™ 

n 1 98 +0 006 

-0.147!»° 33 

0.218!°! 4 ° 
.n 347+0.008 

-0 090 +0 065 

"■"'"-0.065 

054+ 0102 
u.ujt |03 

032 +a087 

-0 092 +0 070 

"•"'^-0.070 

054+ 0120 
-0 079+"" 1 ' 4 

"•" ' '-0.004 

-0 027 +0 045 

"■" z - / -0.046 

-0 030 +0 065 

"■"-'"-0.065 

-o i6r 0018 

"• 1UI -0.018 

-0 007 +0118 

"•""'-0.124 

-0 1 47+0004 

U' 1 ^' -0.004 

-0 457 +0002 

"■ J -0.002 

-0 435 +0033 
-0.080!™ 29 
0.009!- 87 

-o 03r 0118 

"■"•"-0.123 

-0 194+ 0008 
"• ly -0.008 

0.094+- 74 



(6) 

no 7! 1+0-092 

590 +a27S 

"■-""-0.278 

1 255 +a013 

i,iJJ -0.O03 

0. 828!$$ 

i ris:q+0.004 

1, UJ -0.M3 
341+0.054 
" ,J1 -0.034 
fl 17Q+0.022 
"•-"'-0.024 

71 54+0.048 
U,JJ -0.004 

1 9 1 O+0.209 
■ °-0.209 

926 +0 - 257 

-0.257 

1 4,9+0.005 

-0.004 

0.481!- 2 J 

904+0.011 
" iZ -'^-0.011 

o.9io!« : ;° 9 
i.i48!° : | 7 ; 

1 97!7+0.015 
1 - AJI -0.015 

730 +0 016 

"■ ' -'"-0.023 

1 9 4f) +o.08i 
0.348!- 33 

fl 7!7Q+0.034 
" ,J ' '-0.040 

0. 8163*1 

462 +0 - 060 

1 9SQ+0056 

1, ZO '-0.056 
fl 8 1 7+0.067 
U - 81 '-0.067 

043 +0 003 
"■" H -'-o.ooi 

1 949+O.OO6 

-0.006 
n 97JQ+0.028 
"■^-"-0.027 
1 949+0-090 
L - AH -0.090 
| ion+0.021 

59 o+0.075 
"•-"'-0.075 
D Q95+0.149 
"•' z -'-O.I49 

0.642!- 29 
0.438!- 4 " 

1 7i>9+0.238 

-0.238 

0.105+-5? 



A^ckssical 

M Q /yr 
(7) 



*7f 



spec 

io/yr 
(8) 

<0.6 
8 9 +1 1 

°-'-1.2 

9 +0 3 

"•'-0.3 



M spec2 

M /yr 

(9) 

180.1!iJ 9 

<0.5 

35.8! 88 
9 4+0.7 
'• H -0.5 

q 9+0.8 

'•^-0.5 

2 4 +3 - 7 
^•^-1.0 

1 6 +5 2 
1 -"-i.o 

4 7+4.0 

H -'-1.5 

H0 5+107.8 
11 "--'_36.7 

1 96 +006 

1 -'"-0.04 

93.9 + 7 ° 4 3 

i +22 - 9 

4 0+141-4 

133.4 + ^- 5 
16.5!?f 

| 75 ri+10'91.0 

| 9 q+0.6 
lz -'-0.8 

54 6 +4 6 

-"•"_4.8 

<0.4 



kTi ow 
keV 
(10) 
236!*! 

o.o8 + r 

1 00 +0 13 

1 '""-0.15 

0-713$ 

0- 783$ 
0.893* 

0.08J.S- 82 

51 +4 06 

"• J1 -0.43 

1 S7!+0.38 
i.oj _ 038 

72 +0 02 
"■'^-o.oi 

65 +0 09 

"•" J -0.08 

3 10 +033 

-'• 1 "-0.40 

1 15 +305 

1 ' 1J -0.85 

0.293!' 
0.783" 

9 ^9+0.85 
z "-' -1.02 

0.08J 31 

| 97+O.O7 
1 ' z -'-0.07 

o.ii! 79 79 

0.08S 2 
0.08J. 79 - 82 

2 5 5 +77 - 35 

z "- JJ -2.46 
9 Q1+1 91 
z " yl -0.91 

1 20 +0U 

1 - z '"-0.17 
1 471+0.47 
1 -^-'-1.34 

0.583; 47 
1.093»? 

7 l+63.09 
"■ ' 1 -0.63 

0.42!- 2 
0.863" 

1 40 +0 13 

4-^"-0.14 
9 55+O.62 
i - J -'-0.43 

2 30 +L52 

z —'"-2.11 
9 9Q+138 
^•^'-0.94 

62 +0 01 
"■" -0.01 

0.58 + Z 9 ; 32 



,-1/3 



^0 

keV cm 2 

(11) 
12.0! 2 - 4 
246 9 +45 4 
12 4 +L3 

H.2!!| 
7.83J 

104.4!! 2 - 9 

267.3!| 9 - 7 

250.9! 2 ! 9 

14.0 + |-5 

13.6±« 

5 -8-03 

94.0!=^ 
190.8! 99 
13.8 + |' 7 

6 7 +0 i ' 

"•'-0.6 

18.0!}- 4 
260.2!« 2 

8.23 7 
222.0! 2 ; ; | 

z-/ j.z,_ 296 

211.1! 247 
86.2*1" 

276.7! 38 ; 8 
14 2 +a5 

1H ' -0.5 

61 5 +1L4 
296.5!| 32 
12 +0 - 6 

IZ -"-0.6 

246 4 +26 4 
z.tu.t_ 264 

6.53 2 

8-13^ 
19 2 +35 

1 '- z '-3.5 

45 7 +5 - 8 

HJ - '-5.8 

118.7 +97 
364.6! 398 

7 5+0?8 ' 
'• J -0.8 

228.0+"? 



ccc 



(12) 



A0085" 
A0119 
A0133" 
NGC0507" 
A0262" 
A0400 
A0399 
A040F 
A3112 
NGC1399" 
2A0335" 
IIIZw54 
A3158 
A0478 
NGC1550 
EXO0422 
A3266 
A0496" 
A3376" 
A3391 
A3395s 
A0576'" 5 
A0754 /i 
A0780 
A 1060 
A1367" 15 
MKW4 
ZwC11215 7 
NGC4636" 
A3526" 
A 1644" 
A1650 
A1651 
A1656 
NGC5044 
A1736" 



4.06 + 



;+0.28 
-0.28 

0.28!°- 

0. 45!- 4 

67 +a08 

"•"'-0.08 

8.20!|- 34 

1 qq+o.93 
1,7 -0.93 

7 03 +ai7 

'•"-'-0.17 

o 5r - 04 

"■ J1 _0.04 

52 +0 - 02 

"■ J -0.02 

io.o3!i L8 ! 

5 qq+0.87 
J -"-0.87 

5 48 +a59 

-'• HO -0.59 

7 6 +o.i8 

"•'"-0.18 
5 7q+0.27 
J -"-0.27 

24+ 03 

"■^-0.03 

35 +a05 

"■ J -'-0.05 

35 +0 05 

"• J -0.05 

90 +016 

"•'"-0.16 

64 +012 

"•" -0 17 
7 7q +0.29 
' • ' -0.29 

1 15 +a32 

1,1 -'-0.32 

o ir 000 

"• 11 -0.00 

4 36 +a3 ° 

t.ju Q30 

44 +005 

"• H ^-0.05 

2 64 +0 " 
4 5q+0.07 
^•-"-0.07 
t 37+O.9O 
-'-"-0.90 
71 4q+0.65 
-'■^'-0.65 

1 95+O.IO 

1, z "'-0.10 

40 +008 

"• H "-0.08 

3 04+°-47 
0.17+-Z 



131.8 



68.4!" 
7.7! 2 ; 2 

12 r 14 



146.3! 4 ! 2 

0.33 2 " 

267.8!i 2 ; 8 
18.13? 



559 3 +24L7 
9-311;! ' 
42.5! 82 

g 2+40.6 

78.7!fe 



12.9 + 



165.4 + ! 2 - 8 

6-63:8' 



6.831 



1 +0 1 
1 -"-o.i 

23 5 +0 5 

"•-'-0.5 

24 6 +1L9 

^■"-11.9 
qT 7+28.2 
' J - -28.2 

165.8i 2 - 8 



2.9 + 



3+0.2 
-0.2 
[ 1+0 6 

1 6 +1 5 
1 -"-i.o 

4 +2 7 

<1.6 

52 +0 02 

"•-' z '-0.02 

65.8! 37 
<0.2 
7 +a9 

"•'-0.4 

136.7!^| 

3 9+0.5 
J - -0.5 

<0.9 
12 9 +a7 

lz -'-0.8 

4 6 +06 

^■"-0.6 

<0.4 
<0.93 
<0.4 
<0.3 

1+2.1 
-1.2 

14 7 +1 ° 
iH -'-i.o 

0.6 +02 



2.0 + 



J -0.2 

4 +0 2 

"■^-0.2 

<0.1 



2.5 + 



11.0 



+2.5 
2.5 



;+3.2 

0.483| 
5-23! 

9 3+0.7 

<0.7 

7 7+7.0 
'•'-4.8 

<0.3 
0.83; 3 
0.3+25 



<2.4 
<3.1 



20.8 



+ 135.8 
-13.0 

42 1 +17 - 9 

^ z " 1 -13.9 
74 7+14.5 
'^■'-14.7 
4 4+4.0 

6 +0 5 
"•"-0.2 
q 3+19 
'• J -1.5 

6.0 +259 ° 
1.5 + 



5+0.2 
J -0.2 



+0.3 
■0.3 



13.9_ u 

97 7+2.O 
^'•'-1.9 

84.5! 28 ;^ 
69.0!^' 

11 1 +48 
ii.i 38 

27 9 +L1 

^■'•'-0.6 

0.5+S 1 ,- 7 



SCC 
NCC 
SCC 
SCC 
SCC 
NCC f 
NCC 
NCC 
SCC 
SCC 
SCC 
WCC 
NCC 
SCC 
SCC 
SCC 
WCC 
SCC 
NCC 
NCC 
NCC 
WCC* 
NCC 
SCC 
WCC 
NCC 
SCC 
NCC 
SCC 
SCC 
SCC 
WCC* 
WCC 
NCC 
SCC 
NCC 



Table 2. continued. 



Cluster 



(1) 



no 

*4f 10- 2 cm- 3 
(2) 

c+0.97 



CCT 

7i' /2 °y r 

(3) 
1 69 +1 15 
5 15 + ° 72 

J - -0.57 
9 13+0.71 

61 +002 
u.ui _ 002 

55 + ° 06 

"• J -0.05 

10.87!^ 
59 +002 

"■ J7 -0.02 

53 + °° 4 

"• J -0.04 

0.511Z 

0. 86^11 

1 34 + ° 37 

• J -0.24 

2 36 + ° 14 

I 94+O.I6 

1, ?H -0.14 

17.04+ 3 ^ 
9 65 + ° 73 
60 +0 07 

"•""-0.06 

25 + ° 01 
"•■"-o.oi 

1 53 +0 - 27 
11.56+ 24 ? 
20.66^11 
6.14!°- 
0.881°;?° 
1.18!^ 
42 +0 04 

"■^-0.03 

2-681^ 
1.681°;]? 
o.70 + !!!!I 



(4) (5) 

ToF^ -o.083!«r, 

0.95+°-- -0.1 



(6) 

0- 738!";- 

0. 675!^ 

703 +O157 

"•'"-'-0.157 

1 ioc+0.028 

002g 

1 298+ - 013 

1. Z7O_ 0013 

"• JI7 -0.048 
i 171+O.OO6 

1- 1 ' -0.006 
1 i i 9+O.O66 
1-1 ii -0.066 

1 205 +a ° 29 

1 138+0.186 

469+o.oso 

n Q3 5 +0.069 
\j.yjj _ .069 

893+ 0006 

f) 993+0.035 

25 r°- 007 

1 027+ 0101 

1 -0.101 
1 7!Q4+0.008 
•■""-0.007 
fl 548+ -029 
"""-0.029 

234+ - 031 

u ' -0.032 
n 9Q5+0.068 
„ 56 

n 994+O.OO5 
-O.0O5 

66 r°- 188 

"•" ul -0.188 

569+ 0164 

1 QCn+OOK-l 

l.J8U_ 00g4 
240 +0 034 
907+ 022 

"''"'-0.028 

966+ 007 

"•'""-0.007 

639+ 018 



^classical 

M /yr 
(7) 



-15.3 

23 6 +1(w 

"•"-10.9 
43 ,+24.4 

255 4+ 15 -° 

+5.8 
5.8 



M Q /yr 
(8) 

tTXT— 

-1.7 

1 7 +1 1 
1, '-1.0 

6 2 +1 3 
14 5 +u 
1 2 +1 ° 
<0.2 
<7.4 
18.7^ 

i3.o:;: 

9 Q+10 
^•'-1.0 
9 7+1.9 
A ''-1.7 

1.3!- 
25-01;? 

3 i +0.4 
-'■-0.8 

57 3 +27 -' 

' '-^-23 4 
9 3+1.0 
i "'-1.0 

<17.5 
1.7-? 

8-5-^ 
4.0!" 

<o.i 

14 1 +4 1 



M spec2 
M /yr 
(9) 

9 7+2.0 

z "'-1.7 
17 9 +3847031.1 
1 '-^-10.3 
x O+14.0 

222 4 +32 - 8 
21 5 +91 

—5.1 

5 +45 4 

"• J -0.3 

2671.5_ l278 5 
328.9!?|- 7 
33.9!| 8 
40.0!| 7 - 4 
47.2! 8 ? 

5-3! 2 » 
184.1!^ 

3 l + ° 5 

J - -0.8 
87 3+73.0 

/ . J_37 g 

76.9! 4 j 
913 8 +1 191 

us.?": 

8.5! 2 ; 3 
4.0!°;= 

23.5 + 7 6 4 5 
9 6 +3 4 

'•"-3.5 

141. 8 +3 ^ 

1 +0 3 

1 3 n+ 196.8 
1 -'-'-8.4 

<0.2 
38.6! 4 ;? 



kTi ow 
keV 

(10) 

12 +79 - 78 

"' -0.04 

1 81+2 21 
L81 -0.62 
9 J9+0.50 

-0.87 

i.80!°;= 

0.88!°;= 
1.48!] 8 4 42 
3 26 + ° 82 

-'• z "-1.15 
i+0.27 
-0.15 

°- 85 0.20 

1 55 +0 - 45 
1 74+0.52 
1 -'^-0.34 

1 07 + ° 65 

1 '"'-0.50 
9 91+0.59 
■ i " il -0.47 

0-08 + ° 00 78 
1.43!?;° 7 

1 co+0.14 
J -0.08 

1 89 + ° 39 

1 '°'-0.26 
9 43+0.85 

0- 08!°; 4 

o.oso+r 

9 3 i +0.65 

1 - -' 1 -0.54 

61 + ° 09 

"•"-0.08 

1 33+O.28 

1 97+0.17 
1 - i '-0.16 

0. 08 + ° 48 

i OT + 1.74 

1. »J_ 07] 

0.08J.S- 82 

| 34+0.09 



K 

h 71 1/3 keV cm 2 
(11) 
59 4 +15 - 9 

Jy - -15.9 

113.7+ 8 - 8 



ccc 



(12) 



A3558 
A3562'' 
A3571 
A1795 
A3581" 
MKW8 
RXJ1504 
A2029 
A2052" 
MKW3S 
A2065 15 
A2063 
A21427 
A2147 
A2163/ 3 
A2199" 
A2204 15 
A2244 
A2256" 
A2255' 5 ' 5 
A3667^ 
S1101 
A2589 
A2597 
A2634 
A2657 
A4038 5, 
A4059 mS 



2.46! 097 

0.71!°;° 7 
9 4Q+O.6O 
^■^'-0.60 

6 ir 019 

"• 11 -0.19 

2 63 +0 - 23 

—0.23 

29 +0 ° 5 

"•^'-0.05 

7 Q9+0.19 
'■'^-0.19 

9 02 +0 - 64 

'■"' 6 -0.64 

4 62 +018 

H -"^-0.18 

3 46+ - 58 

J '^"-0.58 
9 QQ+0.61 
^•"-0.61 

1 9Q+0.06 

2 65 +0 - 17 

—0.17 

0.23!°;° 3 
0.84!°;° 3 

5 61 +0 - 58 

J - ul -0.58 

16 70 +0 - 65 

lu -'"-0.65 
9 q i +0.40 
71 -0.40 

0.38!°;°« 
0.26!- 7 
0-78!^ 

9 Q-7+0.49 
2 -°'-0.49 
9 49+0.52 
^■•^-0.52 
7 74+0.57 

9 9 i +0.52 
z " z ' 1 -0.52 
1 9|+0.36 
i ' zl -0.36 
1 7Q+0.11 

■^•"-o.n 

3 85+0-35 



J -0.05 

0.92+ 005 



--0.05 

97+0 02 

"•"-0.02 

23 +0 - 01 



-0.05 

95+ 05 -0 

"■ 7J -0.05 " 

50 +0 01 -0 

"• J "-o.oi " 

52 +0 ° 3 -0 

"• JZ -0.03 " 

1 20 +0 ° 6 

1 -^"-0.06 "' 

0.48!°° 4 -0 

43+° :o1 -0 

"■ HJ -o.oi " 

0.304!°;°- -0. 

0.78!°-°5 -0. 

79+ 03 -0 

"•"-0.03 " 

97+0 03 _q 

"•"-0.03 " 

0.70 + °°5 -0. 



-0.05 

11+0.14 q 

0.31!°;°^ -0. 

0.38!°,°, -0. 

93+0 04 _q 

"■ 7 - 3 -0.04 " 

0.63!°;^ -0. 

1.18!°° 9 -0 

0.83!°° 4 -0. 

0.88!°g 4 -0. 

0.93!°°? -0. 

0.57!°;°! -0. 

1.00!°° 7 0. 

1.12^2 o. 



-0.067 
HS7+0.073 
■ U82 -0.073 
1 nq+0.090 
lu7 -0.093 
91-13+O.OO8 
'^"-'-0.008 
3C3+O.O54 
JJJ -0.055 
162 +0.074 
lu -0.076 
91 3+O.OI8 
,z ' i -'-0.018 
9/-19+O.OO9 
-0.010 

QC4+0.018 

■ 7JH -0.018 

Q95+0.015 

1 93+O.O46 
1Z-J -0.047 
Q, 6 +0.047 
■" lu -0.047 

, 149+0.019 

• x ^ 7 -0.019 

l.058!°;° 8 ? 

9CA+0.202 
■ 28U -0.219 
61 Q+0.093 
,ul7 -0.097 
3CQ+0.034 
,JJ -0.033 
A3Q+0.069 
, "-' o -0.071 
943+O.O6O 
,z '^-'-0.061 

l.088!°;!? 9 

1 33+O.O8O 
IJ -'-0.083 
1 no+0.034 
1U8 -0.034 
n7q +0.042 
■"' 7 -0.042 
181 +0 010 
loI -0.010 
006 +0.109 
"""-0.112 
02 3+0.057 
■"^ -0.057 

I 024 + ° 039 

■"^-0.039 
I co 1 +0.028 
UO '-0.028 



31.5+ 



33.5 



2.7+ 



2266.5!|? 97 
324.2! 37 ; 8 
69.6! 38 

79 1 +24.2 
' z - I -24.2 

50.6 + ?|- 3 
28-6!l! 

97c 7+40.4 
i/J -'-40.4 



78.9 + 



H2.0 
-12.0 

637.8! 448 
131.8^ 7 



7 + ° 3 
"•'-0.3 
999 3+120.0 

z,z,z..j_ 1200 

1Q Q+12.3 

ly - y -l23 
774 9+36.4 
Z '^' Z '-36.4 

2 + ° 1 

-0.1 

15 4 +10 - 9 
1J - -10.9 

47 2 +6 3 

^'•^-6.3 

no.6 



0.8 + 



57.7^ 



10.6 



2+0.3 
-0.3 

10.0! 2 ;« 

<0.4 
<0.2 
4 l + ° 5 



i.e. 



77 Q+13-3 
' '-^-13.3 

20.7 + ° 7 

11 6 +i: ° 

177 3 +2l ° 
1 ' '•-'-21.0 

74 9+0.9 

24.9!| 4 

o 7+0.5 

28.6! 3 ; 2 

44 3 +62 
^•-'-6.2 

66. r? 8 

65.8!t° 
212.4!? 3 ;! 
427.6!^ 8 

Q 9+2.0 
y - z --2.0 

11 r°- 5 

57 1+ 5 - 8 

197.6!g- 4 
365.0! 7 °;| 
134.6+ 8 - 7 
24 2 +2 - 9 

-2.9 

43 4 +6J 
^-'•^-6.3 
1 3 9+0.8 

-0.8 

40 7 +6 - 9 
^"- ' -6.9 

75 l +l5 ° 

44 5+2.1 

8.4 + ° 7 



wcc 
wcc 
wcc 
sec 
sec 

NCC 

sec 
sec 
sec 
sec 
wcc* 
wcc 
wcc* 

NCC 
NCC 

sec 
sec 
wcc 

NCC* 
NCC 

wcc* 
sec 
wcc* 
sec 
wcc 
wcc 
wcc 
sec 



One or more of the innermost annuli were fit with a double thermal model which significantly (> 99% confidence according to the F test) improved the fit over a single 
thermal model. 

13 This is one of eight clusters in which the offset between the X-ray peak and BCG is quite large (> 50 /i 71 ' kpc). 

7 This is one of six clusters in which the offset between the X-ray peak and BCG is moderate (> 12 /! 71 ' kpc but < 50 /7 71 ' kpc). 

* The BCG of this cluster is one of five that has a peculiar velocity >50% of the cluster velocity dispersion. 

T Although this cluster is an NCC there is an appreciable amount of cool gas around the emission peak so that the temperature profile is not consistent with a non-negative 
slope. 

* This WCC cluster has an appreciable amount of cool gas around the emission peak so that the temperature profile is not consistent with a non-negative slope. 



